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Plate 36. 

The following paper contains the investigation of the mass-motion of viscous and 
imperfectly elastic spheroids, as modified by a relative motion of their parts, produced 
in them by the attraction of external disturbing bodies ; it must be regarded as 
the continuation of my previous paper/" where the theory of the bodily tides of such 
spheroids was given. 

The problem is one of theoretical dynamics, but the subject is so large and complex, 
that I thought it best, in the first instance, to guide the direction of the speculation 
by considerations of applicability to the case of the earth, as disturbed by the sun 
and moon. 

In order to avoid an incessant use of the conditional mood, I speak simply of the 
earth, sun, and moon ; the first being taken as the type of the rotating body, and the 
two latter as types of the disturbing or tide-raising bodies. This course will be justi- 
fied, if these ideas should lead (as I believe they will) to important conclusions with 
respect to the history of the evolution of the solar system. This plan was the more 
necessary, because it seemed to me impossible to attain a full comprehension of the 
physical meaning of the long and complex formulas which occur, without having 
recourse to numerical values ; moreover, the differential equations to be integrated were 
so complex, that a laborious treatment, partly by analysis and partly by numerical 
quadratures, was the only method that I was able to devise. Accordingly, the earth, 
sun, and moon form the system from which the requisite numerical data are taken. 

It will of course be understood that I do not conceive the earth to be really a 
homogeneous viscous or elastico-viscous spheroid, but it does seem probable that the 
earth still possesses some plasticity, and if at one time it was a molten mass (which is 
highly probable), then it seems certain that some changes in the configuration of the 
three bodies must have taken place, closely analogous to those hereafter determined. 
And even if the earth has always been quite rigid, the greater part of the same effects 
would result from oceanic tidal friction, although probably they would have taken 
place with less rapidity. 

# " On the Bodily Tides of Viscous and Semi -elastic Spheroids," &c, Phil. Trans. 1879, Part I, 
MPCCCLXXIX. 3 M 
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As some persons may wish to obtain a general idea of the drift of the inquiry with- 
out reading a long mathematical argument, I have adhered to the plan adopted in my 
former paper, of giving at the end (in Part III.) a general view of the whole subject, 
with references back to such parts as it did not seem desirable to reproduce. In order 
not to interrupt the mathematical argument in the body of the paper, the discussion of 
the physical significance of the several results is given along with the summary ; such 
discussions will moreover be far more satisfactory when thrown into a continuous 
form than when scattered in isolated paragraphs throughout the paper. I have tried, 
however, to prevent the mathematical part from being too bald of comments, and to 
place the reader in a position to comprehend the general line of investigation. 

Before entering on analysis, it is necessary to give an explanation of how this 
inquiry joins itself on to that of my previous paper. 

In that paper it was shown that, if the influence of the disturbing body be expressed 
in the form of a potential, and if that potential be expressed as a series of solid 
harmonic functions of points within the disturbed spheroid, each multiplied by a simple 
time harmonic, then each such harmonic term raises a tide in the disturbed spheroid, 
which is the same as though all the other terms were non-existent. This is true, 
whether the spheroid be fluid, elastic, viscous, or elastico- viscous. Further, the free 
surface of the spheroid, as tidally distorted by any term, is expressible by a surface 
harmonic of the same type as that of the generating term ; and where there is a 
frictional resistance to the tidal motion, the phase of the corresponding simple time 
harmonic is retarded. The height of each tide, and the retardation of phase (or the 
lag) are functions of the frequency of the tide, and of the constants expressive of the 
physical constitution of the spheroid. 

Each such term in the expression for the form of the tidally distorted spheroid may 
be conveniently referred to as a simple tide. 

Hence if we regard the whole tide-wave as a modification of the equilibrium tide- 
wave of a perfectly fluid spheroid, it may be said that the effect of the resistances to 
relative displacement is a disintegration of the whole wave into its constituent simple 
tides, each of which is reduced in height, and lags in time by its own special amount. 
In fact, the mathematical expansion in surface harmonics exactly corresponds to the 
physical breaking up of a single wave into a number of secondary waves. 

It was remarked in the previous paper,* that when the tide- wave lags the attraction 
of the external tide-generating body gives rise to forces on the spheroid which are not 
rigorously equilibrating. Now it was a part of the assumptions, under which the 
theory of viscous and elastico-viscous tides was formed, that the whole forces which 
act on the spheroid should be equilibrating ; but it was there stated that the couples 
arising from the non-equilibration of the attractions on the lagging tides were pro- 
portional to the square of the disturbing influence, and it was on this account that 
they were neglected in forming that theory of tides. The investigation of the effects 

* "Bodily Tides," &c. Sec, 6, 
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which they produce in modifying the relative motion of the parts of the spheroid, that 
is to say in distorting the spheroid, must be reserved for a future occasion/* 

The effects of these couples, in modifying the motion of the rotating spheroid as a 
whole, affords the subject, of the present paper. 

According to the ordinary theory, the tide-generating potential of the disturbing 
body is expressible as a series of Legendre's coefficients ; the term of the first order 
is non-existent, and the one of the second order has the type fcos 2 — \. Throughout 
this paper the potential is treated as though the term of the second order existed alone, 
but at the end it is shown that the term of the third order (of the type fcos 3 — fcos) 
will have an effect which is fairly negligeable compared with that of the first term. 

In order to apply the theory of elastic, viscous, and elastico-viscous tides, the first 
task is to express the tide-generating potential in the form of a series of solid harmonics 
relatively to axes fixed in the spheroid, each harmonic being multiplied by a simple 
time harmonic. 

Afterwards it will be necessary to express that the wave surface of the distorted 
spheroid is the disintegration into simple lagging tides of the equilibrium tide-wave of 
a perfectly fluid spheroid. 

The symbols expressive of the disintegration and lagging will be kept perfectly 
general, so that the theory will be applicable either to the assumptions of elasticity, 
viscosity, or elastico-viscosity, and probably to any other continuous law of resistance 
to relative displacement. It would not, however, be applicable to such a law as that 
which is supposed to govern the resistance to slipping of loose earth, nor to any Jaw 
which assumes that there is no relative displacement of the parts of the solid, until 
the stresses have reached a definite magnitude. 

After the form of the distorted spheroid has been found, the couples which arise 
from the attraction of the disturbing body on the wave surface will be found, and the 
rotation of the spheroid and the reaction on the disturbing body will be considered. 

This preliminary explanation will, I think, make sufficiently clear the objects of the 
rather long introductory investigations which are necessary. 

PART I. 

§ 1. The tide-generating potential. 

The disturbing body, or moon, is supposed to move in a circular orbit, with a 
uniform angular velocity — /2. The plane of the orbit is that of the ecliptic ; for the 
investigation is sufficiently involved without complicating it by giving the true 
inclined eccentric orbit, with revolving nodes. [I hope however in a future paper to 
consider the secular changes in the inclination and eccentricity of the orbit and the 
modifications to be made in the results of the present investigation.] 

* See tlie next paper " On Problems connected with the Tides of a Viscous Spheroid." Part I. 

3 M 2 
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Let m be the moon/s mass, c her distance, and r=-f^. 

c 6 

Let XYZ (Plate 36, fig. 1) be rectangular axes fixed in space, XY being the ecliptic. 
Let M be the moon in her orbit moving from Y towards X, with an angular 
velocity fl. 

Let ABC be rectangular axes fixed in the earth, AB being the equator. 

Let % t// be the coordinates of the pole C referred to XYZ, so that i is the obliquity 

of the ecliptic, and -~ the precession of the equinoxes. 

Let ?, 9, (j> be the polar coordinates of any point P in the earth referred to ABC, 
as indicated in the figure. 

Let ot* l9 <a 2 , o) 3 be the component angular velocities of the earth about the instan- 
taneous positions of ABC. 

Then we have, as usual, the geometrical equations, 

u'l 

-j- sini=-~ o) x cos ^;— cousin x f* ...... (1) 

dy . chfr 

^~ir+ir cos % = w q 

dt dt ° j 



Let n coseci be the precession of the equinoxes, or -~, so that -£ = n cot i— g> 3 .* 

Now the earth rotates with a negative angular velocity, that is from B to A ; therefore 

if we put -j"=n, n is equal to the true angular velocity of the earth +U cot i But for 

purposes of numerical calculation n may be taken as the earth's angular velocity ; and 
care need merely be taken that inequalities of very long period are not mistaken for 
secular changes. 

Let the epoch be taken as the time when the colure ZC was in the plane of ZX, 
when x was zero an( l the moon on "the equator at Y. It will be convenient also to assume 
later that there was also an eclipse at the same instant. A number of troublesome 
symbols are thus got rid of, whilst the generality of the solution is unaffected. 

tjf ITT" 

Then by the previous definitions we have x =nt ? MN=/2£, NR=- — RD=- — (<£— x)* 

Now if w be the mass of the homogeneous earth per unit volume, then the tide- 
generating gravitation potential V of the moon, estimated per unit volume, at the 
point r, 0, <j> or P in the earth is, by the well-known formula, Y=wr? s (cos <i PM— |-). 

This is the function on which the tides depend, and as above explained, it must be 

* The limit of n cot i is still small when i is zero. In considering the precession with one disturbing 
body only, TI cosee^ is merely the precession due to that body; but afterwards when the effect of the sun 
is added it must be taken as the full precession, 
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expanded in a series of solid harmonics of r, 0, 0, each multiplied by a simple time 
harmonic, which will involve n and fl. 

For brevity of notation nt 9 fit are written simply n 9 fl, but wherever these symbols 
occur in the argument of a trigonometrical term they must be understood to be multi- 
plied by t the time. 
We nave 

cos PM= sin cos MR+ cos sin MR sin MRQ 
and 

cos MR= cos MN cos NR+ sin MN sin NR cos i 

= cos fl sin (0— n)-\- sin fl cos (0— n) cos i 
also 

sin MR sin MRQ= sin MQ= sin fl sin i 
Therefore 



cos PM= sin sin (0— n) cos /2 + sin cos ($—n) sin /2 cos i+ cos sin /2 sin i 
=-J- siii #{sin[0— tyi— fl)~\-\- sin[0— (n-\-fl)]} 

+isin #cosi{sin[0— (^— /2)] — sin[0 — (n+/2)]}+ cos #sin/2 sin i 

Let 



^ 



_p=. cos -, q= sin 



2 



„&. ilUll 



cos PM=p 3 sin sin[0— (^— /2)]+2j3f/ c °s sin /2+g 3 sin sin[0— («+/2)] . (2) 



Therefore 



cos 3 PM=^ 4 sin 3 0{ 1 — cos [20-2(n-/2)]} + 2pV cos 3 0(1 — cos 2/2) 
+^sin 3 0{l — cos [20— 2(ft+/2)]} + 2p 3 ^sin cos 0{eos(0— /&)— cos[0 — (n— 2/2)]} 
+ 2pg 3 sin 0cos {cos [0— (ft +2/2)] — eos (0— n) } +jp 3 g 3 sin 3 0{cos 2/2— cos (20— 2n)} 



Then collecting terms, and noticing that 



£(P*+2*) sin " ^ + 2pV cosS ^=i+i(l-6j9 3 c/)(|- cos 3 0) 



we have 



V 1 

_— \j\JO X -LtJ_ n 



— |sin 3 0{p 4 cos[20— 2(n— /2)]+2jo 3 vf cos[20— 2ft] +g 4, cos [20— 2(ft+/2)]} 

— 2sin0cos0|j9 3 gcos[0— (ft— 2/2)]— pq(p^— g 3 )cos(0— ft) — _pg 3 cos[0— (ft+2/2)]} 

+ (i-co S ^){3pYcos2i2+l(l-6pV)} (3) 



Now if all the cosines involving be expanded, it is clear Miat we have V consisting 



452 MR. G. H. DARWIN ON THE PRECESSION OF A VISCOUS SPHEROID, 

of thirteen terms "which, have the desired form, and a fourteenth which is independent 
of the time. 

It will now be convenient to introduce some auxiliary functions, which may be 
defined thus, 

<p(2n)=±p 4 > cos 2{n— n)+2 j2( l* cos 2n + W cos 2(n+^2) 

ty(n)=2p*q cos (n — 2D) — 2pq(p % —q % ) cosn— 2pg 3 cos (n-\-2fl) y . . . (4) 

X(2Q)=3pY cos2 ^ 

<E> (2n — -|tt) , ^(w— |7r), X(2Q — \tt) are functions of the same form with sines replacing 
cosines. When the arguments of the functions are simply 2n, n, 2/2 respectively, they 

will be omitted and the functions written simply <&, M*, X; and when the arguments 
are simply 2n— \ir. n—^ir, 2fl~^ir, they will be omitted and the functions written <£>', 
M>', X'. These functions may of course be expanded like sines and cosines, e.g., 
^(n — ot) = ^f cos a+M*' sin a and ^{n — a)=*!r' cos a— M> sin a. 

If now these functions are introduced into the expression for V, and if we replace 
the direction cosines sin 6 cos </>, sin 6 sin <j>, cos # of the point P by £, tj, £, we have 

V - 3 =-(f - 7 f)<s>-2£ n &-aV- v &'+Ue +^-2^)[X+Ml-6^ 2 )] . (5) 



wrr 



P— t^ 2 , 2^, ££, 17 £, i(P+V 3 "" 2 £ 2 ) are surface harmonics of the second order, and 
the auxiliary functions involve only simple harmonic functions of the time. Hence we 
have obtained V in the desired form. 

We shall require later certain functions of the direction cosines of the moon referred 
to A B C expressed in terms of the auxiliary functions. The formation of these 
functions may be most conveniently done before proceeding further. 

Let x, y, z be these direction cosines, then 



whence 



cos PM = x£-\- yrj + %t> 

(K)B a PM-i==(a;f+y^+«0 s -t(?+^+O 

= ^(^-*)+^V-*) + W-*) + 2^+2£^+2^^ . (6) 

But from (5) we have on rearranging the terms, 

cos a PM-*= 

P{-*+lX+i(l-6pV)}+7 ? 3 {<I>+lX+l(l-6pY)}4-r{-|X-i(l-6 i >V)} 
-2 v t.%V-2&#Jr-2£ n & . . (5') 
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Then equating coefficients in these two expressions (5') and (6) 

a s -|= -<&+|X+|(l -6pY) 
^_i = $ + p;+£(l-6pY) 



Whence 



also 



2 3_ a .2_ $_X — 1(1 — 6pY) 

= — #' 



y (7) 



The^e six equations (7) are the desired functions of x 9 y, z in terms of the auxiliary 
functions. 

I 2. The form of the spheroid as tidally distorted. 

The tide-generating potential has thirteen terms, each consisting of a solid harmonic 
of the second degree multiplied by a simple harmonic function of the time, viz. : three 
in <3>, three in <!>', three in ^, three in ^ and one in X. The fourteenth term of V 
can raise no proper tide, because it is independent of the time, but it produces a 
permanent increment to the ellipticity of the mean spheroid. 

Hence according to our hypothesis, explained in the introductory remarks, there 
will be thirteen distinct simple tides ; the three tides corresponding to <J>' may 
however be compounded with the three in <£, and similarly the W tides with the 
M? tides. Hence there are seven tides with speeds'* [2n—2S2, 2n, 2^ + 2/2], \ji-2S2, 
n, w+2/2], [2/2], and each of these will be retarded by its own special amount. 

The <& tides have periods of nearly a half-day, and will be called the slow, sidereal, 
and fast semi-diurnal tides, the ^ tides have periods of nearly a day, and will be called 
the slow, sidereal, and fast diurnal tides, and the X tide has a period of a fortnight, 
and is called the fortnightly tide. 

The retardation of phase of each tide will be called the "lag," and the height of 
each tide will be expressed as a fraction of the corresponding equilibrium tide of a 
perfectly fluid spheroid. Then the following schedule gives the symbols to be 
introduced to express lag and reduction of tide : — 

* The useful term " speed " is due, I believe, to Sir William Thomson, and is much, wanted to indicate the 
angular velocity of the radius of a circle, the inclination of which to a fixed radius gives the argument of 
a trigonometrical term. It will be used throughout this paper to indicate v, as it occurs in expressions of 
the type cos (vt+7]). 
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Tide . . . 


Semi-diurnal. 


Diurnal. 


Fortnightly. 
(2/3). 


/— -— '■■- ■■ 

Slow 

(2n-2f2). 


Sidereal 
(2»). 


Fast 
(271+2/2). 


Slow 

(71-2/2). 


Sidereal 


Fast 
(ti + 2/2). 


Height . 


■%i 


E 


R 2 




£' 


£', 


i?" 


JLJclff • . . 


2e x 


2e 


^2 


/ 

€ 




2e" 



The J^'s are proper fractions, and the e's are angles. 

Let r=a-\-cr be the equation to the surface of the spheroid as tidally distorted, a 
being the radius of the mean sphere,- — for we may put out of account the permanent 
equatorial protuberance due to rotation, and to the non-periodic term of V. 

It is a well known result that, if wr*& cos (vt-\-rj) be a tide-generating potential, 
estimated per unit volume of a homogeneous perfectly fluid spheroid of density w 9 

(8 being of the second order of surface harmonics), then the equilibrium tide due to this 

5 a 2 2q 

potential is given by o-=~-~S cos (vt+rj). If we write g=^" 3 this result may be 



S 



2g 



written -=- cos ( vt+rj). 

Now consider a typical term— say one part of the slow semi-diurnal term — of the 
tide-generating potential, as found in (3) : it was 

—wr^r^p* sin s cos 2<£ cos 2(w— /2). 

The equilibrium value of the corresponding tide is found by putting - equal to this 

•^ CO 

expression divided by wr^Q. 

Then if we suppose that there is a frictional resistance to the tidal motion, the tide 
will lag and be reduced in height, and according to the preceding definitions the 
corresponding tide of our spheroid is expressed by 

-= — E^jfi sin 3 9 cos 2</> cos [2(n— 12) — 2e Y ' 
co g 

All the other tides may be treated in the same way, by introducing the proper j^'s 

and e's. 

Thus if we write 



sir 

x e 



E x lp 4 eos (2n—2n-2e 1 )+EpY cos (2n—2e)+E z lq i cos (2n+2J2-2e 2 ) 
--E' l 2p s qcos(n — 2n — e' l )—E'2pq(p*—q*)cos(n—e')—E ii 2pq s cos (n+2J2-e a ) ^ (8) 

■-E"3pYcos(2n~2e") 



and if in the same symbols accented sines replace cosines, then, by comparison with (5), 
we see that 
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g<7 



r a 



(?-v^-2£n&-€&-v!y.+§(?+n i -H*)x< • • • ( 9 ) 



This is merely a symbolical way of writing down that every term in the tide- 
generating potential raises a lagging tide of its own type, but that tides of different 
speeds have different heights and lags. 

This same expression may also be written 



2** 



r a 



= -^ { <|,_IX 6 }_ 7? 2{-<E>_iX e }-C 2 |X-27 ? ^'-2^^-2^' e . . . (9') 



Then if we put- 



It is clear that 



Whence 



S a 



r a 



3/ d 



cl 



C— b: 

a-~c: 
b— a 

d= 



<£ e -X e 
23> 



"3^e 
& e 

.IaJ/ 



k. 



J 



■af a — hrf— c^+2d^^+-2e^+2f^7 . 



or 



~\ 



9iT \V ( j£-£x,h=-{( c - h )vl- d (v 2 -^)- e £v+H£} 



2T\ t df 



'drj/a 

cl \<7 






Of which expressions use will be made shortly. 



(10) 



• (11) 



(12) 



§ 3, TAe couples about the axes A, B, C caused by the moon's attraction. 

The earth is supposed to be a homogeneous spheroid of mean radius a, and mass iv 
per unit volume, so that its mass $f=f ttwcP. When undisturbed by tidal distortion 
it is a spheroid of revolution about the axis C, and its greatest and least principal 
moments of inertia are C, A. Upon this mean spheroid of revolution is superposed 
the tide-wave o\ 

The attraction of the moon on the mean spheroid produces the ordinary precessional 
couples 2r(G—A)yz 9 —2t(C—A)zx } about the axes A, B ? C respectively; besides 

MDCCCLXXIX, 3 N 
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these there are three couples, H, J$t, $L suppose, caused by the attraction on the 
wave surface cr. 

As it is only desired to determine the corrections to the ordinary theory of preces- 
sion,, the former may be omitted from consideration, and the attention confined to the 
determination of ft, JW, j£t 

The moon will be treated as an attractive particle of mass m. 

Now cr as defined by (9) is a surface harmonic of the second order; hence by the 
ordinary formula in the theory of the potential, the gravitation potential of the tide- 

wave at a point whose coordinates referred to A, B, C are rf, rrj, r£ is jttwciI - ) cr or 
¥~rcr. Hence the moments about the axes A, B, C of the forces which act on a 

ny> fyfff [ ft (T Ci (T\ 

particle of mass m, situated at that point, are i~ir(^7?~£~r)> & c -> & a Then if this 

particle has the mass of the moon; if r be put equal to c, the moon's distance; and if 
£ 7], I be replaced in cr by x, y, z (the moon's direction cosines) in the previous expres- 
sions, it is clear that —^Mariy^—z-—), &c, &c, are the couples on the earth caused 

by the moon's attraction. 

These reactive couples are the required %, JW, Ji. 

Hence referring back to (12) and remarking that %Ma^=C, the earth's moment of 
inertia, we see at once that 



— = — [(c— h)yz— d(y s — -z 2 ) — exy-\-£zx] 



S 

2 T 2 
= — [(a— o)zx— e(z*— x 1 ) — fyz+dxy] 

p =-^{(b— ajxy—f^—y 2 ) — dzx+eyz] 



j 



. . . . (13) 



Where the quantities on the right-hand side are defined by the thirteen equations 

(?) and (10). 

I shall confine my attention to determining the alteration in the uniform precession, 
the change in the obliquity of the ecliptic, and the tidal friction; because the nutations 
produced by the tidal motion will be so small as to possess no interest. 

In developing % and ffit I shall only take into consideration the terms with argu- 
ment n, and in J^ only constant terms ; for it will be seen, when we come to the 
equations of motion, that these are the only terms which can lead to the desired end. 

§ 4. Development of the couples % and Jj^t* 
Now substitute from (7) and (10) in the first of (13), and we have 

!~ 2 ^=~M<M-xj>m • ( 14 ) 
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A number of multiplications have now to be performed, and only those terms which 
contain the argument n to be retained. 

The particular argument n can only arise in six ways, viz. : from products of terms 
with arguments 2(n— 12), n — 2/2; 2n, n; 2(n+/2), fi+2/2; n— 2/2, 2/2; n-\-2fl, 2/2 and 
from terms of argument n multiplied by constant terms. 

If <£> and M? 1 , and <£>' and Y' be written underneath one another in the various com- 
binations in which they occur in the above expression, it will be obvious that the 
desired argument can only arise from terms which stand one vertically over the other; 
this renders the multiplication easier. The M?, X products are comparatively easy. 

Then we have 

(a) — J£®ip' = — £— E 1 p >1 q$in(n-2e l ) + 2EpY(p 2 — q l )sm(n — 2e) + E % pq 7 sm(n — 2e 2 )] 

08) +**> = +i[-E\p^qsm (n+c'j) + 2E'p\\p*-(f) sin (n+e*) + E' 2 pql&m (n+ c' a )] 

(y) —2*P<&' =same as (J3) 

(8) + 2^ ^ =same as (a) 

(e) -iX.^ = -i[E"6p\* sin (n-2c") -JE"6p*q* sin (n+2«")] 

(0 +^' e X = +i[JE\6p^ sin (n-O ~iWz 5 *** (n-c',)] 

fo) +i^.(l -6pY) = -W 2 P9:(P 2 -<f)( 1 -6pY) sin (n-e') 

Now put — =Fsinn+Grcosw. Then if the expressions (a), (/J) . . . (0 be added 

up when n=~, and the sum multiplied by — , we shall get F ; and if we perform the 

same addition and multiplication when n=0 9 we shall get G. 

In performing the first addition the terms (a) (8) do not combine with any other, 
but the terms (/3), (y), (£), (??) combine. 

Now 

p*q*(p*-q*) -ipq(p* - q%l - 6pY) = -$pq(p*-q*) (p 4 +2*~ 6p Y) 

— f p 5 2 3 +fp 3 2 5 = — %p 2f (f{p 2 —(f). 
Hence 

2t 3 

F : - = 



^E ± p 7 q cos 2e 1 —Ep 3 q s (p*—q 2 ) cos 2e—\E$>q l cos 2e 2 

"P> W- % 2 ) cob<! ^E'pqipP-q*) (p*+ q *-6pY) cose f ^E' 2 pq\3p^q^) cos^ 
— f E"p*q*(p*— q 2 ) cos 2c" (15) 



3 N 2 



458 MB. a. H. DARWIN ON THE PRECESSION OP A VISCOUS SPHEROID, 
Again for the second addition when n=0, we have 



p^^—f) + to(f> 3 — 2 s ) (i - 6p*f) =ipq(p~-q 2 ) 
TtPq 7 + ip V = ^pq & (3p 2 + <? 2 ) 



So that 



2t 2 



G-= = —\E-tfFq sin 2e 1 +l?p 3 g r3 (p 3 --g r2 ) sin 2e+-|-i£ 2 ^(/ 7 sin 2 



e 



2 



-i^LP 5 2(p 2 +3g 2 ) sin^+I^Mb 2 -^) 3 sin e'44^>2W+? 2 ) sin e\ 
+§ .£" jp 3 g 3 sin 2e" (16) 



And 



p = F sin w+Gcos n ......... (17) 



7T 



To find M it is only necessary to substitute n— - for n 9 and we have 

lift 

=— =— Fcosn+Gsinw ........ (18) 

Now there is a certain approximation which gives very nearly correct results and 
which simplifies these expressions very much. It has already been remarked that the 
three <&-tides have periods of nearly a half-day and the three ^-tides of nearly a day, 
and this will continue to be true so long as SI is small compared with n ; hence it may 
be assumed with but slight error that the semi-diurnal tides are all retarded by the 
same amount and that their heights are proportional to the corresponding terms in the 
tide-generating potential. That is, we may put e 1 =e 2 =c and E X =E 2 =E. The 
similar argument with respect to the diurnal tides permits us to put e\=e^=e and 

Then introducing the quantities P=p % -~(f= cos i, Q=2pq= &mi and observing 
that 

we have, 
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r 3 



F^-=l£P()(l-f <2 2 ) cos 2e-E'PQ(l-W) cos e'-%E"PQ s cos 2e" 






> ( 19 ) 



G-5— = -\EPQ(l-lQ?) sin 2e-f S'PQ 3 sin e'+f^'Q 8 sin 2e" 

§ 5. Development of the couple $L. 

In the couple Jfl about the axis of rotation of the earth we only wish to retain non- 
periodic terms, and these can only arise from the products of terms with the same 
argument. 

By substitution from (7) and (10) in the last of (13) 

^^ 2 ^=2^^-2^\^-^^+i%^ . (20) 

Then as far as we are now interested, 

2<S> e <E>'= — 2<E>' e *=jE' 1 |p 8 sin 2 ejl +Ep\f sin 2e-{-E i ^q s sin 2e a 
-K.^= i%V'=E\± P y sin e'^^W^-^ sin e'+i^Y sin e' a 



Hence 

C * 8 



ij^p 8 sin 2c 1 +jB4pV sin 2e+R 2 cj 8 sin 2€, 2 



+^2^8^+1^ . . . (21) 

If as in the last section we group the semi-diurnal and diurnal terms together and 
put E Y =Ec>=E y &c, and observe that 

2>H4pV+2 8 =(p*+2*) s +2 i >Y=(i-i<2 a ) 3 +i<2*=i >i! +l<2* 

2^Y+2i?Y(p s -2 3 ) 2 + 2 ^ 6 = 4 P 2 2 2 [p 4 '+2*-pY]=^ 3 ( 1 -|^ 2 )> 



then 



-f^=^(2»+tQ*) sin 2e+EQ*(l-lQ*) sin € ' .... (22) 



§ 6. 27^ equations of motion of the earth abouts its centre of inertia. 

In forming the equations of motion we are met by a difficulty, because the axes 
A, B ; C are neither principal axes, nor can they rigorously be said to be fixed in the 
earth. But M. Liouville has given the equations of motion of a body which is 
changing its shape, using any set of rectangular axes which move in any way with 
reference to the body, except that the origin always remains at the centre of inertia. 
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If A ; B ; 0, D ; E, F be the moments and products of inertia of the body about these 
axes of reference at any time ; H 1? H 2 , H 3 the moments of momentum of the motion 
of all the parts of the body relative to the axes ; o> l9 o> 2 , o) 3 the component angular 
velocities of the axes about their instantaneous positions, the equations may be written 

cl 
-- ( Ac^ ~ Fo) 2 — Ew 3 + Hj) + D (o) 2 3 — 0)\) + (C — B) 0) 2 0) 3 + F^gO^ — - Ew 2 o) 1 

+W2IT3— -to 3 H a =L (23) 

and two other equations found from this by cyclical changes of letters and suffixes. * 

Now in the case to be considered here the axes A, B, C always occupy the average 
position of the same line of particles, and they move with very nearly an ordinary 
uniform precessional motion. Also the moments and products of inertia may be 
written A+a', B+b' ? C+c', d', e', f, where a', h', c', d\ e', F are small periodic 
functions of the time and a'+b'+c'^O, and where A, B, C are the principal moments 
of inertia of the undisturbed earth, so that B is equal to A. 

Now the quantities a', b', &c, have in effect been already determined, as may be 
shown as follows : By the ordinary formulat the force function of the moon's action on 

the earth is — — +t( — I l where I is the moment of inertia of the earth about 

the line joining its centre to the moon, and is therefore 

=Aa?+ByHC^HaV+by+cV— 2dV«— 2e^— 2fa?y. 

But the first three terms of I only give rise to the ordinary precessional couples, and 
a comparison of the last six with (11) and (13) shows that 

a b c d e f g" " 

Also in the small terms we may ascribe to co l9 co 2 , co 3 their uniform precessional values, 
viz. : w x =: — II cos n 9 co 3 = — II sin n, co s =—n. 

When these values are substituted in (23), we get some small terms of the form 
all 2 sin n, and others of the form a'Tln sin n ; both these are very small compared to the 
terms in % and JM — the fractions which express their relative magnitude being 

112 Un 

— and — . 

T T 

There is also a term — nH 3 sin n, which I conceive may also be safely neglected, as 
also the similar terms in the second and third equations. 

It is easy, moreover, to show that according to the theories of the tidal motion 
of a homogeneous viscous spheroid given in the previous paper, and according to 

* Routh's ' Rigid Dynamics ' (first edition only), p. 150, or my paper in the Phil. Trans. 1877, Vol. 167, 
p. 272. The original is in Liouville's Journal, 2nd series, vol. iii., 1858, p. 1. 
f Route's ' Rigid Dynamics,' 1877, p. 495. 
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Sir William Thomson's theory of elastic tides, H l5 H 3 , H 3 are all zero, Those theories 
both neglect inertia but the actuality is not likely to differ materially therefrom. 

Thus every term where c^ and w 3 occur may be omitted and the equations reduced to 



A l -L./P 



dt 



— B)o). 2 o) 3 +n—+n 2 d / + -~~ 1 +^H 2 =: % 



dt 
cM 



dt 
dH» 



B^+(A— 0)0)30)!+^-— ?^ 2 e / +~7--7^H 1 = 



C^+(B 



dt 



-AX 



0) 



2 






» » 



. (24) 



As before with the couples, so here, we are only interested in terms with the argument 
n in the small terms on the left-hand side of the first two of equations (24), and in 
non-periodic terms in the last of them. 

Now for each term in the moon's potential, as developed in Section 1, there is (by 
hypothesis) a corresponding co -periodic flux and reflux throughout the earth's mass, 
and therefore the H x , H 3 , H 3 must each have periodic terms corresponding to each 
term in the moon s potential. Hence the only term in the moon's potential to be con- 
sidered is that with argument n, with respect to H x and H 2 in the first two equations ; 
and H 3 may be omitted from the third as being periodic. 

Suppose then that B^ was equal to h cos n-\-li sin n, then precisely as we found 



77- . f ^JJ 

from % by writing n—- for n we have K 2 =h sin n— hi cos n. Thus -j^+riH. 2 = 0, 



2 



dt 



dR 



2 



dt 



■nH 1 =0, and the H's disappear from the first two equations. 



Next retaining only terms in argument n in d' and e', we have from (10) 
e / =G^£J / pq(p 2 -—q^) cos (w — e'), & f =Q~-E / pq(]f—(f) sin (n—e) 

d& d& 

Therefore — +^d'~.0, — --— fte'=0, and these terms also disappear. 

(Jjv Cvu 

Lastly, put B=A, and our equations reduce simply to those of Euleu, viz. : 



A^+(C— A)co 2 co 3 =a 
A-~ ? — (C— A)g) 3 g) 1 =JW 



f^dco% 



-\ 



J 



* « * 



* * 1 4mJ*J f 



Now JB is small, and therefore <w 3 remains approximately constant and equal to — n 
for long periods, and as C— A is small compared to A, we may put a) n = —n in the first 
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two equations. But -when 0— A is neglected compared to 0, the integrals of these 
equations are the same as those of 

clco^ % dco 2 JJH da)^ 



dt~C di C 



^ ? 



— ~~ /-sj » » • • * * « • f^jD/ 



di c 



apart from the complementary function, which may obviously be omitted. The two 
former of (26) give the change in the precession and the obliquity of the ecliptic, and 
the last gives the tidal friction. 

c5 



§ 7. Precession and change of 'obliquity. 
Then by (17), (18), and (26) the equations of motion are 



dw^ 

dt ' 

d(0e> 

A , 

dt ' 



F sin ft+G cos n 



F cos ft+G sin ft 



* » > * 



* « » ft (jU § ) 



and by integration 



0) 



lr- 



n 



F cos ft+G sin ft], co 2 =-[— F sinft— Gcos ft] . . , , (28) 

1% 



But the geometrical equations (1) give 



di 

dt~ 

d'yjr 
lit 



os Y sin ft+a> 3 cos ft 



• 



sin ^=— cdi cos n— -a) Q sin n 

JL fV 



Therefore, as far as concerns non -periodic terms, 



di 
dt 



G dty . . F 
n dt n 



• • • » \ji(t7l 



If we wish to keep all the seven tides distinct (as will have to be done later), we 

fin (i'\Jf 

may write down the result for -rand—- from (15) and (16). 

But it is of more immediate interest to consider the case where the semi-diurnal 
tides are grouped together, as also the diurnal ones. In this case we have by (19) 



di r 3 



'^ P Q(l_ sq^e gin 2e+%PQ*E' sin e'-%QW sin 2e"} 



. . . (30) 



and since sin^=Q 



dty t 9 
dt ~~g?fc 



{±P(l-lQ*)E cos 2e-P(l-WW cos e'-%PQW cos 2e"} . (31) 
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In these equations P and Q stand for the cosine and sine of the obliquity of the 
ecliptic. 

Several conclusions may be drawn from this result. 

If e, e', e" are zero the obliquity remains constant. 

Now if the spheroid be perfectly elastic, the tides do not lag, and therefore the 
obliquity remains unchanged ; it would also be easy to find the correction to the 
precession to be applied in the case of elasticity. 

It is possible that the investigation is not, strictly speaking, applicable to the case 
of a perfect fluid ; I shall, however, show to what results it leads if we make the appli- 
cation to that case. Sir William Thomson has shown that the period of free vibration 
of a fluid sphere of the density of the earth, would be about 1 hour 34 minutes/" And 
as this free period is pretty small compared to the forced period of the tidal oscillation, 
it follows that E, E\ E" 9 will not differ much from unity. Then putting them equal 
to unity, and putting e, e ', e ; zero, since the tides do not lag, we find that the obliquity 
remains constant, and 

d i=-^P(l-IQ^)=-^co8i(l^8mH) ..... (32) 

clt %n z x 2 ^ } z %% \ 2 / \ ) 

This equation gives the correction to be applied to the precession as derived from the 
assumption that the rotating spheroid of fluid is rigid. This result is equally true if 
ail the seven tides are kept distinct. Now if the spheroid were rigid its precession 

would be — cos i, where e is the ellipticity of the spheroid. 

The ellipticity of a fluid spheroid rotating with an angular velocity n is f — or \ — ; 

but besides this, there is ellipticity due to the non-periodic part of the tide-generating 
potential. 

By (3) § 1 the non-periodic part of V is \wn s l^— cos 2 0)(l — 6p 2 cf) ; such a disturb- 
ing potential will clearly produce an ellipticity \ ~(l — fyp s y 2 ). 

If therefore we put ^ =l— , and remember that 6p 2 q 2 =% siirH', we have, 

e=e +ij(l-~f sin 2 *) 
Hence if the spheroid were rigid, and had its actual ellipticity, we should have 

-^-=— cos i+1 — ■ cos i(l — - I sin 2 i). . . . , , . (32') 

clt n 2 Wi v - ' v ; 

* Phil. Trans., 1863, p. 608. 
MDCCCLXXIX. 3 O 



464 MR. G. H. DARWIN ON THE PRECESSION OF A VISCOUS SPHEROID, 
Adding (32'j to (32), the whole precession is 

7 , — — COS v » » . « » » • • » • V^^/ 

at n 

We thus see that the effect of the non-periodic part of the tide-generating potential, 
which may be conveniently called a permanent tide, is just such as to neutralise the 
effects of the tidal action. The result (32") may be expressed as follows \- — 

The precession of a fluid spheroid is the same as that of a rigid one ivhich has an 
ellipticity equal to that due to the rotation of the spheroid. 

From this it follows that the precession of a fluid spheroid will differ by little from 
that of a rigid one of the same ellipticity, if the additional ellipticity due to the non- 
periodic part of the tide-generating influence is small compared with the whole 
ellipticity. 

Sir William Thomson has already expressed himself to somewhat the seme effect 

in an address to the British Association at Glasgow. * 



no & rjr 

Since e =-^~ , the criterion is the umallness of -^. 



It may be expressed in a different form; for ™ is small when —-7- ft is small compared 

7b Tb 

with 6, and — -f- n is the reciprocal of the precessional period expressed in days. Hence 

7\/ 

the criterion may be stated thus : The precession of a fluid spheroid differs by little 
from that of a rigid one of the same ellipticity, when the precessional period of the 
spheroid expressed in terms of its rotation is large compared with the reciprocal of 
its ellipticity. 

In his address, Sir William Thomson did not give a criterion for the case of a fluid 
spheroid without any confining shell, but for the case of a thin rigid spheroidal shell 
enclosing fluid he gave a statement which involves the above criterion, save that the 
ellipticity referred to is that of the shell itself; for he says, "The amount of this 
difference (in precession and nutation) bears the same proportion to the actual precession 
or nutation as the fraction measuring the periodic speed of the disturbance (in terms 
of the period of rotation as unity) bears to the fraction measuring the interior ellipticity 
of the shell. v 

This is, in fact, almost the same result as mine. 

This subject is again referred to in Part III. of the succeeding paper. 



*" See i Nature, 5 September 14, 1876, p. 429. The above statement of results, and the comparison with 
Sir William Thomson's criterion was added to the paper on September 17, 1879. 
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§ 8. The disturbing action of the sun, 

Now suppose that there is a second disturbing body, which may be conveniently 
called the sun.* 

* It is not at first sight obvious how it is physically possible that the sun should exercise an influence 
on the moon-tide, and the moon on the sun-tide, so as to produce a secular change in the obliquity of the 
ecliptic and to cause tidal friction, for the periods of the sun and moon about the earth are different. It 
seems, therefore, interesting to give a physical meaning to the expansion of the tide-generating potential ; 
it will then be seen that the interaction with which, we are here dealing must occur. 

The expansion of the potential given in Section 1 is equivalent to the following statement : — 
The tide-generating potential of a moon of mass m, moving in a circular orbit of obliquity i at a 
distance c, is equal to the tide-generating potential of ten satellites at the same distance, whose orbits, 
masses, and angular velocities are as follows : — 

* 

1. A satellite of mass m cos 4 -, moving in the equator in the same direction and with the same angular 

velocity as the moon, and coincident with it at the nodes. This gives the slow semi-diurnal tide of 
speed 2 (n — 0). 

2. A satellite of mass m sin 4 -, moving in the equator in the opposite direction from that of the moon, 

but with the same angular velocity, and coincident with it at the nodes. This gives the fast semi-diurnal 
tide of speed 2(^ + 0). 

• • 

3. A satellite of mass m2sin 3 -cos 2 -, fixed at the moon's node. This gives the sidereal semi-diurnal 

2 2 

tide of speed 2n. 

♦ • 



4. A repulsive satellite of mass — m.2 sin - cos 3 -, moving in N. declination 45° with twice the moon's 

angular velocity, in the same direction as the moon, and on the colure 90° in advance of the moon, when 
she is in her node. 

c 

5. A satellite of mass m sin i cos 3 -, moving in the equator with twice the moon's angular velocity, and 

in the same direction, and always on the same meridian as the fourth satellite. (4) and (5) give the slow 
diurnal tide of speed n— 20. 

6. A satellite of mass m sin 3 - cos—, moving in N. decimation 45° with twice the moon's angular velocity, 

but in the opposite direction, and on the colure 90° in advance of the moon when she is in her node. 

i ' * 

7. A repulsive satellite of mass — m. ^sin 3 ^-cos-, moving in the equator with twice the moon's angular 

u u u 

velocity, but in the opposite direction, and always on the same meridian as the sixth satellite. (6) and 
(7) give the fast semi-diurnal tide of n + 20. 

8. A satellite of mass m sin i cos i fixed in N. declination 45° on the colure. 

9. A repulsive satellite of mass — m. ~ sin i cos i 9 fixed in the equator on the same meridian as the eighth 

satellite. (8) and (9) give the sidereal diurnal tide of speed n. 

10. A ring of matter of mass m, always passing through the moon and always parallel to the equator. 
This ring, of course, executes a simple harmonic motion in declination, and its mean position is the 
equator. This gives the fortnightly tide of speed 20. 

Now if we form the potentials of each of these satellites, and omit those parts which, being indepen- 
dent of the time, are incapable of raising tides, and add them altogether, we shall obtain the expansion 
for the moon's tide- generating potential used above ; hence this system of satellites is mechanically 

3 O 2 
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II cosec i must henceforth be taken as the full precession of the earth, and the time 
may be conveniently measured from an eclipse of the sun or moon. Let m n c / be the 
sun's mass and distance ; il t the earth's angular velocity in a circular orbit ; and let 



3 m, 



T/ 2 cf 

It would be rigorously necessary to introduce a new set of quantities to give the 
heights and lagging of the seven solar tides : but of the three solar semi-diurnal tides, 
one has rigorously the same period as one of the three lunar semi-diurnal tides 
(viz. : the sidereal semi-diurnal with a speed 27i), and the others have nearly the same 
period ; a similar remark applies to the solar diurnal tides. Hence we may, without 
much error, treat E, e, E\ e as the same both for lunar and solar tides ; but E'" 9 e'" 
must replace E f \ e", because the semi-annual replaces the fortnightly tide. 

Then if new auxiliary functions <&„ M* /? X, be introduced, the whole tide-generating 
potential V per unit volume of the earth at the point rg } rrj y r£ is given by 

If then, as in (10), we put 

c— b=<l> e +X e , &c, c / — b^^+X^, &c, 

the equation to the tidally-distorted earth is r=a+o , +0" / > where 

equivalent to the action of the moon alone. The satellites 1, 2, 3, in fact, give the semi-diurnal or 
<I> terms ; satellites 4, 5, 6, 7, 8, 9 give the diurnal or ^r terms ; and satellite 10 gives the fortnightly or 
X term. 

This is analogous to " GrAUSs's way of stating the circumstances on which * secular' variations in the 
elements of the solar system depend ;" and the analysis was suggested to me by a passage in Thomson and 
Tait's 'Nat. Phil.,' § 809, who there refer to the annular satellite 10. 

It will appear in Section 22 that the 3rd, 8th, and 9th satellites, which are fixed in the heavens and 
which give the sidereal tides, are equivalent to a distribution of the moon's mass in the form of a uniform 
circular ring coincident with her orbit. And perhaps some other simpler plan might be given which 
would replace the other repulsive satellites. 

These tides, here called "sidereal," are known, in the reports of the British Association on tides for 1872 
and 1876, as the K tides. 

In a precisely similar way, it is clear that the sun's influence may be analysed into the influence of 
nine other satellites and one ring, or else to seven satellites and two rings. Then, with regard to the 
interaction of sun and moon, it is clear that those satellites of each system which are fixed in each system 
(viz. : 3, 8, and 9), or their equivalent rings, will not only exercise an influence on the tides raised by 
themselves, but each will necessarily exercise an influence on the tides raised by the other, so as to produce 
tidal friction. All the other satellites will, of course, attract or repel the tides of all the other satellites 
of the other systems ; but this interaction will necessarily be periodic, and will not cause any interaction in 
the way of tidal friction or change of obliquity, and as such periodic interaction is of no interest in the 
present investigation it may be omitted from consideration. In the analysis of the present section, this 
omission of all but the fixed satellites appears in the form of the omission of all terms involving the moon's 
or sun's angular velocity round the earth. 
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8_?__ a fa Sec, §■ ^=-*&, &c. 

T CI T k (I 

Also if x } y, z and a?,, y n % t be the moon's and sun's direction cosines, we have as 
in (7), 

^-f=$+X+l(l-6|)Y), &c., y^—zf—% +X / +l(l-6pY), &c. 

Then using the same arguments as in Section 3, the couples about the three axes in 
the earth may be found, and we have 



o : 



{In, ff \ I (T (T \ In u \ I CT (T \ l 



where in the first term x, y, % are written for £ 77, £ in cr+cr,, and in the second term 
a? y , 2//5 % , are similarly written for £ 77, £. 

Now let % m *, % m ^ %mm y indicate the parts of the couple % which depend on the 
moon's action on the lunar tides, the sun's action on the solar tides, and the moon's 
and sun's action, on the solar and lunar tides respectively, then 



% %m? , %m * f / d d\a t , / d d\ 



{I d d W, , / d 



<T 



C — C + "C 1 T l% 2 %jfl ^ T \ V/ dz Z 'dyja 



Then obviously 



36/mm/ . ^ 



^ — / ==(c—b)y / « / +(c / --b / )y2f+ &c. 



C S 



As before, we only want terms with argument n in % mM/i JJft^, and non-periodic 
terms m "j^% mmj - 

The quantities a, b, &c., as, 3/, 2 with suffixes differ from those without in having 
fL / in place of fl } and it is clear that no combination of terms which involve fl y and 
12 can give the desired terms in the couples. Hence, as far as % mm/9 ffitmnm i$mm, are 
concerned^ the auxiliary functions may be abridged by the omission of all terms 
involving fl or fl /t 

Therefore ; from (4), we now simply have 

<&=® / =p2qZ cos 2n 9 M r =*I> , / = — 2p$(p 2 — q 2 ) cos n, X=X y =0. 

But c— b only differs from c / — -b / in that the latter involves il f instead of S2 P and the 
same applies to yz and yp t . 

Hence, as far as we are now concerned, 

(z-h)y ( z={c-h)yz 

and similarly each pair of terms in % m?// are equal inter se. 
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Ihus 



%mm/ ' 4 ^==(c-b)y»--d(^-^)-e^+fe». 



Comparing with (14), when X is put equal to zero, we have 

al UL.I-TT 



4mm/ . 



o ' s 



4*.^+K. {*+i(l -6PV)1 -£*,*'+£*>• 



This quantity may be evaluated at once by reference to (15), (16), and (17), for it is 
clear that % mm/ is what % m * becomes when E 1 =^E 9/ =0 } E'^E^O, and when 2^ 
replaces r 2 . 

If, therefore, we put-^=F WOT/ .sin n+G^, cos n, and remark that 



C 



^<f{f-f) = l>PQ\ 2pg(p 3 -^)^+2 4 ^6^Y)=:Pg(l-2^), 2pq(p*--qZy=P 2 Q, 
we have by selecting the terms in jE 7 , A 7/ out of (15) and (16), 



TT A 



F m/ ~ ^= -\EPQ* cos 2e-E , PQ(l~-2Q*) cos e' 
8 



V . . . « [dd) 



T7\ 



G mm/ +-^= \EPQ* sin 2e+ET z Q sin e 



It may be shown in a precisely similar way by selecting terms out of (21) that 

^^ T ^^EQHm2e+E / P^Bme / (34) 

It is worthy of notice that (33) and (34) would be exactly the same, even if we did 
not put E^—E^E; E\=E\=E' ; e^e^e; e\=e 2 =6, because these new terms 
depend entirely on the sidereal semi-diurnal and diurnal tides. The new expressions 
which ought rigorously to give the heights and lagging of the solar semi-diurnal and 
diurnal tides would only occur in % m ». 

In the two following sections the results are collected with respect to the rate of 
change of obliquity and with respect to the tidal friction. 



§ 9. The rate of change of obliquity due to both sun and moon. 

An 

The suffixes m 2 , mf, n%m i to — will indicate the rate of change of obliquity due to 

the moon alone, to the sun alone, and to the sun and moon jointly. 

Then writing for P and Q their values, cos i and sin i, we have by (19) and (29), or 
by (30), 
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-f -—=1 s in i cos i(l — § sin 2 i)E sin 2e+f sin 3 i cos ii?' sin e— f sin 3 ijE7" sin 2 



• 3 ^ 



e 



-| ~^~= | sin i cos i(l — f sin 3 i)E sin 2e+f sin 3 i cos iE' sin c'— § sin 3 £/£'" sin 2 



W35) 



T/ 3 * 



e 



and by (33) and analogy with (19) and (29) 



<yyft ff'7 

_a -~pz= — 1 sin 3 i cos -ii£ sin 2e— sin i cos 3 iE' sin e' . . . . (36) 



TT, dt 



The sum of these three values of -- gives the total rate of change of obliquity due 

both to sun and moon, on the assumption that the three semi-diurnal terms may be 
grouped together, as also the three diurnal ones. 

It will be observed that the joint effect tends to counteract the separate effects ; this 
arises from the fact that, as far as regards the joint effect, the two disturbing bodies 
may be replaced by rings of matter concentric with the earth but oblique to the 
equator, and such a ring of matter would cause the obliquity to diminish, as was shown 
in the abstract of this paper (Proe. Roy. Soc, No. 191, 1878), by general considera- 
tions, must be the case. 



§ 10. The rate of tidal friction due to both sun and moon. 

The equation which gives the rate of retardation of the earth's rotation is by (26) 

— -=— ; it will however be more convenient henceforward to replace o) s by — n and 

to regard n as a variable, and to indicate by n the value of n at the epoch from 
which the time is measured. 

Generally the suffix to any symbol will indicate its value at the epoch. 

Then the equation of tidal friction may be written 

*\V — ciwo Gn o G X ' * ^ ' 

Then by (22) and (34), in which the semi-diurnal and diurnal terms are grouped 
together, we have 



'2j)fe= (cos 3 i+l sin 4 i)E sin 2e+ sin 3 i(l - f sin 3 i)E sin /= fe)^f ! 
— ) — — =-|l sin 4 iE sin 2e+ sin 3 i cos 3 iE' sin e 



■ (38) 
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§ 11. The rate of change of obliquity ivhen the earth is viscous, 

In order to understand the physical meaning of the equations giving the rate of 
change of obliquity (viz.: (35) and (36) if there be two disturbing bodies, or (29) if 
there be only one) it is necessary to use numbers. The subject will be illustrated in 
two cases: first, for the sun, moon, and earth with their present configurations; 
and secondly, for the case of a planet perturbed by a single satellite. For the first 
illustration I accordingly take the following data: ^=32*19 (feet, seconds), the earth's 
mean radius a=20*9xl0 6 feet, the sidereal day '9973 m. s. days, the sidereal year 
= 365*256 m. s. days, the moon's sidereal period 27*3217 m. s. days, the ratio of the 
earth's mass to that of the moon p=-82 } and the unit of time the tropical year 365*242 
m. s. days. 



Then we have 



f£ =:27r-f- , 9978 in radians per m. s. day 

^ 5a 

r=f X sV of 47T 3 -" (month) 2 
T/ =f of 47r 3 -r- (sidereal year) 3 . 



Then it will be found that 



— =-6598 degrees per million tropical years 

8^7 C? X At/ 

fmmfmm 

— ^z=-3064 

rtVM >* J5 5? 



• » . . » \0«7f 



These three quantities will henceforth be written u\ u~ 9 uu r 

For the purpose of analysing the physical meaning of the differential equations for 



O 9 



fin // I i)%\ 

— and -r ( — , no distinction will be made between — and — -, &c,, for it is here only 

dt dt\nj %n pi Q 

sought to discover the rates of changes. But when we come to integrate and find the 
total changes in a given time, regard will have to be paid to the fact that both r and 
n are variables. 

For the immediate purpose of this section the numerical values of u% u*, / uu / given 
in (39), will be used. 

I will now apply the foregoing results to the particular case where the earth is a 
viscous spheroid. 

Let p=-T7~~, where v is the coefficient of viscosity. 
r 19v J 
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Then by the theory of bodily tides as developed in my last paper 



E— cos 2e, E= cos e', W— cos 2e", J^ cos 2 c'" 



jW / & • / & / / jW U« y / / ^ m Jul t 

tan 2e= — , tan e = -, tan 2e = — , tan 2e = — 
P P P P 



> . . . . (40) 



Rigorously, we should add to these 

Ej= cos 2e v E%= cos 2e 2 , JS?'^ cos e\, Ef%— cos c'g 

x rt 2(^-/2) , rt 2(^4-/2) , , w-2/2 , , ^ + 2/2 
tan 26!=— -, tan 2c 3 = , tan € x = , tan e 2 = 



> . (40') 



But for the present we classify the three semi-diurnal tides together, as also the 
three diurnal ones. 
Then we have 



di 

di 



[_\ sin i cos i(l — f sin s i) sin 4e+f sin 3 i cos i sin 2e](u*+u*) — t<s sin 3 i ^ n ^ f u % 
— -^ sin 8 i sin 4e"%/— (£ sin 3 i cos i sin 4e+| sin i cos 8 i sin 2e)uu / . 



Now 



| sin { cos i(l — I sin 2 i)=w$ sin 2i(5 + 3 cos 2i)=-^(5 sin 2i+f sin 4i) 
| sin 3 1 cos i= -^ sin 2i(l — cos 2i)=~i € (2 sin 2i— sin 4i) 
-^ sin 3 i=-^f (3 sin i— sin 3i), ^ sin 8 1 cos i=-^(2 sin 2i— sin &i) 
| sin i cos 3 1=| sin 2i(l + cos 2i)=-g^(2 sin 2i+ sin 4i). 



di 



If these transformations be introduced, the equation for — may be written 



clt 



"7 * 

64y = 

dt 



— 9(w* sin 4c // +^ / 2 si» 4 c'") sin i+3(i* a sin 4c"+w / a sin 4c'") sin 3i 
+ [(5 .sin 4e+6 sin 2c / )(w a +w / 2 )'— (4 sin 4e+8 sin 2e')m* / ] sin 2i 
+[(f sin 4c— 3 sin 2e)(v?-\-u?) + (2 sin 4e— 4 sin 2€')tm y ] sin ii 



. . (41) 



Then substituting for u and u t their numerical values (39), and omitting the term 
depending on the semi-annual tide as unimportant, I find 



64 



di 
dt 



MDCCCLXXIX, 



— 5*9378 sin 4e" sin {+1*9793 sin 4 € " sin 3i 
+ {2*7846 sin 4e+2*3611 sin 2e} sin 2% 
+ {1*8159 sin 4e— 3*6317 sin 2e] sin U 

O -p 



V 



(42) 
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The numbers" are such that j- is expressed in degrees per million years. 

The various values which — is capable of assuming as the viscosity and obliquity 

vary is best shown graphically. In Plate 36, figs. 2 and 3, each curve corresponds 
to a given degree of viscosity, that is to say to a given value of c, and the ordinates 

give the values of y as the obliquity increases from 0° to 90°. The scale at the side 

of each figure is a scale of degrees per hundred million years — e.g., if we had e=30° 
and i about 57°, the obliquity would be increasing at the rate of about 3° 45' per 
hundred million years. 

The behaviour of this family of curves is so very peculiar for high degrees of 
viscosity, that I have given a special figure (viz. : Plate 36, fig. 3) for the viscosities 
for which 6=40°, 41°, 42°, 43°, 44°. 

The peculiarly rapid variation of the forms of the curves for these values of e is due 
to the rising of the fortnightly tide into prominence for high degrees of viscosity." 
The matter of the spheroid is in fact so stiff that there is not time in 12 hours or a 
day to raise more than a very small tide, whilst in a fortnight a considerable lagging 
tide is raised. 

For e=44° the fortnightly tide has risen to give its maximum effect (i.e., sin 4e"=l), 
whilst the effects of the other tides only remain evident in the hump in the middle of 
the curve. Between e=44° and 45° the ordinates of the curve diminish rapidly and 
the hump is smoothed down, so that when e=45° the curve is reduced to the 
horizontal axis. 

By the theory of the preceding paper, * the values of e when divided by 15 give 
the corresponding retardation of the bodily semi-diurnal tide — e.g., when e=30° the 
tide is two hours late. Also the height of the tide is cos 2e of the height of the 
equilibrium tide of a perfectly fluid spheroid — e.g., when e=30° the height of tide is 
reduced by one-half. In the tables given in Part I., Section 7, of the preceding paper, 
will be found approximate values of the viscosity corresponding to each value of e. 

The numerical work necessary to draw these figures was done by means of Orelle's 
multiplication table, and as to fig. 2 in duplicate mechanically with a sector ; the ordi- 
nates were thus only determined with sufficient accuracy to draw a fairly good figure. 

For the two figures I found 108 values of each of the seven terms of— (nine values 

do 

of i and twelve of e), and from the seven tables thus formed, the values corresponding 
to each ordinate of each member of the family were selected and added together. 

From this figure several remarkable propositions may be deduced. When the 
ordinates are positive, it shows that the obliquity tends to increase, and when 
negative to diminish. Whenever, then, any curve cuts the horizontal axis there is a 
position of dynamical equilibrium ; but -when the curve passes from above to below, it 

* " On the Bodily Tides of Viscous and Semi-elastic Spheroids," &c, Phil, Trans., 1879, Part I. 
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is one of stability, and when from below to above, of instability. It follows from tliis 
that the positions of stability and instability must occur alternately. When e=0 or 
45° (fluidity or rigidity) the curve reduces to the horizontal axis, and every position of 
the earth's axis is one of neutral equilibrium. 

But in every other case the position of 90° of obliquity is not a position of equi- 
librium, but the obliquity tends to diminish. On the other hand, from e~0° to about 
30° (infinitely small viscosity to tide retardation of two hours), the position of zero 
obliquity is one of dynamical instability, whilst from then onwards to rigidity it 
becomes a position of stability. 

For viscosities ranging from e=0° to about 42|° there is a position of stability which 
lies between about 50° to 87° of obliquity ; and the obliquity of dynamical stability 
diminishes as the viscosity increases. 

For viscosities ranging from €=30° nearly to about 42|°, there is a second position 
of dynamical equilibrium, at an obliquity which increases from 0° to about 50°, as the 
viscosity increases from its lower to its higher value. But this position is one of 
instability. 

From e= about 42^° there is only one position of equilibrium, and that stable, viz. : 
when the obliquity is zero. 

If the obliquity be supposed to increase past 90°, it is equivalent to supposing the 
earth's diurnal rotation reversed, whilst the orbital motion of the earth and moon 
remains the same as before ; but it did not seem worth while to prolong the figure, as 
it would have no applicability to the planets of the solar system. And, indeed, the 
figure for all the larger obliquities would hardly be applicable, because any planet 
whose obliquity increased very much, must gradually make the plane of the orbit of 
its satellite become inclined to that of its own orbit, and thus the hypothesis that the 
satellite's orbit remains coincident with the ecliptic would be very inexact. 

It follows from an inspection of the figure that for all obliquities there are two 
degrees of viscosity, one of which will make the rate of change of obliquity a maximum 
and, the other minimum. A graphical construction showed that for obliquities of about 
5° to 20°, the degree of viscosity for a maximum corresponds to about €=17-^°^ 
whilst that for a minimum to about €=40°. In order, however, to check this con- 
clusion, I determined the values of e analytically when {=15°, and when the 
fortnightly tide (which has very little effect for small obliquities) is neglected. I 
find that the values are given by the roots of the equation 

x*+l0x*+i3'660x — 20-412 = 0, where x=3 cos 4e 
This equation has three real roots, of which one gives a hyperbolic cosine, and the 

* I may here mention that I found when e=17J°, that it would take about a thousand million years for 
the obliquity to increase from 5° to 23^°, if regard was only paid to this equation of change of obliquity. 
The equations of tidal friction and tidal reaction will, however, entirely modify the aspects of the case, 

3 P 2 
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other two give e= 1 8° 1 5' and €=41° 37'. This result therefore confirms the geometrical 
construction fairly well. 

It is proper to mention that the expressions of dynamical stability and instability 
are onlv used in a modified sense, for it will be seen when the effects of tidal friction 
come to be included, that these positions are continually shifting, so that they may be 
rather described as positions of instantaneous stability and instability. 

* I will now illustrate the case where there is only one satellite to the planet, and 
in order to change the point of view, I will suppose that the periodic time of the satellite 
is so short that we cannot classify the semi-diurnal and diurnal terms together, but 
must keep them all separate. 

Suppose that n=5I2; then the speeds of the seven tides are proportional to the 
following numbers, 8, 10, 12 (semi-diurnal); 3, 5, 7 (diurnal) ; 2 (fortnightly). 

These are all the data which are necessary to draw a family of curves similar to 
those in Plate 36, figs. 2 and 3, because the scale, to which the figure is drawn, is 
determined by the mass of the satellite, the mass and density of the planet, and the 
actual velocity of rotation of the planet. 

Then by (16) and (29) we have 



— zzz—^ptq sin Ae l —p s q s (p^—tf) sin 4e— ^pq 1 sin 4e 2 — |j9 3 g 3 sin 4c" 

+^p 5 g'(_p 2 +3g 2 ) sin %t\-- : kp(l{p % -- ( ff sin 2c' — \p<f {%P* + (f) sin 2c' a ] 

where p= cos - and q= sin - 

This equation may be easily reduced to the form 

- = -4 T sin i\ [10 sin 4e 1 — 10 sin 4e 2 +16 sin 2e' 1 — 16 sin 2€ / a — 12 sin 4c"] 

dt %n l 

+ cos i[l 5 sin 4e x — 4 sin 4e+ 1 5 sin 4e 2 + 1 8 sin 2e\ — 24 sin 2e'+ 1 8 sin 2e 2 

+ cos 2i[6 sin Ae l — 6 sin 4^+12 sin 4e"] 

+ cos 3i[sm 4e x +4 sin 4e+ sin 4e 2 — -2 sin 2e / 1 — 8 sin 2e'— 2 sin 2c / 2 ] > 

m 

which is convenient for the computation of the ordinates of the family of curves which 

d% 
illustrate the various values of ~ for various obliquities and viscosities. 

(jjjj 

In Plate 36, fig. 4, the lag (e) of the sidereal semi-diurnal tide is taken as the 
standard of viscosity. The abscissae represent the various obliquities of the planet's 

equator to the plane of the satellite's orbit ; the ordinates represent the values of — 

the actual scale depending on the value of — ) ; and each curve represents one degree 
of viscosity, via. : when €=10° 9 20°, 30, 40° and 44°. 

* From here to the end of the section was added July 8, 1879. 
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The computation of the ordinates was done by Crelle's three-figure multiplication 
table, and thus the figure does not profess to be very rigorously exact. 

This family of curves differs much from the preceding one. For moderate 
obliquities there is no degree of viscosity which tends to make the obliquity diminish, 
and thus there is no position of dynamically unstable equilibrium of the system 
except that of zero obliquity. Thus we see that the decrease of obliquity for small 
obliquities and large viscosities in the previous case was due to the attraction of the 
sun on the lunar tides and the moon on the solar tides. 

In the present case the position of zero obliquity is never stable, as it was before. 
The dynamically stable position at a large obliquity still remains as before, but in 
consequence of the largeness of the ratio Sl-r-n (^th instead of -^-th), this obliquity of 
dynamical stability is not nearly so great as in the previous case. As the ratio Sl-r-n 
increases, the position of dynamical stability is one of smaller and smaller obliquity, 
until when Sl-r-n is equal to a half, zero obliquity becomes stable, — as we shall see 
later on. 

§ 12. Rate of tidal friction when the earth is viscous. 

If in the same way the equations (37) and (38) be applied to the case where the 
earth is purely viscous, when the semi-diurnal and diurnal tides are grouped together, 
we have 

~dt\^J = ( M *+ u ?)lU GO ** *+f sin4 l ) sin 4e +i sinS i( l "I sinH) sin 2e'] 1 
"\- uu [a s ^ i s i n ^ € +i sin 3 i cos 2 i sin 2c 7 ] J 

d f i% \ 

Plate 36, fig, 5, exhibits the various values of — — for the various obliquities and 

degrees of viscosity, just as the previous figures exhibited ~. The calculations were 

done in the same way as before, after the various functions of the obliquity were 
expressed in terms of cos 2i and cos Ai. 

The only remarkable point in these curves is that, for the higher degrees of 
viscosity, the tidal friction rises to a maximum for about 45° of obliquity. The tidal 
friction rises to its greatest value when e=22^° nearly; this is explained by the fact 
that by far the largest part of the friction arises from the semi-diurnal tide, which has 
its greatest effect when sin 4e is unity. 

§ 13. Tidal friction and apparent secular acceleration of the moon. 

I now set aside again the hypothesis that the earth is purely viscous, and return to 
that of there being any kind of lagging tides. 

I shall first find at what rate the earth is being retarded when it is moving with its 
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present diurnal rotation, and when the moon is moving in her present orbit, and no 
distinction will be made between n and n Q ; all the secular changes will be considered 
later. 

The numerical data of Section 11 are here used, and the obliquity of the ecliptic 
i=2#° 28'; then u and u f being expressed in radians per tropical year, I find 



di /DOO J-.1 * r\ I 0±~tO y-j/ * f 

— = --—— jfa 8U1 2e+ — ~7-Jb sm e 
C 10° 10° 

M 1-1978 „ . . . -2669^, . , 

^i = ~io^ Esm2e+ ~i¥' jE smc 



. . (44) 



Then integrating the equation (37) and putting n=n , when t=-0 

Integrating a second time, we find that a fixed meridian in the earth has fallen 
behind the place it would have had, if the rotation had not been retarded, by 

1^ seconds of arc. And at the end of a century it is behind time 

C 7T 

1900*27.£ T sin 2e+423'49i£ / sin e m. s. seconds of time. 

If the earth were purely viscous, and when e=l7° 30'* (which by Section 11 
causes the rate of change of obliquity to be a maximum), I find that at the end of a 
century the earth is behind time in its rotation by 17 minutes 5 seconds. 

By substitution from the second of (44), equation (45) may be written in the form 



1-1978 . -2669 

To 8 "* sm2c io* 



n=n Q ( 1— -~^~tE sin. 2e— ~^-tE sine') . . . . . . (46) 



which in the supposed case of pure viscosity when e=17° 30' becomes 

n=n Q [l yp-tj (47) 

All these results would, however, cease to be even approximately true after a few 

millions of years. 

The effect of the failure of the earth to keep true time is to cause an apparent 
acceleration of the moon's motion ; and if the moon's motion were really unaffected by 

* This calculation was done before I perceived that I had not chosen that degree of viscositj which 
makes the tidal friction a maximum, but as all the other numerical calculations have been worked out for 
this degree of viscosity I adhere to it here also. 
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the tides in the earth, there would be an apparent acceleration of the moon in a 
century of 

1O43"-28J0 sin 2e+232"-50-#' sine 7 (48) 

for the moon moves over // *5490 of her orbit in one second of time. 

This apparent acceleration would however be considerably diminished by the effects 
of tidal reaction on the moon, which will now be considered. 



14. Tidal reaction on the moon. 



The action of the tides on the moon gives rise to a small force tangential to the 
orbit accelerating her linear motion. The spiral described by the moon about the 
earth will differ insensibly from a circle, and therefore we may assume throughout 
that the centrifugal force of the earth's and moon's orbital motion round their common 
centre of inertia is equal and opposite to the attraction between them. 

We shall now find the tangential force on the moon in terms of the couples which 
we have already found acting on the earth. Those couples consist of the sum of three 
parts, viz. : that due (i) to the moon alone, (ii) to the sun alone, and (hi) to the action 
of the sun on the lunar tides and of the moon on the solar tides, the latter two being 
equal inter se. 

Now since action and reaction are equal and opposite, therefore the only parts of 
these couples which correspond with the tangential force on the moon are those which 
arise from (i), and one-half those which arise from (iii). 

We may thus leave the sun out of account if we suppose the earth only to be acted 
on by the couples % m *+^%nm, 9 JW^+^JWL,,, i£L*+ii$L^; these couples will be 
called %', ffil', ffi 9 and the part of the change of obliquity which is due to %', ffit' 

di' 

will be called -r. 

at 

Let r and —ft be the moon's distance, and angular velocity at any time, and v the 
ratio of the earth's mass to the moon's. 

Let T be the force which acts on the moon perpendicular to her radius vector, in 
the direction of her motion. 

From the equality of action and reaction, it follows that Tr must be equal to the 
couple which is produced by the moon's action on the tides in the earth, acting in the 
direction tending to retard the earth's diurnal rotation about the normal to the ecliptic. 
Referring to Plate 36, fig. 1, we see that the direction cosines of this normal are 
— sin i cos n, — sin i sin n, cos i; hence 

Tr= — sin i(%' cos n+ffil' sin n)-\-ffi cos i. 

* This section has been partly rewritten and rearranged since the paper was presented. (Dec. 19, 1878.) 
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But by (17) and (18) 



^ — \E m* m T"%£ mm,) ^^^i\\^m^i~2^mm / ) ^^S 71 



c 

G 



= ~ (F w9 + PU V ) cos w+(G w9 +JG» w »,) sin w. 



Hence 



Thus 



™ " OOk5 /t'"T" ™ Dxll /i — \Jjm! "t" q VJm — Iv 7 • 

C At ({J 



Tr=C^J — cos i+n sin ^— ^ ....... . (49) 



In order to apply the ordinary formula for the motion of the moon, the earth must 
be reduced to rest, and therefore T must be augmented by the factor (M+iyij-t-M. 
Then if 3- be the moon's longitude, the equation of motion of the moon is 

m l( T 2 d &)-^J n T (50) 

m dty dt)~ M X1 K V) 

But since the orbit is approximately circular — =,0. 



Also rn=C-T-iva 2 , and 



dt 
M + m l-\-v 
v 



Therefore by (49) and (50) 

- 7/ -=-fycr— — i —pr cos ^+7^ sm ?,— \ 
dt * ^ [ G tftf J 

Now let £=(-77) , whence n?=n 2 +£; 6 . 

The suffix to fl indicates the value of fl when the time is zero, and no confusion 
will arise by this second use of the symbol £ 

But since the centrifugal force is equal to the attraction between the two bodies, 
and the orbit is circular, therefore flh°=M~\-m. 

So that fl h*=(M+m)g Q . 

Tb erefore 

r s =(M"+m) f £ 4 /2 -*, and th*=(M+m)*n Q - i (; 
and hence 

But M+m=ga 2 , because M and m are here measured in astronomical units of 

v 

mass* 
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Therefore our equation may be written 



ga 



1 +V 



V 



^^j t =^K l + v ) 



~ eosi+nsinf^- 
(J at 



JN ow Jet 



r> I J* 



J 1^(1+ v) 



1 % 

and let sn^o*^-, and let N= — . 



/* 



n, 



o 



And we have 






<i£ Gn 



o 



cfa 



• • » 



(51) 



(52) 



It is not hard to show that the moment of momentum of the orbital motion of the 
two bodies is C-f-6\f2% and that of the earth's rotation is obviously On. Hence 
mil* is the ratio of the two momenta, and /x, is the ratio of the two momenta at 
the fixed moment of time, which is the epoch. 

In the similar equation expressive of the rate of change in the earth's orbital motion 
round the sun, it is obvious that the orbital moment of momentum is so very large 
compared with the earth's moment of momentum of rotation, that /z, is very large and 
the earth's mean distance from the sun remains sensibly constant (see Section 19). 

Then by (16) and (29), remembering that 



% 



• V lb Off, '2 

p= cos -, q= BID. -, — 



®'"- nnrl AT- % 



n 



n, 



o 



we have 



2 



Nsmi J ^-=~-2pq[E 1 p 7 q sin 2e l —E2p^q s (p i, —q i ) sin 2e—E%p(f sin 2e 3 



dt %n, 



o 



+E\p 5 q(pZ+3qZ) sin e\-E , pq{p^-q % f sin e'-£>g 5 (3p 3 +g a ) sin e' a 

J~J O hJ Li ul IX £i C I • . * • . . . . • . , . « . . I Jy I 



And by (21) 



3 



cos i ~ = —(iP—tf)[E l p* sin 2e 1 +EipY sin 2e+,% 8 sin 2e 3 



C% 8% 



+ E\2pY sin e\+E'2p*g*(p*-tf)* sin e'+^'s W sin e' 3 ] . (54) 



By (33) and (34), and remembering to take the halves of ^& mMf and $$ mm/ , and that 
sin 2=^, cos i=P 



XV Dill m o 7 / — — 



rr--j~ 



-^Q[PP() 3 sin 2e+P"P«^ sin e'] 



B»i 







(55) 



cos i£^=— 'PJ^G 1 sin 2e+p v P 2 (? 3 sin e' 






3 Q 



(56) 
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dP 

Now to obtain fi-j. we have to add the last four expressions together, and we 

dP . 
observe that the last two cut one another out, so that the expression for — is inde- 
pendent of the solar tides; also the terms in sin 2e, sin e cut one another out in the 

dP 
sum of the first two expressions, and hence it follows that -j- is independent of the 

Cvv 

sidereal semi-diurnal and diurnal terms. 
Thus we have 

^=— [^sii^-^^ (57) 

This equation will be referred to hereafter as that of tidal reaction.* From its form 
we see that the tides of speeds 2(n+/2), n+2/2, and 2/2 tend to make the moon 
approach the earth, whilst the other tides tend to make it recede. 

Then if, as in previous cases, we put JE 1 =JS 2 =JE; E\=E^E / ; e^e^e; 6 / 1 =€ / a =e / 
(which is justifiable so long as the moon's orbital motion is slow compared with that of 
the earth's rotation), we have, after noticing that 

p 8 — . q 8 = (p 2 — (f) (p 4 + g 4 ) = cos i( 1 — ^ sin 3 i) 
4p 6 g 2 -—4p 2 ^ 6 =4p 2 5 2 (j9 s -~g' 2 )=: sin 2 i cos i 
6p 4 g^=-| sin 4 i 

dP t^ 
fi—= — [cos i(l— \ sin 2 i)E sin 2e+ sin 2 i cos %E sin e'— f sin 4 %E f sin 2e"~] . (58) 

Now if the present values of n 9 12, i be substituted in this equation (58) (i.e., with 
the present day, month, and obliquity), and if the tropical year be the unit of time, it 
will be found that 

10 10 ^=~(24'27^sin 2e+4'18^' sin e'--27lE' sin 2c") 
at g 

£ 12 enters into this equation because r varies as fl 2 and therefore as f~ 6 . 

But we may here put £=1, because at present we only want the instantaneous rate 
of increase of 12. 

Now — = — \DT*£l£— -= — - ^-zr —r~ when f2=f2 ; hence multiplying the equation by 

Ctt Ctt OiiZr\ Ctv "* 

3/2 we have at the present time 



dO 

^10 10 ~=6115E sin 26+1053^ sin e'-68'28E" sin 2e" . . . (59) 
in radians per annum. 

* In a future paper on the perturbations of a satellite revolving about a viscous primary, I snail obtain 
this equation by the method of the disturbing function. 
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Then if for the moment we call the right-hand of this equation k, we have 

t 
S2=I2 —hT~~ . Integrating a second time, we find that the moon has fallen behind her 

proper place in her orbit -J# ~-~^ seconds of arc in the time t. Put t equal a 

century, and substitute for h 9 and it will then be found that the moon lags in a century 

63O7J0 sin 2e+108'6j£ / sin c' — 7-042J2" sin 2e f seconds of arc . . (60) 

But it was shown in Section 13 (48) that the moon, if unaffected by tidal reaction, 
would have been apparently accelerated 1043'3li/sin 2e+232*5i£ / s in e seconds of arc 
in a century. 

Hence taking the difference of these two, we find that there is an apparent 
acceleration of the moon's motion of 

412-6^ sin 2e+123-9£ / sin e'+7'042.#" sin 2e" ..... (61) 

seconds of arc in a century. 

Now according to Adams and Delaunay, there is at the present time an unex- 
plained acceleration of the moon's motion of about 4" in a century. For the present 
I will assume that the whole of this 4" is due to the bodily tidal friction and reaction, 
leaving nothing to be accounted for by ocean tidal friction and reaction, to which the 
whole has hitherto been attributed. Then we must have 

412-6^siii 2€+123-9^T sin e'+7-042^" sin 2c /7 =4 .... (62) 

This equation gives a relation which must subsist between the heights E, E\ E", of 
the semi-diurnal, diurnal, and fortnightly bodily tides, and their retardations e, e, e" \ 
in order that the observed amount of tidal friction may not be exceeded. But no 
further deduction can be made, without some assumption as to the nature of the 
matter constituting the earth, 

I shall first assume then that the matter is purely viscous, so that E= cos 2e, 

E'= cos e', E"= cos 2e", and tan 2e= — , tan e =-, tan 2e =— . The equation then 

becomes 

412*6 sin 4e+123'9 sin 2e / +7'042 sin 4e / =8 (63) 

If the values of c, e, e" be substituted, we get an equation of the sixth degree for p, 
but it will not be necessary to form this equation, because the question may be more 
jsimply treated by the following approximation. 



3 



Q Zi 
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There are obviously two solutions of the equation/ one of which represents that the 
earth is very nearly fluid, and the other that it is very nearly rigid. 

In the first case, that of approximate fluidity, e, e, e are very small, and therefore 

f2 4 

sin 4e=4e, sin 2e / =2e / =2e, sin 4e"~4e // =4— e=~— ~e 

TV & ( 'o/t 

Hence 

(1650 + 248 + 37^82 of 7-04)€=8 
whence 

That is to say, the semi-diurnal tide only lags by the small angle 14'. But this is 
not the solution which is interesting in the case of the earth, for we know that the 
earth does not behave approximately as a fluid body. 

In the other solution, 2e and e approach 90°, so that p is small ; hence 



4np p . ^ , 2npZ 2p -, ., . , „ 4flp 

~ — — -t— =9 sin 2e ==- =~ very nearly, and sin 4e —~-z — ~ 



Hence we have 

41 2 -6(e)- H 23- 9 (^) + 704 2 ^=8 

Put ~=x, so that x= cot 2e" : then substituting for — its value .^ — . we have 
2/2 & 7i 27-32' 

27-32 ' ^ 2 + l 

whence 

x 3 — '1655x 2 +l-2921tf-- 1655 = 

This equation has two imaginary roots, and one real one, viz. : '12858. Hence the 
desired solution is given by cot 2e // ^= , 12858 ; and 2e // =-|7r— 7°.20 / , and the corres- 
ponding values of 2e and e are 2€=j7r-— 16', and e'=i^r— 32'. If these values for e, 
c, c" be used in the original equation (63), they will be found to satisfy it very closely ; 
and it appears that there is a true retardation of the moon of 3 //# l in a century, whilst 
the lengthening of the day would make an apparent acceleration of 7' / 'l,— the difference 
of the two being the observed 4". 

With these values the semi-diurnal and diurnal ocean-tides are, according to the equi- 
librium theory of ocean-tides, sensibly the same as those on a rigid nucleus, whilst the 
fortnightly tide is reduced to sin 2e" or '992 of its theoretical amount; and the time 

it e" 

of high tide is accelerated by 77.— tt, or 6 -J hours in advance of its theoretical time/'" 



# In the abstract of this paper (Proc. Roy. Soc; ? No. 191, 1878) the height and lag of the bodily tid 
were accidentally given instead of the height and acceleration of the ocean tide. 
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If these values be substituted in the equation giving the rate of variation of the 
obliquity, it will be found that the obliquity must be decreasing at the rate of '00197° 
per million years, or 1° in 500 million years. Thus in 100 million years it would ouly 
decrease by 12'. So, also, it may be shown that the moon's sidereal period is being 
increased by 2 hours 20 minutes in 100 million years. 

Lastly, the earth considered as a clock is losing 13 seconds in a century. 

There is another supposition as to the physical constitution of the earth, which will 
lead to interesting results. 

If the earth be elastico- viscous, then for the semi-diurnal and diurnal tides it might 
behave nearly as though it were perfectly elastic, whilst for the fortnightly tide it 
might behave nearly as though it were perfectly viscous. With the law of elastico- 
viscosity used in my previous paper, * it is not possible to satisfy these conditions very 
exactly. But there is no reason to suppose that that law represents anything but an 
ideal sort of matter ; it is as likely that the degradation of elasticity immediately after 
straining is not so rapid as that Jaw supposes. I shall therefore take a limiting case, 
and suppose that, for the semi-diurnal and diurnal tides, the earth is perfectly elastic, 
whilst for the fortnightly one it is perfectly viscous. This hypothesis, of course, will 
give results in excess of what is rigorously possible, at least without a discontinuity in 
the law of degradation of elasticity. 

It is accordingly assumed that the semi-diurnal and diurnal bodily tides do not lag, 
and therefore €=c / =0 ; whilst the fortnightly tide does lag, and if = cos 2e\ 

Thus by (38) there is no tidal friction, and by (60) there is a true acceleration of the 
moon's motion of ^ of 7*042 sin4e // seconds of arc in a century. Then if we take the 
most favourable case, namely, when €"=22° 30', there is a true secular acceleration of 
3"*521 per century. 

It follows, therefore, that the whole of the observed secular acceleration of the moon 
might be explained by this hypothesis as to the physical constitution of the earth. 
On this hypothesis the fortnightly ocean tides should amount to sin 22° 30', or *38 
of its theoretical height on a rigid nucleus, and the time of high water should be 

accelerated by 1 day 17 hours. Again, by (35) — = — ^-t^smA*, from whence it may be 

shown that the obliquity of the ecliptic would be decreasing at the rate of 1° in 
128 million years. 

The conclusion to be drawn from all these calculations is that, at the present time, 
the bodily tides in the earth, except perhaps the fortnightly tide, must be exceedingly 
small in amount ; that it is utterly uncertain how much of the observed 4" of acce- 
leration of the moon's motion must be referred to the moon itself, and how much to 

# Namely, that if tlie solid be strained, the stress required to maintain it in the strained configuration 
diminishes in geometrical progression as the time, measured from the epoch of straining, increases in 
arithmetical progression. See Section 8 of the paper on " .Bodily Tides," &c, Phil, Trans., Part I., 1879, 
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the tidal friction, and accordingly that it is equally uncertain at what rate the day is 
at present being lengthened ; lastly, that if there is at present any change in the 
obliquity to the ecliptic, it must be very slowly decreasing. 

The result of this hypothesis of elastieo-viscosity appears to me so curious that I 
shall proceed to show what might possibly have been the state of things a very long 
time ago, if the earth had been perfectly elastic for the tides of short period, but 
viscous for the fortnightly tide. 

There will now be no tidal friction, and the length of day remains constant. The 
equation of tidal reaction reduces to 

£.„_* ^ 3 * 4 * * A // 
LL ~~~r~ —— "~™" 7~To ~t~q Sill v Sill 4 6 



clt 

Here v? is a constant, being the value of — at the epoch ; and ^ 2 4-^ 13 is the value 

of - — at the time t. 

The equation giving the rate of change of obliquity becomes 

cli it? 



, 12 yq sir ^ sm 4e 



clt 

Dividing the latter by the former, we have* 

sin idi— /xd£ 
And by integration 

cosi=ccm — /*(£— 1) 

If we look back long enough in time, we may find f=l-0l, and /x being 4*007, we 

have 

cos i=cos^ — * 04007 

Taking v=23° 28', we find ^=28° 40'. 

This result is independent of the degree of viscosity. When, however, we wish to 
find how long a period is requisite for this amount of change, some supposition as to 
viscosity is necessary. The time cannot be less than if sin4e // =l, or e'^22 30', and 
we may find a rough estimate of the time by writing the equation of tidal reaction 

^t 3 ^ * 4 T 



clt 
where I is constant and equal to 24°, suppose. Then integrating we have 



or 



*=-if^cosec 4 I(P-l). 



U 2 



* Concerning the legitimacy of this change of variable, see the following section, 
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When £=1 - 01, we find from this that — £ = 720 million years, and that the length 
of the month is 28*15 m. s. days. Hence, if we look back 700 million years or more, 
we might find the obliquity 28° 40', and the month 28*15 m. s. days, whilst the length 
of day might be nearly constant. It must, however, be reiterated, that on account 
of our assumptions the change of obliquity is greater than would be possible, whilst 
the time occupied by the change is too short. In any case, any change in this direction 
approaching this in magnitude seems excessively improbable. 



PART II. 

§ 15. Integration of the differential equations for secular changes in the variables in 

the case of viscosity/ 7 

It is now supposed that the earth is a purely viscous spheroid, and I shall proceed to 
find the changes which would occur in the obliquity to the ecliptic and the lengths of 
the day and month when very long periods of time are taken into consideration. 

I have been unable to find even an approximate general analytical solution of the 
problem, and have therefore worked the problem by a laborious arithmetical method, 
when the earth is supposed to have a particular degree of viscosity. 

The viscosity chosen is such that, with the present length of day, the semi-diurnal 
tide lags by 17° 30'. It was shown above that this viscosity makes the rate of change 
of obliquity nearly a maximum, t It does not follow that the whole series of changes 
will proceed with maximum velocity, yet this supposition will, I think, give a very 
good idea of the minimum time, and of the nature of the changes which may have 
occurred in the course of the development of the moon-earth system. 

The three semi-diurnal tides will be supposed to lag by the same amount and to be 
reduced in the same proportion ; as also will be the three diurnal tides. 

There are three simultaneous differential equations to be treated, viz. : those giving 
(1) the rate of change of the obliquity of the ecliptic, (2) the rate of alteration of the 
earth's diurnal rotation, (3) the rate of tidal reaction on the moon. They will be 
referred to hereafter as the equations of obliquity, of friction , and reaction respectively. 

To write these equations more conveniently a partly new notation is advantageous, 
as follows : — 

The suffix to any symbol denotes the initial value of the quantity in question. 

9 9 

if- " T T" T" 

Let « 3 =— , t«. 3 =— , 1111,=— ; these three quantities are constant. 

* This section has been rearranged, partly rewritten, and recomputed since the paper was presented. 
The alterations were made on December 19, 1878. 

f If I had to make the choice over again I should choose a slightly greater viscosity as being more 
interesting. 
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Since the tidal reaction on the sun is neglected, r / is a constant, and since r varies as 
/2 3 (and therefore as f"~ 6 ) ; hence 



t 3 % i$ t/ n rr / n Q uu t 



pb n £ 12 %n % 



y>r> 



% n n % 



8 



2gaio 



Let p be equal to 1Q , where v is the coefficient of viscosity of the earth. Then 
according to the theory developed in my paper on tides 4 '' 



tan 2e= — , tan e =-, tan 2e = — 
P P P 



• * * » 



. . . (64) 



To simplify the work, terms involving the fourth power of the sine of the obliquity 
will be neglected. 

Now let 



sin 3 i 



■^ 



P=i lo 8'io e > Q=l sin3 * lo gio e > R=-fcr7 lo gio e=iQ sec i 



cos % 



= % sur i log 10 e, V = 



i 



f COS" t 

1— I sin 3 i 



; loglO e 



W=i cos 2 i, X=-| sin 3 i cos i, Z=-| sin 3 i cos 3 i 



2 



. (65) 



* — 5 



1/2 

Also let sft ,f2 *=-, —=N; and it may be called to mind that £ 

2 



az> 



L\# 



(l+i>) 



The terms depending on the semi-annual tide will be omitted throughout. 
With this notation the equation of obliquity (35) and (36) may be written, 



cli 
di 



log 10 e ~ = sin i cos i(l — § sin 2 i) 



w 



13 



f-t^ 3 )(P sin 4e+Q sin 2e) 



' — '(Usin 4e+ V sin 2e) — -----It sin 4 



// 



The equation (43) of friction becomes 



1vlv> 



dN / V? \ /TTr . -x T . A WWy „ . 

"y- = (^+^ / )(W sm 4c+Xsin 2e)+^- Z sin 2e 



And by (58), Section 14, the equation of reaction becomes 



'W 



• » • * • * l VJ V7 1 



• (67) 



u-~ = -r7"(W sin 4e+X sin 2e) . . 



k • s 



(68) 



* Phil. Trans., 1879, Part I. 
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This is the third of the simultaneous differential equations which have to be treated. 
The four variables involved are i, N, £, t, which give the obliquity, the earth's rotation, 
the square root of the moon's distance and the time. Besides where they are involved 
explicitly, they enter implicitly in Q, R, U, V, W, X, Z, sin 4e, sin 2e', sin 4e". 

Q, R, &c, are functions of the obliquity i only, but P is a constant. Also 

a ^ n P 4:n pJSf . , 2n pJSf . „ . 4f2 0/ o£ 3 T , , 

sin 4e= . ' - =- — r-ri: :, sin 2e = -r-^ -, sin 4e = . _ * ft . I made several 

4^ 2 + p 2 An^N* + /> % ^ + p 4O 2 -f p 2 ^ 6 

attempts to solve these equations by retaining the time as independent variable, and 
substituting for £ and JV approximate values, but they were all unsatisfactory, because 
of the high powers of £ which occur, and no security could be felt that after a con- 
siderable time the solutions obtained did not differ a good deal from the true one. 
The results, however, were confirmatory of those given hereafter. 

The method finally adopted was to change the independent variable from t to c~. 
A new equation was thus formed between iV and £ which involved the obliquity i 
only in a subordinate degree, and which admitted of approximate integration. This 
equation is in fact that of conservation of moment of momentum, modified by the 
effects of the solar tidal friction. Afterwards the time and the obliquity were found 
by the method of quadratures. As, however, it was not safe to push this solution 
beyond a certain point, it was carried as far as seemed safe, and then a new set of 
equations were formed, in which the final values of the variables, as found from the 
previous integration, were used as the initial values. A similar operation was carried 
out a third and fourth time. The operations were thus divided into a series of periods, 
which will be referred to as periods of integration. As the error in the final values in 
any one period is carried on to the next period, the error tends to accumulate ; on this 
account the integration in the first and second periods was carried out with greater 
accuracy than would in general be necessary for a speculative inquiry like the present 
one. The first step is to form the approximate equation of conservation of moment of 
momentum above referred to. 

Let A=W sin 4e+X sin 2c', B=Z sin 2e'. 

Then the equations of friction (67) and reaction (68) may be written, 

-^f^+^A+^B (69) 



«o»*§=piA (70) 

We now have to consider the proposed change of variable from t to £ 

dN . dP 

The full expression for ~ contains a number of periodic terms ; ~ also contains 
terms which are co-periodic with those in — . Now the object which is here in view 

Civ 

MDCCCLXXIX. 3 R 
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is to determine the increase in the average value of N per unit increase of the average 

value of £ The proposed new independent variable is therefore not £, but it is the 

average value of £ ; but as no occasion wilharise for the use of $ as involving periodic 

terms, I shall retain the same symbol. 

In order to justify the procedure to be adopted, it is necessary to show that, if f(t) 

be a function of t, then the rate of increase of its average value estimated over a 

period T, of which the beginning is variable, is equal to the average rate of its increase 

estimated over the same period. Now the average value of f(t) estimated over the 

1 f t+T 
period T, beginning at the time t is - f{f)dt } and therefore the rate of the increase 

d 1 [ t+T . . 1 [ t+T , 

of the average value is -j — f(t)dt, ■ which, is equal to - f\t)dt) and this last 

expression is the average rate of increase of f(t) estimated over the same period. This 
therefore proves the proposition in question, 

Now suppose we have — = — M+ periodic terms, where M varies very slowly; 

Cvv 

then — M is the average value of the rate of increase of N estimated over a period 
which is the least common multiple of the periods of the several periodic terms. Hence 
by the above proposition — M is also the rate of increase of the average value of N 
estimated over the like period. 

dP 
Similarly if ~ = X + periodic terms, X is the rate of increase of the average value 

(XL 

of £ estimated over a period, which will be the same as in the former case. 

But the average value of N is the proposed new dependent variable, and the average 
value of $ the new independent variable. Hence, from the present point of view, 

— = — — . This argument is, however, only strictly applicable, supposing there are not 

dN dP 
periodic terms in — or -r of incommensurable periods, and supposing the periodic terms 

are rigorously circular functions, so that their amplitudes and frequencies are not func- 
tions of the time. 

It is obvious^ however, that if the incommensurable terms do not represent long 
inequalities, and if M and X vary slowly, then the theorem remains very nearly true. 
"With respect to the variability of amplitude and frequency, it is only necessary to pos- 
tulate that the so-called periodic terms are so nearly true circular functions that the 
integrals of them over any moderate multiple of their period is sensibly zero, to apply 
the argument. 

Suppose, for example, \fi(t) cos (^+x(0) were one of the periodic terms, then we have 
only to suppose that \fj(t) and ^(0 vary so slowly that they remain sensibly constant 

during a period — or any moderately small multiple of it, in order to be safe in 



2rr 

assuming 

Jo 



\jj(t) cos (vt-^xifydt as sensibly zero. Now in all the inequalities in iVand <f 
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it is a question of days or weeks, whilst in the variations of the amplitudes and 
frequencies of the inequalities it is a question of millions of years. Hence the above 
method is safely applicable here. 

It is worthy of remark that it has been nowhere assumed that the amplitudes of the 
periodic inequalities are small compared with the non-periodic parts of the expression. 

A precisely similar argument will be applicable to every case where occasion will 
arise to change the independent variable. The change will accordingly be carried out 
without further comment, it being always understood that both dependent and inde- 
pendent variable are the average values of the quantities for which their symbols 
would in general stand."* 

Then dividing (69) by (70) we have 

-^r 1+ ffi n+ l£y (71) 

Now — = — ^-~i = sin 2 i . approximately. This approximation will be suffi- 

sin 2e' 
ciently accurate, because the last term is small and is diminishing. For the same 
reason, only a small error will be incurred by treating it as constant, provided the 
integration be not carried over too large a field — a condition satisfied by the proposed 
" periods of integration." Attribute then to % e, e average values, and put 



»=+£)' r=%*n%£ (72) 



and integrate. Then we have 

i\r=l+^{(W)+/3(l-P)+ r (l-H} (73) 

This is the approximate form of the equation of conservation of moment of momentum, 
and it is very nearly accurate, provided £ does not vary too widely. 

By putting /3=0, y=0, we see that the equation is independent of the obliquity, if 
there be only two bodies, the earth and moon, provided we neglect the fourth power of 
the sine of the obliquity. 

The equation of reaction (68) may be written 

--r—^-f.— (Wsin4€+Xsin2€ / ) (7i) 

* In order to feel complete confidence in my view, I placed the question before Mr. E. J. Routh, and 
with great kindness he sent me some remarks on the subject, in which he confirmed the correctness of 
my procedure, although he arrived at the conclusion from rather a different point of view* 

3 r 2 
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dt 

Also, multiplying the equation of obliquity (66) by -77, we liave 



log 10 e di 1 dt 



to 2 



12 



+u y 2 ) (P sin 4e+Q sin 2e) 



sin i cos i (1 — | sin 3 i) c££ -A^f 

^(U sin 4e+ V sin 2e') — ~ R sin 4e 



WW, ,~r-r • . . -s— r • ~ /\ t-6 -.-„ • , ff 



dP . 

Now by far the most important term in ~r is that in which W occurs, and therefore 

^t~7, only depends on the obliquity in its smaller term. Then, since 2W=cos 2 £, 
therefore 

(JjTj X / n L -r-r-r wv 

— —- 2 W — 

d% cos 3 ^\ dg, 

Also 

cos 3 i , ._ ^ 1 

sin i cos i (1 — f sin 3 ^) ' ^ \/\ — | sin 3 i 

= <$ . log,, tan i (1 — \ sin 3 i) 

when the fourth power of sin { is neglected. 
Hence the equation may be written 



cos 3 i -».-,-, sin i 



%i? 



^log 10 tani(l— isin 3 i)=-( 2W— \ I ~^+u? )( P sin 4e +Q sin 2e 



tV -»~^ • . »/ tt'tt' 



— ^Rsin4e' / ~^(IIsin4€+Vsin2e / ) (75) 

Now the term in P (which is a constant) is by far the most important of those 

dt 
within brackets [ ] on the right-hand side, and 2W— has been shown only to involve 

i in its smaller term. Hence the whole of the right-hand side only involves the 
obliquity to a subordinate degree, and, in as far as it does so, an average value may 
be assigned to i without producing much error. 

In the equation of tidal reaction (68) or (74) also, I attribute to i in W and X an 
average value, and treat them as constants. As the accumulation of the error of 
time from period to period is unimportant, this method of approximation will give quite 
good enough results. 

We are now in a position to track the changes in the obliquity, the day, and the 
month, and to find the time occupied by the changes by the method of quadratures. 

First estimate an average value of i and compute Q, R . . . Z, /3, y. Take seven 
values of £ viz. : 1, '98, '96 ... '88, and calculate seven corresponding values of N\ 
then calculate seven corresponding values of sin 4e, sin 2e', sin 4e". Substitute these 

values in — , and reciprocate so as to get seven equidistant values of — . 
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Combine these seven values by Weddle's rule, viz. : 



\ 



u x dx = £jh[u Q + u 2 + u B + 1%+ u Q + 5 (uy + u 3 + u B ) ] 



and so find the time corresponding to £="88. It must be noted that the time is 

negative because dg is negative. 

clt 
In the course of the work the values of —corresponding to £=1, # 96, '92, *88 

tig 

have been obtained. Multiply them by 2W ; these values, together with the four 

values of sin 4e, sin 2c', sin 4c" and the four of N, enable us to compute four of 

d 
— log 10 tani(l — \ sin 2 i), as given in (75). 

Combine these four values by the rule 



u x dx = — -[u + u s + 3 ( w x + 1%) ] 

and we get 

tan^(l~|-sm 3 ^) 

^ 10 tan i {l — J sin 3 i ) 

from which the value of i corresponding to f=*88 may easily be found. It is here 
useless to calculate more than four values, because the function to be integrated does 
not vary rapidly. 

We have now got final values of i 9 N, t corresponding to £=*88. 

Since the earth is supposed to be viscous throughout the changes, therefore its 
figure must always be one of equilibrium, and its ellipticity of figure e=^N 2 e . 

Also since £= f -~ J = a/ -, where c is the moon's distance from the earth, therefore 
„ = p/_o \ ^ich gives the moon's distance in earth's mean radii. 

The fifth and sixth column of Table IV. were calculated from these formulas. 

The seventh column of Table IV. shows the distribution of moment of momentum in 
the system; it gives /x the ratio of the moment of momentum of the moon's and earth's 
motion round their common centre of inertia to that of the earth's rotation round its 
axis, at the beginning of each period of integration. 

Table I. shows the values of e, e', c" the angles of lagging of the semi-diurnal, 
diurnal, and fortnightly tides at the beginning of each period. 

Tables II. and III. show the relative importance of the contributions of each term 

clt d 

to the values of ~ and — log 10 tani(l — -|shrH) at the beginning of each period. 

The several lines of the Tables II. and III. are not comparable with one another, 
because they are referred to different initial values of fl and n in each line. 

I will now give some details of the numerical results of each integration. The 
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computation as originally carried out* was based on a method slightly different 
from that above explained, but I was able to adapt the old computation to the 
above method by the omission of certain terms and the application of certain 
correcting factors. For this reason the results in the first three tables are only given 
in round numbers. In the fourth table the length of day is given to the nearest 
five minutes, and the obliquity to the nearest five minutes of arc. 

The integration begins when the length of the sidereal day is 23 hrs. 56 min., 
the moon's sidereal period 27*3217 m. s. days, the obliquity of the ecliptic 23° 28', 
and the time zero. 

First period. — Integration from £=1 to '88 ; seven equidistant values computed 
for finding the time, and four for the obliquity. 

For the obliquity the integration was not carried out exactly as above explained, in 

cl cl 

as far as that — log 10 tan i was found instead of — log 10 tan i(l — -|sin%'), but the differ- 

cl 

ence in method is very unimportant. The result marked'* in Table III. is — log 10 tan i. 

The estimated average value of i was 22° 15'. 
The final result is 

iV=r550, ^=20° 42', — £ = 46,301,000 

Second period. — Integration from £=l to '76 ; seven values computed for the time, 
and four for the obliquity. 

The estimated average for i was 19°. 
The final result 

#=1-559, i=l7°2l\ -£=10,275,000 

Third period. — -Integration from. f=l to 76 ; four values computed. 
The estimated average for i was 16° 30'. 
The final result 

#=1-267, i=15° 30', —£=326,000 

Fourth period. — Integration from £=1 to '76 : four values computed. 

The estimated average for i was 1 5°. The small terms in /3 and y were omitted in 
the equation of conservation of moment of momentum. All the solar and combined 
terms, except that in V in the equation of obliquity, were omitted. 

The final result 

#=1-160, i=14° 25 r , -£=10,300 



* I have to thank Mr. E. M. Langley, of Trinity College, for carrying out the laborious computations. 
The work was checked throughout by myself. 
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Table I. — Showing the lagging of the several tides at the beginning of each period. 





Semi-diurnal 


Diurnal 

(O- 


Fortnightly 


I. 


17F 


19i° 


0° 44' 


II. 


m° 


281° 


1° 5' 


III. 


29| 


40° 


2° 27' 


IV. 


321° 


46|° 


5° 30' 



Table II. — Showing the contribution of the several tidal effects to tidal reaction 



i.e., to -~ ) at the beginning of each period. The numbers to be divided by 10 1 



o 





Semi-diurnal. 


Diurnal. 


I. 


12- 


1-2 


II. 


69- 


6-3 


III. 


2200- 


200- 


IV. 


70000- 


6100- 



Table III. — Showing the contributions of the several tidal effects to the change of 



d 



obliquity (i.e., to j- log 10 tan i(l—jf sin 2 i)) at the beginning of each period. 

tig J * ' ""* 



*I. 


Lunar 
semi-diurnal. 


Lunar 
diurnal. 


Solar 
semi-diurnal. 


Solar 
diurnal. 


Combined 
semi-diurnal. 


Combined 
diurnal. 


Fortnightly. 


d 

-7 log tan i (1 — • | sin 2 i) . 
tit, 


•82 


•13 


•18 


•03 


-•06 


-•48 


-•006 


•60* 


II. 


•44 


•06 


•02 


• • 


-•01 


-•16 


-•003 


•36 


III. 


•22 


•03 


• • 


• • 


• • 


-•02 


-•003 


•23 


IV. 


•13 


•02 


• • 


• • 


• • 


• • 


-•004 


14 
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Table IV. — Showing the physical meaning of the results of the integration. 





Time 

(-0- 


Sidereal 

day in m. s. 

hours. 

h. m. 

23 56 


Moon's side- 
real period 
in m. s. days. 


Obliquity of 
ecliptic 

(0. 


Reciprocal 

of ellipticity 

of figure. 


Moon's 

distance in 

earth's mean 

radii. 


Ratio of m. 
of m. of orbi- 
tal motion to 
m. of m. of 
earth's 
rotation. 


Heat gene- 
rated (i?ee 
Section 16). 


Initial 
state. 


Years. 




d. 

27-32 


23° 28' 


232 


60-4 


4-01 


Degrees Fahr. 

0° 

225° 


I. 


46,300,000 


15 30 


18-62 


20° 40' 


96 


46*8 


2-28 


II. 


56,600,000 


9 55 


817 


17° 20' 


40 


27-0 


1-11 


760° 


III. 


56,800,000 


7 50 


3-59 


15° 30'* 


25 


15-6 


-67 


1300° 


IV. 


56,810,000 


6 45 


1-58 


14° 25'* 


18 


9-0 


•44 


1760° 



The whole of these results are based on the supposition that the plane of the lunar 
orbit will remain very nearly coincident with the ecliptic throughout these changes. 
I now (July, 1879), however, see reason to believe that the secular changes in the 
plane of the lunar orbit will have an important influence on the obliquity of the 
ecliptic. Up to the end of the second period the change of obliquity as given in 
Table IV. will be approximately correct, but I find that during the third and fourth 
periods of integration there will be a phase of considerable nutation. The results in 
the column of obliquity marked (*) have not, therefore, very much value as far as 
regards the explanation of the obliquity of the ecliptic ; they are, however, retained as 
being instructive from a dynamical point of view. 



§ 16. The loss of energy of the system. 

It is obvious that as there is tidal friction the moon-earth system must be losing 
energy, and I shall now examine how much of this lost energy turns into heat in the 
interior of the earth. The expressions potential and kinetic energy will be abbreviated 
by writing them p.e. and he. 

The he. of the earth's rotation is \MaV. 

The he. of the earth's and moon's orbital motion round their common centre of 
inertia is 

1 TiJll mr Y^9 ,1 / MT Y^9 1 HJI 9 ^ 



2 



m + M, 



m 4- Mi 



2 



l-\-v 



(' fl\ ^1 -4- V Xz^T CJCfj 

- ) , so that — — = — . Hence the 
TV 14- z/ 



vv 



whole he. of the moon-earth system is 
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MaHhh*' 



i£± 



2 



v r 



The p. e. of tlie system is 



Mm 



M gd 
v r 



3 



Therefore the whole energy E of the system is 



Ma^ 



\rft 



1 I 
2v r 



and in gravitation units 



E—Ma 



\ 



5 



rfta 



1 <x 



# 2z> r 



Now since the earth is supposed to be plastic throughout all these changes, there- 
fore its ellipticity of figure 



e=% 



rfia 



and 



9 



E=MaL\e-~~ - 
[ *° 2v r 



If 6, e+Ae and r, r+Ar be the ellipticity of figure, and the moon's distance at 
two epochs, if J be Joule's equivalent, and or the specific heat of the matter con- 
stituting the earth ; then the loss of energy of the system between these two epochs 
is sufficient to heat unit mass of the matter constituting the earth 



Ma\ 



a. la 

J * 2v r 



degrees, 



and is therefore enough to heat the whole mass of the earth 



-j^|f 5 -Ae-^A-| degrees. 



It must be observed that in this formula the whole loss of he. of the earth's 
rotation, due both to solar and lunar tidal friction, is included, whilst only the gain of 
the moon's p.e. is included, and the effect of the solar tidal reaction in giving the 
earth greater potential energy relatively to the sun is neglected. 

In the fifth and sixth columns of Table IY. of the last section the ellipticity of figure 
and the moon's distance in earth's radii are given ; and these numbers were used in 
calculating the eighth column of the same table. 

I used British units, so that 772 foot-pounds being required to heat 1 lb. of water 
1° Fahr., J= 772; the specific heat of the earth was taken as -Jth, which is about that of 
iron, many of the other metals having a still smaller specific heat ; the earth's radius was 

MDCCCLXXIX. 3 S 
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taken, as before, equal to 20*9 million feet. Then the last column states that energy- 
enough has been turned into heat in the interior of the earth to warm its whole mass 
so many degrees Fahrenheit within the times given in the first column of the same 
table. 

The consideration of the distribution of the generation of heat and the distortion of 
the interior of the earth must be postponed to a future occasion. 

In the succeeding paper I have considered the bearing of these results on the secular 
cooling of the earth, and in a subsequent paper ( e Proceedings of the Royal Society/ 
No. 197, June 19, 1879, p. 168) the general problem of tidal friction is considered by 
the aid of the theory of energy. 

§ 17. Integration in the case of small variable viscosity * 

In the solution of the problem which has just been given, where the viscosity is 
constant, the obliquity of the ecliptic does not diminish as fast as it might do as we 
look backwards. The reason of this is that the ratio of the negative terms to the 
positive ones in the equation of obliquity is not as small as it might be ; that ratio 

principally depends on the fraction . - , which has its smallest value when e is very 

Sill 4t6 

small. 

I shall now, therefore, consider the case where the viscosity is small, and where it 
so varies that e always remains small. 

This kind of change of viscosity is in general accordance with what one may 
suppose to have been the case, if the earth was a cooling body, gradually freezing as 
it cooled. 

The preceding solution is moreover somewhat unsatisfactory, inasmuch as the three 
semi-diurnal tides are throughout supposed to suffer the same retardation, as also are 
the three diurnal tides ; and this approximation ceases to be sufficiently accurate 
towards the end of the integration. 

In the present solution the retardations of all the lunar tides will be kept distinct. 

By (40) and (40'), Section 11, 

2(n—f2) 2n 2(n + fl) „ 212 

tan 2€,= — — , tan 2e=— , tan 2e«= — , tan 2e = — , 

P P P P 

. n—2f2 . f n , n-\-2I2 
tanei = , tan e =-, tane 9 = 

1 p P P 

for the lunar tides. 

For the solar tides we may safely neglect il i compared with n, and we have 

* This section has been partly rewritten and rearranged, and wholly recomputed since the paper was 
presented. The alterations are in the main dated December 19, 1878. 
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tan2e=— -, tan e'=- for the semi-diurnal and diurnal tides respectively. The semi- 
P P 



annual tide will be neglected. 



n 



Then if the viscosity so varies that all the es are always small, and if we put — =X, 



we have 



sin 4e T 
sin4e 

sin 2e' 



1- 



sin 4< 



i_i 



=1— X 



- Dill TtCft 

X, — ~- 
sin 2e' 



1+X, 



sln4e 



1 

2> 



sin 4e" 
sin 4e 

sin JjG i 



=x 



2. 



sin 4e s 



l+X 



• * * » 



(76) 



By means of these equations we may expresss all the sines of the e's in terms of 
sin 4e. 

Then, remembering that the spheroid is viscous, and that therefore Ey= cos 2e 
E\= cos e' l9 &c., we have by Sections 4 and 7, equations (16) and (29), 



i? 



di 



vtv 



1 T 3 



dt N Qn 



p^q sin 4€ 1 —p B q s (p*-~q 2 ) sin 4e— |pg 7 sin 4e 2 — f p 3 $ 3 sin 4c' 7 



i/n7 
.2 



fi^(^ + 32 3 )sin2€' 1 -^(p 3 -^) 3 sin2e / -lp 2 5 (3p 3 + 2 3 )sin2 e ' 3 ] . (77) 



diK 



m 



a T 



3 



^ fi^i 



— LiP 8 sin 4e 1 +%>V sin 4e+-|2 8 sin 4e 3 

8% 



+pV sin 2^ +pV(p*-g a ) s sin 2€'+# Y sin 2 € ' a ] . (78) 



And by (57), Section 14, 

d>l* t 3 
/ji-~ = — [^p B sin 4c!— ^q s sin 4e s — Sp 4 ^ sin 4e // +2p 6 g' 3 sin 2e\— 2p 2 g 6 sin 26^] • (79) 

The first two of these equations only refer to the action of the moon on the lunar 
tides, but the last is the same whether there be solar tides or not. 

Then if we substitute from (76) for all the e's in terms of sin 4c, and introduce 
cos i==P==p s — ■ g' 2 , sin i=Q=2pq, we find on reduction 



T * 



m* 



1 T S 



dt JV g% 



C-v-X V * 



.2 



w 



IV i 



sin4€[|P$— |X$] 
sin 4c[l — \Q % — XP] 



~N 



M*"~tT~~ "s sin 4fc€ jl ■"■ *" A J 

O S -<y 



^ 2 8^i 







• • • * • »\OU) 
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di dN 

The parts of •— and — which, arise from the attraction of the sun on the solar tides 

dt dt 

fl 

may be at once written down by symmetry, and X y =— may be considered as a small 

fraction to be neglected compared with unity. Thus we have 



^=4 ^ sin 4e J-P$ 
dt N%n Q 4 



dt 



=|^ S in4e(l-^) 



(81) 



Lastly as to the terms due to the combined action of the two disturbing bodies, it 
was remarked that they only involved e and e', which are independent of the orbital 
motions. 

Thus by (33) we have 



~dT' 



1 TT 



iVg^ 



- sin ie.^PQ 



dK 



mm/ 



rrr f 



dt 



~n, 



' sin 4e.i() 3 



(82) 



~J 



Then collecting results from the last three sets of equations and substituting cos i 



n 



and sin i for P and 0, arid — for X, we have 

n 



1 sin4e 



[JL 



dt JSF pi Q 

^iV__ 1 sin4e 
' dt ' 

d± 

dt' 



^ sin % cos i 



T "-f-7", 



2, 



2/2 



TT, 



t 3 sec i 



n 



2 



'7b, 







(1 — J sin 2 i)(r s +T / 2 )+ i^tt, sin 2 i— ■— r 3 cos i 



/2 

i - — / 

71 



, sin4e . 



8^. 







n 



sec & 



. (83) 



These are the simultaneous equations which are to be solved. 

Subject to the special hypothesis regarding the relationship between the retarda- 



/jS L 



7h 



V s sec i in the 



tions of the several tides, and except for the neglect of a term • 
first of them, and of — ~rf cos % in the second, they are rigorously true. 

7b 

We will first change the independent variable in the first two equations from t to £ 
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Dividing the first and second equations by the third, and observing that 



we have 



— log tan 3 - 



fldf; 



2^ 7 1 x 2 % 

- — = d lo g tan^ - 



i+l?)' -®-T 

Nil sec % 

n 



sec£ 



l-—\ sin^i 






cos * 



l+(^ 






4-I-] sin i tan i— 



2 



/2 



/2 



sec£ 



n 



(84) 



If there be only one disturbing body, which is an interesting case from a theoretical 
point of view, the equations may be found by putting t=Q, and may then be written 



log tan^ - 



fldl; 






1 

'N 



cos i— 



2/2 



w 



cos^' 



n 

n 



l oir^a 



. n 



dN_ 
fidf 



2 



suri 



n 



cos& 



. /2 

C0S£ 

n 



/^"T^i sin 4e. — cos & 

t# * g?i \ n 



-J 



• • • • 



• » 



• (85) 



From these equations we see that so long as f2 is less than n cos i, the satellite 
recedes from the planet as the time increases, and the planet's rotation diminishes, 

because the numerator of the second equation maybe written cos if cos i )+-|sin 3 ^ 

which is essentially positive so long as fl is less than n cos i. But the tidal friction 

1 JL (3f)S % 1 -4- POS % 

vanishes whenever f2=n— -. The fraction —r - is however necessarily greater 

2 cos % 2 cos % J ° 

than unity, and therefore the tidal friction cannot vanish, unless the month be as 
short or shorter than the day. The obliquity increases if 12 be less than \n cos i, 
but diminishes if it be greater than \n cos i. Hence the equation 12= ^n cos i gives 
the relationship which determines the position and configuration of the system for 
instantaneous dynamical stability with regard to the obliquity (compare the figures 
2, 3, 4, Plate 36). From this it follows that the position of zero obliquity is one of 
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dynamical stability for all values of n between f2 and 2f2, but if n be greater than 2/2, 
this position is unstable. * 

n 

We will now return to the problem regarding the earth. We may here regard — as 



a small fraction, and i as sufficiently small to permit us to neglect ■§■ sin 4 i ; also I — sec i ) , 

— — sec L [-A — sec i will be neglected. 

r n \t J n 

* Added on September 25, 1879. — The result in the text applies to the case of evanescent viscosity. 
If the viscosity be infinitely large the sines of twice the angles of lagging will be inversely instead of 
directly proportional to the speeds of the corresponding tides (compare p. 482). Thus we must here invert 
the right-hand sides of the six equations (76). If the obliquity be very small (77), (78), (79) become 

£! = — i s in&'sin4ei 1 + -4 — s -~— 2(1— X) 

dt N g# L 1—2X v 'J 

1 T 9 ! . . . , /1 + 2X- 4X 3 \ i , ft „, N 

dN dP t 3 i . . 
— = w-s. = — i sin 4c-, 



When 2X=1, — - apparently becomes infinite; but in this case the viscosity must be infinitely large in 

order to make the tide of speed n—2Q lag at all, and if it be infinitely large sin 4e t is infinitely small. If 
the viscosity be large but finite, then when 2X=1, the slow diurnal tide of speed n— 20 is no longer a 

true tide, but is a permanent alteration of figure of the spheroid. Thus e' 1 =0 and — depends on 
[sii^ej— sin2 6 '] which is equal to sin4e 1 [l— 2(1 —X)] when the viscosity is large, and vanishes when 

•t * 

2X=1. Thus when the viscosity is very large (not infinite) -- vanishes when 2Q-~w=l, as it does when 

the viscosity is very small. 

a/ k _l \ d% 

When 1 -h2X— 4X 3 =0, that is, when X= — - — =1 — 1*236, — vanishes; and it is negative if X be a little 

greater, and positive if a little less than 14-1*236. And 1 — 2X is negative if X be greater than f. 

Hence it follows that for large viscosity of the planet, zero obliquity is dynamically unstable, if the satellite's 

period be less than 1*236 of the planet's period of rotation; is stable if the satellite's period be between 1*236 

and 2 of the planet' s period ; and is unstable for longer periods of the satellite. 

N d i l-f-2X— 4X 3 
If the viscosity be very large — -—log tan 3 -==--- - — ■, but if the viscosity be very small the same 

p> ag & 1 — 2X 

i 2\ 

expression =- — --. For positive values of X, less than 1 and greater than -6910 or l-~ 1*447, the former 

is less than the latter, and if X be less than 1 -—1*447 and greater than the former is greater than the 
latter. 

Hence if there be only a single satellite, as soon as the month is longer than two days, the obliquity of 
the planet's axis to the plane of the satellite's orbit will increase more, in the course of evolution, for 
large than for small viscosities. This result is reversed if there be two satellites, as we see by comparing" 
figs. 2 and 4, Plate 36, 
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Then our equations are 



log, tan 2 - 



1 + 



2 



2/2 



sec ^ 



n 



fid% 



2 



N 1 



n 



sec^ 



n 



— 7L=1 + (~I +i-sini tan^+ — (seci— 1) 



(86) 



The experience of the preceding integration shows that i varies very slowly com- 
pared with the other variables N and £; hence in integrating these equations an 
average value will be attributed to i, as it occurs in small terms on the right-hand 
sides of these equations. 

The second equation will be considered first. 

t It \ t 

We have t=-~, so that if we put /3=-i- 3 -( -M , y=-~^~ sin i tan % and omit the last 

S \ T 0/ T o 

term, we get by integrating from 1 to N and from 1 to £ 



i\r = i + ^i_£ +/3(1 _p )+r(1 _£7 )} 



(87) 



as a first approximation. This is the form which was used in the previous solution, 
for, by classifying the tides in three groups as regards retardation of phase, we virtually 
neglected ft compared with n. 



n 



This equation will be sufficiently accurate so long as — is a moderately small frac- 

tion ; but we may obtain a second approximation by taking account of the last term. 

Now 

1 



— (sec i— 1 )=i sir % — - • -r—. 



=^ sin 2 i 



very nearly 



n 


1 






1 + fjb 


-f 



by substituting an approximate value for N. 

A more correct form for the equation of conservation of moment of momentum will 
be given by adding to the right-hand side of equation (87) the integral of this last 
expression from 1 to £ and multiplying it by /jl. And in effecting this integration i 
may be regarded as constant. 
14- />& 



Let &=■ 



V* 



Then since 



L . J_ . JL i x 

7„ $-3 I 7„2 *-9 I 7„R f-"T" 



1 



f O-f) &f ' Pf 3 ' #£ ' /fc 3 (&-£) 
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Therefore 



= _j^A_ 1 \A + l+!M + /^\ 3 lo /i+Ki-?) 



Hence the second approximation is 

^l+ ( ,{(l-^^(l-f)+ r (l-a)+| S in^»^(---l|+^ 



H sin 3 1 -° f-^- Ylog 



1 + /*(!-?)■ 



(88) 



It would no doubt be possible to substitute this approximate value of N in terms of 
£ in the equation which gives the rate of change of obliquity, and then to find an 
approximate analytical integral of the first equation. But the integral would be very 
long and complicated, and I prefer to determine the amount of change of obliquity by 
the method of quadratures. 

In the present case it is obviously useless to try to obtain the time occupied by 
the changes, without making some hypothesis with regard to the law governing the 
variations of viscosity ; and even supposing the viscosity small but constant during 
the integration, the time would vary inversely as the coefficient of viscosity, and would 
thus be arbitrary. The only thing which can be asserted is that if the viscosity be 
small, the changes proceed more slowly than in the case which has been already solved 
numerically. 

To return, then, to the proposed integration by quadratures : by means of the 
equation (88) we may compute four values of N (corresponding, say, to £=1, '96, # 92, 

*88); and since r=-£, and — = — r— ^, we may compute four equidistant values of all 

the terms on the right-hand side of the first of equations (86), except in as far as i is 
involved. Now i being only involved in small terms, we may take as an approximate 
final value of i that which is given by the solution of Section 15, and take as the four 

q «_ q (q O \ 

corresponding values i Q9 v+^T^ ^o"^ — 7 )~ L > ^ 

Hence four equidistant values of the right-hand side may be computed, and com- 

rU 372, 

bined by the rule Uxdx=— [^ +^3+ 3 (^i+%)l which will give the integral of the 

* a 

q n 

right-hand side from <f to 1 ; and this is equal to log tan 3 ~— log tan 2 ^. 

The integration was divided into a number of periods, just as in the solution of 
Section 15. The following were the results : 
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First period. From f=l to -88; /x=4'0074; i=20° 28'; #= 1-5478. The term 
in — in the expression for JV added '0012 to the value of JV. 



™0 



Second period. From £=1 to 76; i u,=2 , 2784; i=l7° 4'; iV=l'5590. The term 

SI 

in — ° added '0011 to the value of N. 

n 

Third period. From f=l to 76; /x=11107; {=15° 22'; i^=l'2677. The term 

•T2 

in — added *0007 to the value of N. 

ft 

reaches its minimum about the end of the period, During the rest of the integration 
it increases. If we neglect the solar action and the obliquity, it is easy to find the 

minimum value of—. For — = — 9 —~z and reaches its minimum when — = — \ but 

n n n N% 6 cl% 



— = — p. Therefore 2V=^-. Now N=1-\-[l(1 — f) 5 an d hence f=f . If /*=4, 

fz=i|-=: , 9375. This value of £ is passed through at near the end of the first period 
of integration. At this period there are 19*2 mean solar hours in the day; 22-| mean 
solar days in the sidereal month; and 28y rotations of the earth in the sidereal 
month. This result of 28y is, of course, only approximate, the true result being 
about 29.* 

The physical meaning of these results is given in a table below. 

At the end of the third period of integration the solar terms (those in — ) have 

become small in all the equations, and as they are rapidly diminishing they may be 
safely neglected. To continue the integration from this point a slight variation of 
method will be convenient. 

Our equations may now be written approximately 

m 

.. . _j i- ■"■"* sec (/ 

—Mr io s tan 2 = n 7Z^~ 

n 

In order to find how large a diminution of obliquity is possible if the integration 
be continued, we require to stop at the point where n cos i=2S2. 
Now the equation N=l+fi(l — £) may be written 



n 



n 



M i -^\ 



* The subject is referred to from a more general point of view in a paper on the " Secular Effects of 
Tidal Friction," see ' Proc. Roy. Soc.,' No. 197, 1879. 
MDCCCLXXIX. 8 T 
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If therefore we put x=yi2, we must stop the integration at the point where 
n=2x 2, sec i, x being given by the equation 



r sec % 

- = 1+jU, 



n 



•t \Uh 



And if we assume i=14°, x is given by 



^ 4 ""i n o cos 14°(1+/a)$+ — cos 14°=0 

because /x= 1 -r-sw /2 i 

Now at the end of the third period of integration, which is the beginning of the 
new period, I found 

log w =3'84753, log /a=9'82338-10, and log $=5'39378 — 10 

The unit of time being the present tropical year. 

Hence the equation is 

a?*— 5690a;+19586 = 

The required root is nearly ^/5690, and a second approximation gives x=I2 l — 16703 
(16*51 would have been more accurate). 

But if2 *=8*616. Hence we desire to stop the integration when 

g ~\nJ ~ 16-703 - 

Now /x=-6659 ; hence when |=*516, iV"= 1-322. 

In order to integrate the equation of obliquity by quadratures, I assume the four 
equidistant values, 

iV=l-000, 1-107, 1-214, 1*321 

TV"— 1 
And by means of the equation £=1 — -— — = 1 — (iV— l)(r502) the corresponding 

vj vJ tj tj 

values of £ are found to be 

1-000, -8393, -6786, -5179 

Then by means of the formula— ==— — *, the corresponding values of — are found 

to be 

'0909, '1388, '2395, '4951 

I assumed conjecturally four values of i lying between i = 15° 22' and i=14°, which 
I knew would be very nearly the final value of i ; and then computed four equidistant 

values of — — log 10 tan -. 

The values were 

•19381, '16230, '11882, —'00684. 

The fact that the last value is negative shows that the integration is carried a little 
beyond the point when n cos i=2fl, but this is unimportant. 
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Combining these values by the rules of the calculus of finite differences, I find 
{=13° 59'. 

This final value of f (viz.: '5179) makes the moon's sidereal period 12 hours, and 
the value of N (viz.: 1*321) makes the day 5 hours 55 minutes. 

These results complete ijie integration of the fifth period. 

The physical meaning of the results for all five periods is given in the following 
table : — 



Sidereal day in m.s. 
hours and minutes. 


Moon's sidereal period 
in m.s. days. 


Obliquity of 
ecliptic. 


h. m. 

Initial 23 56 

15 28 

9 55 

7 49 
Final 5 55 


27*32 days 
18-62 „ 
o* ±7 ,, 

o*oy 5» 

12 hours 


23° 28' 
20° 28' 
17° 4' 
15° 22' * 
14° 0' * 



It is worthy of notice that at the end of the first period there were 28*9 days of 
that time in the then sidereal month ; whilst at the end of the second period there 
were only 19*7. It seems then that at the present time tidal friction has, in a sense, 
done more than half its work, and that the number of days in the month has passed 
its maximum on its way towards the state of things in which the day and month are 
of equal length — as investigated in the following section. 

In the last column of the preceding table the last two results in the column giving 
the obliquity of the ecliptic (which are marked with asterisks) cannot safely be 
accepted, because, as I have reason to believe, the simultaneous changes of inclination 
of the lunar orbit will, after the end of the second period of integration, have begun 
to influence the results perceptibly. 

For this same reason the integration, which has been carried to the critical point 

where n cos {=2/2, and where —changes sign, will not be pursued any further. Never- 

Cvv 

theless we shall be able to trace the moon's periodic time, and the length of day to 
their initial condition. It is obvious that as long as n is greater than /2, there will 
be tidal friction, and n will continue to approach /2, whilst both increase retrospectively 
in magnitude. 

I shall now refer to a critical phase in the relationship between n and /2, of a totally 
different character from the preceding one, and which must occur at a point a little 
more remote in time than that at which the above integration stops. 

This critical phase occurs when the free nutation of the oblate spheroid has a fre- 
quency equal to that of the forced fortnightly nutation. 

In the ordinary theory of the precession and nutation of a rigid oblate spheroid, the 
fortnightly nutation arises out of terms in the couples acting about a pair of axes 
fixed in the equator, which have speeds n — 2/2 and w+2/2. If C and A be the 
greatest and least principal moments of inertia, then on integration these terms are 

3 T 2 
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divided by — ~— n -{-n' = f2f2 and give rise to terms in ~ and ~~- sin i of speed 2/2, When 

2/2 is neglected compared with n 9 we obtain the formula, given in any work on physical 

astronomy, for the fortnightly nutation. 

C— A 

Now it is obvious that if — - — n-\-n=2f2, the former of these two terms becomes 

C— -A 
infinite. Since in our case the spheroid in homogeneous — - — =e the ellipticity of the 



A 



f? 



spheroid ; and since the spheroid is viscous e=|— . Therefore the critical relationship 

8 



7V> 



is A- — [-71=212. 
2 8 

When this condition is satisfied the ordinary solution is nugatory, and the true 

solution represents a nutation the amplitude of which increases with the time. 

2/2 
The critical point where the above integration stops is given by — =cos i, and this 

7b 

2/2 n 2 fl 

critical point by — =l + g— ; it follows therefore that — is little larger in the second 



n ~ g n 

case than in the first. Therefore this critical point has not been already reached where 
the integration stops, but will occur shortly afterwards. 

It is obvious that the amplitude of the nutation cannot increase for an indefinite 
time, because the critical relationship is only exactly satisfied for a single instant. 
In fact, the problem is one of far greater complexity than that of ordinary disturbed 
rotation. The system is disturbed periodically, but the periodic time of the disturb- 
ance slowly increases, passing through a phase of equality to the free periodic time ; 
the problem is to find the amplitude of the oscillations when they are at their maximum, 
and to find the mean configuration of the system some time before and some time 
after the maximum, when the oscillations are small. This problem does not seem to 
be soluble, unless we take into account the slow variation of the argument in the 
periodic disturbing term; and when the argument varies, the disturbing term is not 
strictly a simple time harmonic. 

In the case of the viscous spheroid, the question would be further complicated by 
the fact that when the nutation becomes large, a new series of bodily tides is set up 
by the effects of inertia. 

I have been unable to make a satisfactory examination of this problem, but as far as 
I have gone it appeared to me probable that the mean obliquity of the axis of the 
spheroid would not be affected by the passage of the system through a phase of large 
nutation ; and although I cannot pretend to say how large the nutation might be, yet 
I consider it probable that the amplitude would not have time to increase to a very 
wide extent.* 

* I believe that I shall be able to show in an investigation, as yet incomplete, that when this critical 
phase is reached, the plane of the lunar orbit is nearly coincident with the equator of the earfch. As 
the amplitude of this nutation depends on the sine of the obliquity of the equator to the lunar orbit, it 
seems probable that the nutation would not become considerable. — June 30, 1879, 
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Throughout all the preceding investigations, the periodic inequalities have been 
neglected. Now a full development of the couples H, jflft, iHt which are due to the 
tides, shows that there occur terms of speeds n— 2/2, and n — 4/2 in the first two, and 
of speeds 2/2 and 4/2 in the last. The terms in n — 2/2 in % and 'ffil will clearly 
give rise to an increasing nutation at the critical point which we are considering, 
but they will be so very much smaller than those arising out of the attraction on 
the permanent equatorial protuberance that they may be neglected. The terms in 
n— 4/2 are multiplied by very small quantities, and I think it may safely be assumed 

that the system would pass through the critical phase where \ — |-w=4/2 with sufficient 

rapidity to prevent the nutation becoming large. 

If we were to go to higher orders of approximation in the disturbing forces, it is 
clear that we should meet with an infinite number of critical phases, but the coefficients 
representing the amplitudes of the resulting nutations would be multiplied by such 
small quantities that they may safely be neglected. 



§18. The initial condition of the earth and moon. 






It is now supposed that, when the earth's rotation has been tracked back to where 
it is equal to twice the moon s orbital motion, the obliquity to the plane of the lunar 
orbit has become zero. Then it is clear that, as long as there is any relative motion 
of the earth and moon, the tidal friction and reaction must continue to exist, and n 
and /2 must tend to an equality, The previous investigation shows also that for small 
viscosity, however nearly n approaches /2, the position of zero obliquity is dynamically 
stable. 

As n is approaching /2, the changes must have taken place more and more slowly in 
time. For if the earth was a cooling spheroid, it is unreasonable to suppose that the 
process of becoming less stiff in consistency (which has hitherto been supposed to be 
taking place, as we go backwards in time) could ever have been reversed ; and if it 
were not reversed, then the lunar tides must have lagged by less and less, as more and 
more time was given by the slow relative motion of the two bodies for the moon's 
attraction to have its full effect. Hence the effects of the sun's attraction must again 
become sensible, after passing through a phase of insensibility — a phase perhaps short 
in time, but fertile in changes in the system. I shall not here make the attempt to 
trace the reappearance of these solar terms. 

It is, however, possible to make a rough investigation of what must have been the 
initial state from which the earth and moon started the course of development, which 
has been tracked back thus far. To do this, it is only necessary to consider the equa- 
tion of conservation of moment of momentum. 

* For further consideration of this subject, see a paper on the " Secular Effects of Tidal Friction," 
i Proc. Roy. Soc.,' No. 197, 1879. The arithmetic of this section has been recomputed since the paper 
was presented. 
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n - ■ fi /AT" 



When the obliquity is neglected, that equation may be written — = 1 +/x^ l-i n 

and it is proposed to find what values of n would make n equal to fl. 

In the course of the above investigation four different starting points were taken, 
viz.: those at the beginning of each period of integration. There are objections to 
taking any one of these, to give the numerical values required for the solution of the 
above equation ; for, on the one hand, the errors of each period accumulate on the 
next, and therefore it is advantageous to take one of the early periods ; whilst, on the 
other hand, in the early periods the values of the quantities are affected by the sensi- 
bility of the solar terms, and by the obliquity of the ecliptic. The beginning of the 
fourth period was chosen, because by that time the solar terms had become insigni- 
ficant. At that epoch I found log ?i — 3*84753, when the present tropical year is the 
unit of time, and ^ = '6659, /x being the ratio of the orbital moment of momentum 
to the earth's moment of momentum; also log 5= 5*39378 — 10, s being a constant. 
Now put x 3 z=n=f2, and we have 

x 4> — (1 -\-fjL)n x + -=0 

Then substituting the numerical values, 

as* — 11727x+40385~ 

This equation has two real roots, one of which is nearly equal to y / 11727, and the 
other to 40385 — 11727. By Horner's method these roots are found to be 21*4320 
and 3'4559 respectively. These are the two values of the cube root of the earth's 
rotation, for which the earth and moon move round as a rigid body. 

The first gives a day of 5 hours 36 minutes, and the second a day of about 
55^ m. s. days. 

The latter is the state to which the earth and moon tend, under the influence of 
tidal friction (whether of oceanic or bodily tides) in the far distant future. For this case 
Thomson and Tait give a day of 48 of our present days;* the discrepancy between 
my value and theirs is explicable by the fact that they are considering a heterogeneous 
earth, whilst I treat a homogeneous one. Since on the hypothesis of heterogeneity the 
earth's moment of inertia is about ^-Ma 3 , whilst on that of homogeneity it is |- Ma 3 , and 
since the f- which occurs in the quantity s enters by means of the expression for the 
earth's moment of inertia, it fallows that in my solution /x has been taken too small in 
the proportion 5 : 6. Hence if we wish to consider che case of heterogeneity, we must 
solve the equation x^— 12664^+48462 = 0. r pj :ie ^ w0 r oots of this equation are such 
that they give as the corresponding lengths of the day, 5 hours 16 minutes and 40*4 days 
respectively. The remaining discrepancy (between 40 and 48) is doubtless due in part 

* 'Nat. Phil.,' § 276. They say: — " It is probable that the moon, in ancient times liquid or viscous in 
its outer layer or throughout, was thus brought to turn always the same face to the earth/' In the new 
edition (1879) the ultimate effects of tidal friction are considered, 
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to the crude method of amending the solution, but also to the fact that they partly in- 
clude the obliquity in one way, whilst I partly include it in another way, and I include 
a large part of the solar tidal friction whilst they neglect it. It is interesting to note 
that the larger root, which gives the shorter length of day, is but little affected by the 
consideration of the earth's heterogeneity. 

With respect to the second solution (56 days), it must be remarked that the sun's 
tidal friction will go on lengthening the day even beyond this point, but then the 
lunar tides will again come into existence, and the lunar tidal friction will tend in part 
to counteract the solar. The tidal reaction will also be reversed, so that the moon 
will again approach the earth. Thus the effect of the sun is to make this a state of 
dynamical instability. 

The first solution, where both the day and month are 5 hours 36 minutes long, is 
the one which is of interest in the present inquiry, for this is the initial state towards 
which the integration has been running back. 

"* This state of things is one of dynamical instability, as may be shown as follows : — 

First consider the case where the sun does not exist. Suppose the earth to be 
rotating in about 5^ hours, and the moon moving orbitally around it in a little less 
than that time. Then the motion of the moon relatively to the earth is consentaneous 
with the earth's rotation, and therefore the tidal friction, small though it be, tends to 
accelerate the earth's rotation ; the tidal reaction is such as to tend to retard the moon's 
linear velocity, and therefore increase her orbital angular velocity, and reduce her 
distance from the earth. The end will be that the moon falls into the earth. 

This subject is graphically illustrated in a paper on the " Secular Effects of Tidal 
Friction," read before the Royal Society on June 19, 1879. 

Secondly, take the case where the sun also exists, and suppose the system started 
in the same way as before. Now the motion of the earth relatively to the sun is rapid, 
and such that the solar tidal friction retards the earth's rotation ; whilst the lunar 
tidal friction is, as before, such as to accelerate the rotation. 

Hence if the viscosity be very large the earth's rotation may be accelerated, but if it 
be not very large it will be retarded. The tidal reaction, which depends on the lunar 
tides alone, continues negative, and the moon approaches the earth as before. Thus 
after a short time the motion of the moon relatively to the earth is more rapid than 
in the previous case, whatever be the ratio between solar and lunar tidal friction. 
Hence in this case the moon will fall into the earth more rapidly than if the sun did 
not exist, and the dynamical instability is more marked. 

If, however, the day were shorter than the month, the moon must continually recede 
from the earth, until it reaches the outer limit of a day of 56 m. s. days. 

There is one circumstance which might perhaps decide that this should be the 
direction in which the equilibrium would break down ; for the earth was a cooling 

* From here to the end of the section a good many alterations have been made since the paper was 
presented. — July 5, 1879. 
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body, and therefore probably a contracting one, and therefore its rotation would tend 
to increase. Of course this increase of rotation is partly counteracted by the solar 
tidal friction, but on the present theory, the mere existence of the moon seems to show 
that it was not more than counteracted, for if it had been so the moon must have been 
drawn into and confounded with the earth. 

This month of 5 hours 36 minutes corresponds to a lunar distance of 2*52 earth's 
mean radii, or about 10,000 miles; the month of 5 hours 16 minutes corresponds to 
2*39 earth's mean radii ; so that in the case of the earth's homogeneity only 
6,000 miles intervene between the moon's centre and the earth's surface, and even 
this distance would be reduced if we treated the earth as heterogeneous. This small 
distance seems to me to point to a break-up of the earth-moon mass into two bodies at 
a time when they were rotating in about 5 hours ; for of course the precise figures 
given above cannot claim any great exactitude (see also Section 23). 

It is a material circumstance in the conditions of the breaking-up of the earth into 
two bodies to consider what would have been the ellipticity of the earth's figure when 
rotating in 5^ hours. Now the reciprocal of the ellipticity of a homogeneous fluid or 
viscous spheroid varies as the squaxe of the period of rotation of the spheroid. The 
reciprocal of the ellipticity for a rotation in 24 hours is 232, and therefore the reciprocal 
of the ellipticity for a rotation in 5-| hours is (||) 3 of 232= 2 X 3 - - 4 -X 232=12*2. 

Hence the ellipticity of the earth when rotating in 5^ hours is Yg-th. 

The conditions of stability of a rotating mass of fluid are as yet unknown, but 
when we look at the planets Jupiter and Saturn, it is not easy to believe that an 
ellipticity of x^-th is sufficiently great to cause the break-up of the spheroid. 

A homogeneous fluid spheroid of the same density as the earth has its greatest 
ellipticity compatible with equilibrium when rotating in 2 hours 24 minutes.* 

The maximum ellipticity of all fluid spheroids of the same density is the same, and 
their periods of rotation multiplied by the square root of their densities is a function of 
the ellipticity only. Hence a spheroid, which rotates in 4 hours 48 minutes, will be in 
limiting equilibrium if its density is (~) 3 or \ of that of the earth. If this latter 
spheroid had the same mass as the earth, its radius would be lf\ or 1*59 of that of 
the earth. If therefore the earth had a radius of 6,360 miles, and rotated in 4 hours 
48 minutes, it would just have the maximum ellipticity compatible with equilibrium. 
It is, however, by no means certain that instability would not have set in long before 
this limiting ellipticity was reached. 

In Part III. I shall refer to another possible cause of instability, which may perhaps 
be the cause of the break-up of the earth into two bodies. 

It is easy to find the minimum time in which the system can have passed from this 
initial configuration, where the day and month are both h\ hours, down to the present 

* Pratt's ; Fig', of Earth,' 2nd edition., Arts. 68 and 70. 
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condition. If we neglect the obliquity of the ecliptic, the equation (57) of tidal 
reaction, when adapted to the case of a viscous spheroid, becomes 

u- =f — sin 4€-, 
r At *%n Q l 

Now it is clear that the rate of tidal reaction can never be greater than when 
sin 4€ x =l, when the lunar semi-diurnal tide lags by 22|°. Then since t=-~, we shall 
obtain the minimum time by integrating the equation 



Whence 

Now £= I ~~ ) , and we have found by the solution of the biquadratic that the initial con- 
dition is given by /2*=21 , 4320 ; also with the present value of the month /2 ' l =4 , 38, 
the present year being in both cases the unit of time. Hence it follows that £ is very 

Oil fr/yt 

nearly '2, and £ 13 may be neglected compared with unity. Thus—-£= — — j. 
Now /x= 4*007 and ^° is 86,844,000 years. 

T o 

Hence —^=53,540,000 years. 

Thus we see that tidal reaction is competent to reduce the system from the initial 
state to the present state in something over 54 million years. 

The rest of the paper is occupied with the consideration of a number of miscellaneous 
points, which it was not convenient to discuss earlier. 

§ 19. The change in the length of year. 

The effects of tidal reaction on the earth's orbit round the sun have been neglected ; 
I shall now justify that neglect, and show by how much the length of the year may 
have been altered. 

It is easy to show that the moment of momentum of the orbital motion of the moon 

O 
and earth round their common centre of inertia is -777, where C is the earth's moment of 



av x 3 



, , (1 + ") 



inertia, and s=- 

The moment of momentum of the earth's rotation is obviously Cn. The normal 
to the lunar orbit is inclined to the earth's axis at an angle i. Hence the resultant 
moment of momentum of the moon and earth is 

C\ ft S +7~7^+~7T l COsi 
MDCCCLXXIX. 3 U 
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The change in this quantity from one epoch to another is the amount of moment 
of momentum of the moon-earth system which has been destroyed by solar tidal 
friction. This destroyed moment of momentum reappears in the form of moment of 
momentum of the moon and earth in their orbital motion round the sun. 

Now at the beginning of the integration of Section 17, that is to say at the present 
time, I find that when the present year is taken as the unit of time, the resultant 
moment of momentum of the moon and earth is 11369 C. 

At the end of the third period of integration (after which the solar terms were 
neglected), and when the obliquity has become 15° 22', I find the same quantity to be 

1 1 U iio \j* 

Hence the loss of moment of momentum is 256 C, or 102*4 Md z * 

Now at the present time the moment of momentum of the moon and earth in their 

orbit is (A£-]-m)f2c^=Ma % . (-) 12 / ; 7- is clearly the sun's parallax, and with the 

present unit of time 12, is 2tt\ 

Hence the loss of moment of momentum is equal to the present moment of 

momentum of orbital motion multiplied by — — (sun's parallax) 2 . 

But the moment of momentum of the earth's and moon's orbital motion round the 

sun varies as f2~*; hence the loss of moment of momentum corresponding to a change 

Sf2 
of f2 / to 12,-^-812, is the present moment of momentum multiplied bj^-~ 9 whence it 

is clear that 

Sf2, 102*4 v / ; n \o 

7r= 3 ^TTr/ x (sun s P arallax ) 2 - 



812 

But the shortening of the year is -—^ of a year ; taking therefore the sun's parallax 

all , 

as 8"*8, we find that at the end of the third period of integration the year was shorter 
than at present by 

102-4 82 / 8-87T Y 3 o/. t ^r rsr> A ^ 

X ~^ X 83 X \648~OOo) X 365 " 25 X 86 > 400 seconds, 

which will be found equal to 2*77 seconds. 

Thus the solar tidal reaction had only the effect of lengthening the year by 2f seconds, 
since the epoch specified as the end of the third period of integration. The whole 
change in the length of year since the initial condition to which we traced back the 
moon would probably be very small indeed, but it is impossible to make this assertion 
positively, because, as observed above, the solar effects must have again become sensible, 
after passing through a period of insensibility. 
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§ 20. Terms of the second order in the tide-generating potential. 

The whole of the previous investigation has been conducted on the hypothesis that 
the tide-generating potential, estimated per unit volume of the earth's mass, is 
W7T 2 (cos 2 PM— •§■),* but in fact this expression is only the first term of an infinite 
series. I shall now show what kind of quantities have been neglected by this treat- 
ment. According to the ordinary theory, the next term of the tide-generating 
potential is 

V 2 =ii^Y(f cos 3 PM-f cos PM) 

Although for my own satisfaction I have completely developed the influence of this 
term in a similar way to that exhibited at the beginning of this paper, yet it does not 
seem worth while to give so long a piece of algebra ; and I shall here confine myself to 
the consideration of the terms which will arise in the tidal friction from this term, in 
the potential, when the obliquity is neglected. A comparison of the result with the 
value of the tidal friction, as already obtained, will afford the requisite information as 
to what has been neglected. 

Now when the obliquity is put zero (see Plate 36, fig. 1), 

cos PM=sin 6 sin(<£— co) 

where co is written for n—tl for brevity. Then 

cos 3 PM= J sin 3 8 sin (<£ — co) — \ sin 3 6 sin 3 (<£ — co) 
and 

cos 3 PM-f cos PM=/ - sin 0(1 - 5 cos 3 9) sin (<£-a>) -J sin 3 6 sin 3(^-a>). 



Then since 



m/r\ 3 5 r 3 5 

W 7[c) 2= WT 7 3 



therefore 

Y 2 -7-iv T -r s = —1-2" sin 3 ^ sin 3(<£— co) +i sin 6(1 — 5 cos 2 6) sin (cft—co) 

G 

If sin 3(</>-—a)) and sin (<[)—- co) be expanded, we have V s in the desired form, viz. : a 
series of solid harmonics of the third degree, each multiplied by a simple time har- 
monic. Now if w 3 S 3 cos (vt-\-rj) be a tide-generating potential, estimated per unit 
volume of a homogeneous perfectly fluid spheroid of density w, S 3 being a surface har- 
monic of the third order, then the equilibrium tide due to this potential is given by 

7 a- 3 cr 7a 

or=-r-S 3 cos (vt-\-rj), or ~= : 1 r^~ $ 3 cos(i>£+??). Hence just as in Section 2, the tide- 

* See Section 1 . 
3 u 2 



514 MR, Q. H. DARWIN COT THE PRECESSION OF A VISCOUS SPHEROID, 

generating potential of the third order due to the moon will raise tides in the earth, 
when there is a frictional resistance to the internal motion, given by 



<r 7 ra 

a lOgc 



■iV^sin 3 & sin 3 (<£— <o+f) +W S] * n 0(1 — 5 cos 2 6) sin (<f)—a)+f) 



Now or is a surface harmonic of the third order, and therefore the potential of this 
layer of matter, at an external point whose coordinates are r, 0, <£, is 

4 /a\ 4 3 Ma? 

Hence the moment about the earth's axis of the forces which the attraction of the 

3 Jyfh'ii ft ci u 

distorted spheroid exercises on a particle of mass m, situated at r, 9, <j>, is - — — — , 

» IVl *Y)'l fl* / rr R f»- 

Now if this mass be equal to that of the moon, and r=c, then z— ^— =--ikra 2 =--C, 

IT I C ( C 

where, as before, C is the moment of inertia of the earth. 

Hence the couple jpt 3 , which the moon's attraction exercises on the earth, is given 

K ry ft ff '77* f7T 

by j$&2— — f7~OT7i where after differentiation we put Q-=- - and ^>=-+<y. 

Now 

— 77— to - - * [f .Fsin 3 cos 3(<£-~ o>+jQ —\F' sin 0(1 — 5 cos 2 #) cos (<£—«+/')] 

TT 
JLXo.ll.Ofc5 



c ' s W — \ 2 / V 2 

-sl^sin Sf+i-F" sin/' 



In the case of viscosity 



Therefore 



i^=cos 3f, F'=cosf 



3 ~ / T (A sin 6/+ ^ sin -2/' 



.W 8\ 
Now if the obliquity had been neglected, the tidal friction J^ l5 due to the term of 



T 3 



the first order in the tide-generating potential, would be given by — p = — A sin4e 1 . 

L gj 

Hence 

?__ i AAV^ s i n Q/+ s i n 2/ /N 



*| \ / \ bill 4tfc-i 



That is to say, this is the ratio of the terms neglected previously to those included. 
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Now according to the theory of viscous tides, 






tan 3/==— 7T— ^--fv=ff(3ft))i 



19u 



#W£ v 7 \2gvja, 

where v is the coefficient of viscosity. 

But throughout the previous work we have written p=——. 

X. J V 
o 

Hence tan 3/=|-§-— , and similarly tan/ =ft~« 

Also tan 2e,= — . 

I will now consider two cases :— 

1st. Suppose the viscosity to be small, then f,f, €-, are all small, and 

sin 6/ = tan 3/ = 2 2 8 sin 2f _ tan / ' __ 2 2 

8^46! tan2e x 19 2 ' sn^! tan 2^ 19 2 
Therefore 

2nd. Suppose the viscosity very great, then 3f 9 f, 2e l are very nearly equal to 

~, and tan (|— s /)= 22^ tan (f — /'Wit^ tan yj— 2eJ=-^, so that we have approxi- 
mately 



and similarly 



So that 



sin 6/ __ sin (tt--6/) _^ 19y2 
sin4e 1 sin(7r— 4e x ) 



22X3 



sin 2 /' - 1 9 
sin 4e x 



— JLil v 9 



a \ 1 w 19/10 1 o\ JL9./^ 

J 8 A2"2V 3 "T"^; — S3l c 



Hence it follows that the terms of the second order may bear a ratio to those of the 
first order lying between xf (-] , or 1*16 ( -j , and -g-f (-) . or '576 l- 



Now at the end of the fourth period of integration in the solution of Section 15, - or 

Co 

the moon's distance in earth s mean radii was 9 ; hence the terms of the second 
order in the equation of tidal friction must at that epoch lie in magnitude between 
- 7 - -th and -fixst of those of the first order. It follows, therefore, that even at that 
stage, when the moon is comparatively near the earth, the effect of the tides of the 
second order (i.e., of the third degree of harmonics) is insignificant, and the neglect of 
them is justified. 

In the case of those terms of this order, which affect the obliquity, a very similar 
relationship to the terms of the lower order would be found to hold good. 

* " Bodily Tides," &c, Phil, Trans., 1879, Part I., Section 5, 
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5 21. On certain other small terms. 

It will be well to advert to certain terms, the neglect of which might be suspected 
of vitiating my results. 

According to the hypothesis of the plastic nature of the earth's mass, that body 
must have been a figure of equilibrium at every time throughout the series of 
changes which are to be followed out. In consequence of tidal friction the earth's 
rotation is diminishing, and therefore its ellipticity (which by the ordinary theory is 

7vCt\ 

■£ — J is also diminishing ; this change of figure might be supposed to exercise a 

material influence on the results, but I will now show that in one respect at least its 

effects are unimportant. 

C—A 
In a previous paper * I showed that, neglecting — — compared with unity, when 

the earth's figure changed symmetrically with respect to the axis of rotation, 

di t + t, . . . d fnA . x 

•-= — — — - sin i cos ^T;(vJ — A) 

dt Gn 2 dtr J 

Now if e be the ellipticity of figure 

C-A=fMa»e 

So that 

Gd~r '~dt~2Ydt~~~~% C~ 

and therefore 

di t+t, . . . jg," 

~= sm % cos % — 

dt pi 

Now numerical calculation shows that at present -==——, and since 7^- sin i cos i is 

1 g 10 7 Gn 

of the same order of magnitude as — > — (on which the changes of obliquity have been 

shown to depend), it follows that this term is fairly negligeable compared with those 
already included in the equations. As far as it goes, however, this term tends in the 
direction of increasing the obliquity with the time.t 

* " On the Influence of Geological Changes," &c, Phil. Trans, Yol. 167, Part I., page 272, Section 8. 
The notation is changed, and the equation presented in a form suitable for the present purpose. 

t In a paper in the ' Phil. Mag.,' March, 1877, I suggested that the obliquity might possibly be due to 
the contraction of the terrestrial nebula in cooling ; I there neglected tidal friction and assumed the con- 
servation of moment of momentum to hold good for the earth by itself, so that the ellipticity was con- 
tinually increasing with the time. I did not at that time perceive that this increase of ellipticity was 
antagonistic to the effects of contraction. Though the work of that paper is correct, as I believe, yet the 
fundamental assumption is incorrect, and therefore the results are not worthy of attention, 
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[It will however appear, I believe, that this secular change of ellipticity of the 
earth's figure will exercise an important influence on the plane of the lunar orbit and 
thereby will affect the secular change in the obliquity of the ecliptic. The investiga- 
tion of this point is however as yet incomplete.]"* 

The other small term which I shall consider arises out of the ordinary precession, 
together with the fact that the tide-generating force diminishes with the time on 
account of the tidal reaction on the moon. 

The differential equations which give the ordinary precession are in effect (compare 
equations (26)) 

— i= t—^ — sin % cos % sm n 

at C 

CvGOq \j "~A ... 

-77= — t— — sm % cos % cos n 

dt C 

and they give rise to no change of obliquity if r be constant, but 

T=5=r.fl-«(fw 



when t is small. 

Also - =^=—— =i~. Hence as far as regards the change of obliquity the 

C 4# g 

equations may be written 

dco-, 3r n % 2 /^|?\ . . , . 

dt g \dt ) 



dt g \dt 



sm * cos ^ t cos ^ 



Then if we regard all the quantities, except t> on the ridxfc-hand sides of these 
equations as constants and integrate, we have 

a)i = — ( -^ sm % cos &{ n£ cos n— sm w * 

1 8W 






£o 3 = — s i ~ i sin i cos i{nt sin ?i+ cos w} 



And if these be substituted in the geometrical equations (1) we have 

cli 3r n . . JdF 
-7 = — sm t cos & -if 
a£ g \cw 

* Added July 3, 1879. 
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Now by comparing this with the small term due to the secular change of figure of 
the earth, we see that it is fairly negligeable, being of the same order of magnitude as 
that term. As far as it goes, however, it tends to increase the obliquity of the ecliptic. 



§ 22. The change of obliquity and tidal friction due to an annular satellite. 

Conceive the ring to be rotating round the planet with an angular velocity il, let its 

radius be c, and its mass per unit length of its arc——, so that its mass is m. Let cl be 

the length of the arc measured from some point fixed in the ring up to the element 
c8l ; and let fit be the longitude of the fixed point in the ring at the time t. Let SV 

be the tide-generating potential due to the element —3^. Then we have by (5) 

Where the suffixes to the functions indicate that ft-\-l is to be written for fl. Then 
integrating all round the ring from 1=0 to 1=2it it is clear that 

V 

: —p 2 q 2 sin 2 6 cos 2(<£— n)-\-2pq(]f-™(f) sin 6 cos 6 cos {$—n) 



+ (i- cos 2 0)1(1 -6pV) 

which is the tide-generating potential of the ring. 

Hence, as in Section 2, the form of the tidally-distorted spheroid is given by (9), 
save that E ± , E 2) JE\ } E' 9J E" are all zero. Also, as in that section, the moments of 
the forces which the tidally- distorted spheroid exerts on the element of ring are 

/ QTb Xmi-CCi d(T u(T\ 

f(— 81) ~ s \^~Tf-- t~T~ )> &c., &c, where £r, rjr P £r are put equal to the rectangular 

coordinates of the element of ring, whose annular coordinate is I. 

Now if x, y, z are the direction cosines of the element, equations (7) are simply 
modified by f2 being written f2-\-l. Hence the couples due to one element of ring 
may be found just as the whole couples were found before, and the integrals of the 
elementary couples from Z=0 to 2n are the desired couples due to the whole ring. 
Now a little consideration shows that the results of this integration may be written 
down at once by putting JE l9 E%, E\, E \, E ff zero in (15), (16), and (21). Thus in 
order to determine the change of obliquity and the tidal friction due to an annular 
satellite, we have simply the expressions (33) and (34), save that rr / must be replaced 

by** 

It thus appears that an annular satellite causes tidal friction, in its planet, and that 
the obliquity of the planet's axis to the ring tends to diminish, but both these 
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effects are evanescent with the obliquity. Since this ring only raises the tides which 
are called sidereal semi-diurnal and sidereal diurnal, and since we see by (57), 
Section 1 4, that tidal reaction is independent of those tides, it follows that there is 
no tangential force on the ring tending to accelerate its linear motion. If, however, 
the arc of the ring be not of uniform density, there is a slight tendency for the lighter 
parts to gain on the heavier, and the heavier parts become more remote from the 
planet than the lighter. 

§ 23. Double tidal reaction. 

Throughout the whole of this investigation the moon has been supposed to be 
merely an attractive particle, but there can be no doubt but that, if the earth was 
plastic, the moon was so also. To take a simple case, I shall now suppose that both 
the earth and moon are homogeneous viscous spheres revolving round their common 
centre of inertia, and that the moon is rotating on her own axis with an angular 
velocity a), and that their axes are parallel and perpendicular to the plane of their 
orbit. Then the whole of the argument with respect to the earth as disturbed by 
the moon, may be transferred to the case of the moon as disturbed by the earth. 

All symbols which apply to the moon will be distinguished from those which apply 
to the earth by an accent. 

Then from (21) or (43) we have 



/o 



p/ 2 ft' > DLXX ^^ 1 



and the equation which gives the lunar tidal friction is 



j% 



Now 



at ~ g v ' 



T" ii. -— — jjfm,— — — 1 ff 

— g, __— \/ i — | 

2 c 3 y/a' s 



and 



So that 



it — *^ -s — ■ —• — • — - — , —— ■» - — -ft 

** a f ba w w* 



t " / wa 4 \ 6 T" 



Also 



g' \vfa'*/ g 

C irfaf* 

C vjcfi 



(90) 



and therefore 



t 3 vj 2 a 



= JL — qJij, Ac 

C 2 a w'*a! l 



g w 4 a 
MDCCCLXXTX. 3 X 
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Now the force on the moon tangential to her orbit, results from a double tidal 
reaction. By the method employed in Section 14, the tangential force due to the 
earth's tides is 

T$ri Or". 
= — = - sin 461 

r 2r g L 

and similarly the tangential force due to the moon's tides is 



r 2r g w f *a! 



T / zzz — ^-= — — sin 4e\ 



and the whole tangential force is (T+T"). 

Hence following the argument of that section, the equation of tidal reaction becomes 



2 



d£ , t 



?#% 



sin 4 ex +-777 sin ie\ 






Then taking the moon's apparent radius as 16', and the ratio of the earth's mass to 
that of the moon as 82, we have - = 3*567 and — =1*806 (so that taking %v as 5^, the 



%VOb 

specific gravity of the moon is 3), and hence ~^- ;= 11*64. 

At first sight it would appear from this that the effect of the tides in the moon was 
nearly twelve times as important as the effect of those in the earth, as far as concerns 
the influence on the moon's orbit, and hence it would seem that a grave oversight has 
been made in treating the moon as a simple attractive particle ; a little consideration 
will show, however, that this is by no means the case. 

Suppose that v 9 v are the coefficients of viscosity of the moon and earth respec- 
tively ; then the only tides which exist in each body being those of which the speeds 
are 2(a)— fl), 2(n—S2) in the moon and earth respectively, 



But 



and hence 



J rt , 19i/(©-/2) , , rt 19v(n-I2) 

tan 2e , = -—7-7-; — and tan 2e^ =- — 

1 g aw ga/w 



q aiv =gaw[ — 



tan 2e , = - —77 tan 2ci 

1 n — il v\wa J 



ivct> \ * 

It will be found that (-7-,) =41*10. It is also almost certain that v must for a 

1 w a / 
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long time be greater than v, because the moon being a smaller body must have 
stiffened quicker than the earth. Hence unless co — 12 is very much less than n — 12, 
e\ must be larger than c 2 . Therefore if in the early stages of development the earth 
had a small viscosity, it is probable that the effects of the moon s tides on her own 
orbit must have had a much more important influence than had the tides in the earth. 
I shall now show, however, that this state of things must probably have had so short 
a duration as not to seriously affect the investigation of this paper. By (89) and (90) 
we have, as the equation which determines the rate of tidal friction reducing the 
moon's rotation round her axis, 

dco , tV wo? \ 3 . . 

— =— • h~~ \-TTs>) sui 4c , 
dt * g \w f a!*j J 

—7-7- ) = 12,148 ; and hence, for the same values of e\ and e,, the moon's rotation 

round her axis is reduced 12,000 times as rapidly as that of the earth round its axis, 
and therefore in a very short period the moon's rotation round her axis must have 
been reduced to a sensible identity with the orbital motion. As a) becomes very 
nearly equal to /2, sin Ae\ becomes very small. Hence the term in the equation of tidal 
reaction dependent on the moon's own tides must have become rapidly evanescent. 
Now while this shows that the main body of our investigation is unaffected by the 
lunar tide, there is one slight modification of them to which it leads. 

In Section 18 we traced back the moon to the initial condition, when her centre 
was 10,000 miles from the earth's centre. If lunar tidal friction had been included, 
this distance would have been increased ; for the coefficient of x in the biquadratic 

on n. ° QJO CL 

(viz. : 11,727) would have to be diminished by — -(a>-— g> ). Now — - is very nearly 

ToVotK and the unit of time being the year, it follows that we should have to suppose 
an enormously rapid primitive rotation of the moon round her axis, to make any 
sensible difference in the configuration of the two bodies when her centre of inertia 
moved as though rigidly connected with the earth's surface. 

The supposition of two viscous globes moving orbitally round their common centre 
of inertia, and one having a congruent and the other an incongruent axial rotation, 
would lead to some very curious results. 



24. Secular contraction of the earth. 



If the earth be contracting as it cools, it follows, from the principle of conservation 
of moment of momentum, that the angular velocity of rotation is being increased. Sir 
William Thomson has, however, shown that the contraction (which probably now 
only takes place in the superficial strata) cannot be sufficiently rapid to perceptibly 
counteract the influence of tidal friction at the present time, 

* Rewritten in July, 1879. 
o X £, 
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The enormous height of the lunar mountains compared to those in the earth seems, 
however, to give some indications that a cooling celestial orb must contract by a 
perceptible fraction of its radius after it has consolidated/'" Perhaps some of the 
contraction might be due to chemical combinations in the interior, when the heat had 
departed, so that the contraction might be deep-seated as well as superficial. 

It will be well, therefore, to point out how this contraction will influence the initial 
condition to which we have traced back the earth and moon, when they were found 
rotating as parts of a rigid body in a little more than 5 hours. 

Let C, C be the moment of inertia of the earth at any time, and initially. Then 
the equation of conservation of moment of momentum becomes 



C % \ \ fl 



And, the biquadratic of Section 18 which gives the initial configuration becomes 



C ^Cs 



a A_( 1+ rf^ x+ b± =0 



The required root of this equation is very nearly equal to 



(i+/^?~ 



Now 



C n 
x s =f2 ; hence fl is nearly equal to (l+/x)-~— . But in Section 18, when C was equal 

to C , it was nearly equal to (l+/x)?7 . Therefore on the present hypothesis, the value 



* Suppose a sphere of radius a to contract until its radius is a + $a, but that, its surface being incom- 
pressible, in doing so it throws up n conical mountains, the radius of whose bases is b, and their height h, 



and let b be large compared with h. The surface of such a cone is 7rb^/Ifi + b i2 ='7r(b 2 +^~h?), Hence the 
excess of the surface of the cone above the area of the base is -jp/^, and 4>7ra 2 =4t7r(a+Say + ^n7rh 2 . 

Therefore =-«|-| . 

* a loyv 

Then suppose we have a second sphere of primitive radius a', which contracts and throws up the same 

number of mountains : then similarlv — — = — i — I and 1 — =( — S . ISTow let these two spheres be 

J a' 16\a/ a' a \ha'J 

the earth and moon. The height of the highest lunar mountain is 23,000 feet (Grant's * Physical Astron.,' 

p. 229), and the height of the highest terrestrial mountain is 29,000 feet; therefore we may take 

h' a 1 lied 

-=|f. Also - = -2729 (Heesohel's 'Astron.,' Section 404). Therefore — =ff of '2729=-344. and 

fo cl ha 

P la '\ oo A' A 2 __ 8 a' $a 

\TT ) =*lloo or I — I =8*45. Hence — -4- — = 8| ; whence it appears that, if both lunar and terrestrial 

mountains are due to the crumpling of the surfaces of those globes in contraction, the moon's radius has 
been diminished by about eight times as large a fraction as the earth's. 

This is, no doubt, a very crude way of looking at the subject, because it entirely omits volcanic action 
from consideration, but it seems to justify the assertion that the moon has contracted much more than 
the earth, since both bodies solidified, 



AND ON THE REMOTE HISTORY OF THE EARTH. 523 


of fl as given in that section must be multiplied by -~ ; and the periodic time must be 

\j 


multiplied by — . But in this initial state C is greater than C ; hence the periodic 

time when the two bodies move round as a rigid body is longer, and the moon is 
more distant from the earth, if the earth has sensibly contracted since this initial 
configuration. 

If, then, the theory here developed of the history of the moon is the true one, as I 
believe it is, it follows that the earth cannot have contracted since this initial state by 
so much as to considerably diminish the effects of tidal friction, and it follows that 
Sir William Thomson's result as to the present unimportance of the contraction must 
have always been true. 

If the moon once formed a part of the earth we should expect to trace the changes 
back until the two bodies were in actual contact. But it is obvious that the data at 
our disposal are not of sufficient accuracy, and the equations to be solved are so com- 
plicated, that it is not to be expected that we should find a closer accordance, than has 
been found, between the results of computation and the result to be expected, if the 
moon was really once a part of the earth. 

It appears to me, therefore, that the present considerations only negative the 
hypothesis of any large contraction of the earth since the moon has existed. 



PART III. 



Summary and discussion of results. 



The general object of the earlier or preparatory part of the paper is sufficiently 
explained in the introductory remarks. 

The earth is treated as a homogeneous spheroid, and in what follows, except where 
otherwise expressly stated, the matter of which it is formed is supposed to be purely 
viscous. The word " earth " is thus an abbreviation of the expression " a homogeneous 
rotating viscous spheroid ; " also wherever numerical values are given they are taken 
from the radius, mean density, mass, &c, of the earth. 

The case is considered first of the action of one tide-raising body, namely, the moon. 
To simplify the problem the moon is supposed to move in a circular orbit in the 
ecliptict — that plane being the average position of the lunar orbit with respect to the 

* This part has been altered in accordance with the several additions and alterations occurring above. 
The results of subsequent investigations have modified the interpretation to be put on several of the results 
here obtained. I have, moreover, had the advantage of discussing several points with Sir William 
Thomson.— July 9, 1879. 

f The effect of neglecting the eccentricity of the moon's orbit is, that we underestimate the efficiency 
of the tidal effects. Those effects vary as the inverse sixth power of r the radius vector, and if T be the 
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earth's axis. The case becomes enormously more complex if we suppose the moon to 
move in an inclined eccentric orbit with revolving nodes. The consideration of the 
secular changes in the inclination of the lunar orbit and of the eccentricity will form 
the subject of another investigation. 

The expression for the moon's tide-generating potential is shown to consist of 13 
simple tide-generating terms, and the physical meaning of this expansion is given in 
the note to Section 8. The physical causes represented by these 13 terms raise 13 simple 
tides in the earth, the heights and retardations of which depend on their speeds and on 
the coefficient of viscosity. 

The 13 simple tides may be more easily represented both physically and analytically 
as seven tides, of which three are approximately semi-diurnal, three approximately 
diurnal, and one has a period equal to a half of the sidereal month, and is therefore 
called the fortnightly tide. 

Then by an approximation which is sufficiently exact for a great part of the investi- 
gation, the semi-diurnal tides may be grouped together, and the diurnal ones also. 
Hence the earth may be regarded as distorted by two complex tides, namely, the semi- 
diurnal and diurnal, and one simple tide, namely, the fortnightly. The absolute heights 
and retardations of these three tides are expressed by six functions of their speeds and 
of the coefficient of viscosity (Sections 1 and 2). 

When the form of the distorted spheroid is thus given, the couples about three axes 
fixed in the earth due to the attraction of the moon on the tidal protuberances are 
found. It must here be remarked that this attraction must in reality cause a tan- 
gential stress between the tidal protuberances and the true surface of the mean 
oblate spheroid. This tangential stress must cause a certain very small tangential 
flow/' 5 " and hence must ensue a very small diminution of the couples. The diminution 
of couple is here neglected, and the tidal spheroid is regarded as being instantaneously 
rigidly connected with the rotating spheroid. The full expression for the couples on 
the earth are long and complex, but since the nutations to which they give rise are 
exceedingly minute, they may be much abridged by the omission of all terms except 
such as can give rise to secular changes in the precession, the obliquity of the ecliptic, 
and the diurnal rotation. The terms retained represent that there are three couples 
independent of the time, the first of which tends to make the earth rotate about an 
axis in the equator which is always 90° from the nodes of the moon's orbit : this 
couple affects the obliquity to the ecliptic ; second, there is a couple about an axis in 

1 1 f T dt 

periodic time of the moon, the average value of ^ is jj — 6 . If c be the mean distance and e the eecen- 

1 l + 3e 3 +^-e 4 
tricity of the orbit, this integral will be found equal to -$ ,-. -yrf . If the eccentricity be small the 

1 1 / 15 \ 1 54 1 

average value of -^ is ^ ( l+^e 2 ); if e is ^g this is ^ of ~ 6 . There are obviously forces tending to 

modify the eccentricity of the moon's orbit. 
* See Part I, of the next paper. 
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the equator which is always coincident with the nodes : this affects the precession ; 
third, there is a couple about the earth's axis of rotation, and this affects the length of 
the day (Sections 3, 4, and 5). All these couples vary as the fourth power of the moon s 
orbital angular velocity, or as the inverse sixth power of her distance. 

These three couples give the alteration in the precession due to the tidal movement, 
the rate of increase of obliquity, and the rate at which the diurnal rotation is being 
diminished, or in other words the tidal friction. The change of obliquity is in reality 
due to tidal friction, but it is convenient to retain the term specially for the change of 
rotation alone. 

It appears that if the bodily tides do not lag, which would be the case if the earth 
were perfectly fluid or perfectly elastic, then there is no alteration in the obliquity, nor 
any tidal friction (Section 7). The alteration in the precession is a very small fraction 
of the precession due to the earth considered as a rigid oblate spheroid. I have some 
doubts as to whether this result is properly applicable to the case of a perfectly fluid 
spheroid. At any rate, Sir William Thomson has stated, in agreement with this 
result, that a perfectly fluid spheroid has a precession scarcely differing from that of a 
perfectly rigid one. Moreover, the criterion which he gives of the negligeability of the 
additional terms in the precession in a closely analogous problem appears to be almost 
identical with that found by me (Section 7). I am not aware that the investigation on 
which his statement is founded has ever been published. The alteration in the pre- 
cession being insignificant, no more reference will be made to it. This concludes the 
analytical investigation as far as concerns the effects on the disturbed spheroid, where 
there is only one disturbing body. 

The sun is now (Section 8) introduced as a second disturbing body. Its independent 
effect on the earth may be determined at once by analogy with the effect of the moon. 
But the sun attracts the tides raised by the moon, and vice versa. Now notwith- 
standing that the periods of the sun and moon about the earth have no common 
multiple, yet the interaction is such as to produce a secular alteration in the position 
of the earth's axis and in the angular velocity of its diurnal rotation. A physical 
explanation of this curious result is given in the note to Section 8. I have dis- 
tinguished this from the separate effect of each disturbing body, as a combined effect. 

The combined effects are represented by two terms in the tide-generating potential, 
one of which goes through its period in 12 sidereal hours, and the other in a sidereal 
day* ; the latter being much more important than the former for moderate obliquities 
to the ecliptic. Both these terms vanish when the earth's axis is perpendicular to the 
plane of the orbit. 

As far as concerns the combined effects, the disturbing bodies may be conceived to be 

* These combined effects depend on the tides which are designated as K x and K 3 in the British Asso- 
ciation's Report on Tides for 1872 and 1876, and which I have called the sidereal semi-diurnal and 
diurnal tides. For a general explanation of this result see the abstract of this paper in the ' Proceedings 
of the Royal Society,' No. 191, 1878. 



526 MR, G. H. DARWIN ON THE PRECESSION OF A VISCOUS SPHEROID 

replaced by two circular rings of matter coincident with their orbits and equal in mass 
to them respectively. The tidal friction due to these rings is insignificant compared 
with that arising separately from the sun and moon. But the diurnal combined effect 
has an important influence in affecting the rate of change of obliquity. The combined 
effects are such as to cause the obliquity of the ecliptic to diminish, whereas the 
separate effects on the whole make it increase — at least in general (see Section 22). 

The relative importance of all the effects may be seen from an inspection of Table III., 
Section 15. 

Section 11 contains a graphical analysis of the physical meaning of the equations, 
giving the rate of change of obliquity for various degrees of viscosity and obliquity. 

Plate 36, figures 2 and 3, refer to the case where the disturbed planet is the earth, 
and the disturbing bodies the sun and moon. 

This analysis gives some remarkable results as to the dynamical stability or 
instability of the system. 

It will be here sufficient to state that, for moderate degrees of viscosity, the position 
of zero obliquity is unstable, but that there is a position of stability at a high obliquity. 
For large viscosities the position of zero obliquity becomes stable, and (except for a 
very close approximation to rigidity) there is an unstable position at a larger obliquity, 
and again a stable one at a still larger one.** 

These positions of dynamical equilibrium do not rigorously deserve the name, since 
they are slowly shifting in consequence of the effects of tidal friction ; they are rather 
positions in which the rate of change of obliquity becomes of a higher order of small 
quantities. 

It appears that the degree of viscosity of the earth which at the present time would 
cause the obliquity of the ecliptic to increase most rapidly is such that the bodily semi- 
diurnal tide would be retarded by about 1 hour and 10 minutes; and the viscosity 
which would cause the obliquity to decrease most rapidly is such that the bodily semi- 
diurnal tide would be retarded by about 2f hours. 

The former of these two viscosities was the one which I chose for subsequent 
numerical application, and for the consideration of secular changes in the system. 

Plate 36, fig. 4 (Section 11), shows a similar analysis of the case where there is only 
one disturbing satellite, which moves orbitally with one-fifth of the velocity of rotation 
of the planet. This case differs from the preceding one in the fact that the position of 
zero obliquity is now unstable for all viscosities, and that there is always one other, 
and only one other position of equilibrium, and that is a stable one. 

This shows that the fact that the earth's obliquity would diminish for large viscosity 
is due to the attraction of the sun on the lunar tides, and of the moon on the solar 

tides. 

It is not shown by these figures, but it is the fact that if the motion of the satellite 

* For a general explanation of some part of these results, see the abstract of this paper in the 
* Proceedings of the .Royal Society,' No. 191, 18 78. 
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relatively to the planet be slow enough (viz. : the month less than twice the day), the 
obliquity will diminish. 

This result, taken in conjunction with results given later with regard to the evolu- 
tion of satellites, shows that the obliquity of a planet perturbed by a single satellite 
must rise from zero to a maximum and then decrease again to zero. If we regard the 
earth as a satellite of the moon, we see that this must have been the case with the moon. 
Plate 36, fig. 5 (Section 12), contains a similar graphical analysis of the various 
values which may be assumed by the tidal friction. As might be expected, the tidal 
friction always tends to stop the planet's rotation, unless indeed the satellite's period 
is less than the planet's day, when the friction is reversed. 

This completes the consideration of the effect on the earth, at any instant, of the 
attraction of the sun and moon on their tides; the next subject is to consider the 
reaction on the disturbing bodies. 

Since the moon is tending to retard the earth's diurnal rotation, it is obvious that 
the earth must exercise a force on the moon tending to accelerate her linear velocity, 
The effect of this force is to cause her to recede from the earth and to decrease her 
orbital angular velocity. Hence tidal reaction causes a secular retardation of the 
moon s mean motion. 

The tidal reaction on the sun is shown to have a comparatively small influence on 
the earth's orbit and is neglected (Sections 14 and 19). 

The influence of tidal reaction on the lunar orbit is determined by finding the dis- 
turbing force on the moon tangential to her orbit, in terms of the couples which have 
been already found as perturbing the earth's rotation ; and hence the tangential force 
is found in terms of the rate of tidal friction and of the rate of change of obliquity. 

It appears that the non-periodic part of the force, on which the secular change in 
the moon's distance depends, involves the lunar tides alone. 

By the consideration of the effects of the perturbing force on the moon's motion, an 
equation is found which gives the rate of increase of the square root of the moon's 
distance, in terms of the heights and retardations of the several lunar tides 
(Section 14). 

Besides the interaction of the two bodies which affects the moon's mean motion, 
there is another part which affects the plane of the lunar orbit ; but this latter effect 
is less important than the former, and in the present paper is neglected, since the moon 
is throughout supposed to remain in the ecliptic. The investigation of the subject will 
however, lead to interesting results, since a complete solution of the problem of the 
obliquity of the ecliptic cannot be attained without a simultaneous tracing of the 
secular changes in the plane of the lunar orbit. 

It appears that the influence of the tides, here called slow semi-diurnal and slow 
diurnal, is to increase the moon's distance from the earth, whilst the influence of the 
fast semi-diurnal, fast diurnal, and fortnightly tide tends to diminish the moon's dis- 
tance ; also the sidereal semi-diurnal and diurnal tides exercise no effects in this 

MDCCCLXXIX. 3 Y 
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respect. The two tides which tend to increase the moon s distance are much larger 
than the others, so that the moon in general tends to recede from the earth. The 
increase of distance is, of course, accompanied by an increase of the moon s periodic 
time, and hence there is in general a true secular retardation of the moon's motion. 
But this change is accompanied by a retardation of the earth's diurnal rotation, and a 
terrestrial observer, taking the earth as his clock, would conceive that the angular 
velocity of an ideal moon, which was undisturbed by tidal reaction, was undergoing a 
secular acceleration. The apparent acceleration of the ideal undisturbed moon must 
considerably exceed the true retardation of the real disturbed moon, and the difference 
between these two will give an apparent acceleration. 

It is thus possible to give an equation connecting the apparent acceleration of the 
moon's motion and the heights and retardations of the several bodily tides in the earth. 

Now there is at the present time an unexplained seeukr acceleration of the moon of 
about 4" per century, and therefore if we attribute the whole of this to the action of 
the bodily tides in the earth, instead of to the action of ocean tides, as was done by 
Adams and Delaunay, we get a numerical relation which must govern the actual 
heights and retardations of the bodily tides in the earth at the present time. 

This equation involves the six constants expressive of the heights and retardations of 
the three bodily tides, and which are determined by the physical constitution of the 
earth. No further advance can therefore be made without some theory of the earth's 
nature. Two theories are considered. 

First, that the earth is purely viscous. The result shows that the earth is either 
nearly fluid — which we know it is not — -or exceedingly nearly rigid. The only traces 
which we should ever be likely to find of such a high degree of viscosity would be in 
the fortnightly ocean tide ; and even here the influence would be scarcely perceptible, 
for its height would be *9 9 2 of its theoretical amount according to the equilibrium theory, 
whilst the time of high water would be only accelerated by six hours and a half. 

It is interesting to note that the indications of a fortnightly ocean tide, as deduced 
from tidal observations, are exceedingly uncertain, as is shown in a preceding paper/'" 
where I have made a comparison of the heights and phases of such small fortnightly tides 
as have hitherto been observed. And now (July, 1879) Sir William Thomson has 
informed me that he thinks it very possible that the effects of the earth's rotation may- 
be such as to prevent our trusting to the equilibrium theory to give even approximately 
the height of the fortnightly tide. He has recently read a paper on this subject 
before the Royal Society of Edinburgh. 

With the degree of viscosity of the earth, which gives the observed amount of secular 
acceleration to the moon, it appears that the moon is subject to such a true secular 
retardation that at the end of a century she is 3 //# l behind the place in her orbit which 
she would have occupied if it were not for the tidal reaction, whilst the earth, considered 
as a clock, is losing 1 3 seconds in the same time. This rate of retardation of the earth 

* See the Appendix to my paper on the " Bodily Tides/' &c,, Phil. Trans., Part L, 1879, 
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is such that an observer taking the earth as his clock would conceive a moon, which 
was undisturbed by tidal reaction, to be 7"'l in advance of her place at the end of a 
century. But the actual moon is 3"'l behind her true place, and thus our observer 
would suppose the moon to be in advance 71 — 3*1 or 4" at the end of the century. 
Lastly, the obliquity of the ecliptic is diminishing at the rate of 1° in 500 million 
years. 

The other hypothesis considered is that the earth is very nearly perfectly elastic. In 
this case the semi-diurnal and diurnal tides do not lag perceptibly, and the whole of the 
reaction is thrown on to the fortnightly tide, and moreover there is no perceptible tidal 
frictional couple about the earth's axis of rotation. From this follows the remarkable 
conclusion that the moon may be undergoing a true secular acceleration of motion of 
something less than 3"*5 per century, whilst the length of day may remain almost un- 
affected. Under these circumstances the obliquity of the ecliptic must be diminishing 
at the rate of 1° in something like 130 million years. 

This supposition leads to such curious results, that I investigated what state of 
things we should arrive at if we look back for a very long period, and I found that 
700 million years ago the obliquity might have been 5° greater than at present, whilst 
the month would only be a little less than a day longer. The suppositions on which 
these results are based are such that they necessarily give results more striking than 
would be physically possible. 

The enormous lapse of time which has to be postulated renders it in the highest 
degree improbable that more than a very small change in this direction has been taking 
place, and moreover the action of the ocean tides has been entirely omitted from 
consideration. 

The results of these two hypotheses show what fundamentally different interpreta- 
tions may be put to the phenomenon of the secular acceleration of the moon. 

Sir William Thomson also has drawn attention to another disturbing cause in the 
fall of meteoric dust on to the earth. * 

Under these circumstances, I cannot think that any estimate having any pretension 
to accuracy can be made as to the present rate of tidal friction. 

Since the obliquity of the ecliptic, the diurnal rotation of the earth, and the moon's 
distance change, the whole system is in a state of flux ; and the next question to be 
considered is to determine the state of things which existed a very long time ago 
(Part II.). This involved the integration of three simultaneous differential equations; 
the mathematical difficulties were, however, so great, that it was found impracticable 
to obtain a general analytical solution. I therefore had to confine myself to a 
numerical solution adapted to the case of the earth, sun, and moon, for one particular 
degree of viscosity of the earth. The particular viscosity was such that, with the 
present values of the day and month, the time of the lunar semi-diurnal tide was 
retarded by 1 hour and 10 minutes ; the greatest possible lagging of this tide is 

* i Glasgow Geological Society/ Vol. III. Address "On Geological Time." 

3 y 2 
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3 hours, and therefore this must be regarded as a very moderate degree of visco- 
sity. It was chosen because initially it makes the rate of change of obliquity a 
maximum, and although it is not that degree of viscosity which will make all the 
changes proceed with the greatest possible rapidity, yet it is sufficiently near that 
value to enable us to estimate very well the smallest time which can possibly have 
elapsed in the history of the earth, if changes of the kind found really have taken 
place. This estimate of time is confirmed by a second method, which will be referred 
to later. 

The changes were tracked backwards in time from the present epoch, and for con- 
venience of diction I shall also reverse the form of speech— -e.g., a true loss of energy 
as the time increases will be spoken of as a gain of energy as we look backwards. 

I shall not enter at all into the mathematical difficulties of the problem, but shall 
proceed at once to comment on the series of tables at the end of Section 15, which 
give the results of the solution. 

The whole process, as traced backwards, exhibits a gain of kinetic energy to the 
system (of which more presently), accompanied by a transference of moment of 
momentum from that of orbital motion of the moon and earth to that of rotation of 
the earth. The last column but one of Table IV. exhibits the fall of the ratio "of the 
two moments of momentum from 4*01 down to *44. The whole moment of momentum 
of the moon-earth system rises slightly, because of solar tidal friction. The change is 
investigated in Section 19. 

Looked at in detail, we see the day, month, and obliquity all diminishing, and the 
changes proceeding at a rapidly increasing rate, so that an amount of change which at 
the beginning required many millions of years, at the end only requires as many thou- 
sands. The reason of this is that the moon s distance diminishes with great rapidity ; 
and as the effects vary as the square of the tide-generating force, they vary as the 
inverse sixth power of the moon's distance, or, in physical language, the height of the 
tides increases with great rapidity, and so also does the moon s attraction. But there 
is a counteracting principle, which to some extent makes the changes proceed slower. 
It is obvious that a disturbing body will not have time to raise such high tides in a 
rapidly rotating spheroid as in one which rotates slowly. As the earth's rotation 
increases, the lagging of the tides increases. The first column of Table I. shows the 
angle by which the crest of the lunar semi-diurnal tide precedes the moon ; we see 
that the angle is almost doubled at the end of the series of changes, as traced back- 
wards. It is not quite so easy to give a physical meaning to the other columns, 
although it might be done. In fact, as the rotation increases, the effect of each tide 
rises to a maximum, and then dies away ; the tides of longer period reach their maxi- 
mum effect much more slowly than the ones of short period. At the point where I 
have found it convenient to stop the solution (see Table IV.), the semi-diurnal effect has 
passed its maximum, the diurnal tide has just come to give its maximum effect, whilst 
the fortnightly tide has not nearly risen to that point. 
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As the lunar effects increase in importance (when we look backwards), the relative 
value of the solar effects decreases rapidly, because the solar tidal reaction leaves the 
earth's orbit sensibly unaffected (see Section 19), and thus the solar effects remain nearly 
constant, whilst the lunar effects have largely increased. The relative value of the 
several tidal effects is exhibited in Tables II. and III. 

Table IV. exhibits the length of day decreasing to a little more than a quarter of its 
present value, whilst the obliquity diminishes through 9°. But the length of the 
month is the element which changes to the most startling extent, for it actually falls 
to iVth of its primitive value. 

It is particularly important to notice that all the changes might have taken place in 
57 million years; and this is far within the time which physicists admit that the earth 
and moon may have existed. It is easy to find a great many verce causce for changes 
in the planetary system ; but it is in general correspondingly hard to show that they 
are competent to produce any marked effects, without exorbitant demands on the 
efficiency of the causes and on lapse of time. 

It is a question of great interest to geologists to determine whether any part of 
these changes could have taken place during geological history. It seems to me that 
this question must be decided by whether or not a globe, such as has been considered, 
could have afforded a solid surface for animal life, and whether it might present a 
superficial appearance such as we know it. These questions must, I think, be answered 
in the affirmative, for the following reasons. 

The coefficient of viscosity of the spheroid with which the previous solution, deals is 

given by the formula —-tan 35° (see Section 11, (40)), when gravitation units of force 

are used. This, when turned into numbers, shows that 2*055 X 10 7 grams weight 
are required to impart unit shear to a cubic centimeter block of the substance in 
24 hours, or 2,055 kilogs. per square centimeter acting tangentially on the upper 
face of a slab one centimeter thick for 24 hours, would displace the upper surface 
through a millimeter relatively to the lower, which is held fixed. In British units 
this becomes, — 13-J tons to the square inch, acting for 24 hours on a slab an inch thick, 
displaces the upper surface relatively to the lower through one-tenth of an inch. It 
is obvious that such a substance as this would be called a solid in ordinary parlance, 
and in the tidal problem this must be regarded as a rather small viscosity. 

It seems to me, then, that w r e have only got to postulate that the upper and cool 
surface of the earth presents such a difference from the interior that it yields with 
extreme slowness, if at all, to the weight of continents and mountains, to admit the 
possibility that the globe on which we live may be like that here treated of. If, 
therefore, astronomical facts should confirm the argument that the w T orld has really 
gone through changes of the kind here investigated, I can see no adequate reason for 
assuming that the whole process was pre-geological. Under these circumstances it 
must be admitted that the obliquity to the ecliptic is now probably slowly decreasing; 
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that a long time ago it was perhaps a degree greater than at present, and that it was 
then nearly stationary for another long time, and that in still earlier times it was 
considerably less/ 7 " 

The violent changes which some geologists seem to require in geologically recent 
times would still, I think, not follow from the theory of the earth's viscosity. 

According to the present hypothesis (and for the moment looking forward in time), 
the moon-earth system is, from a dynamical point of view, continually losing energy 
from the internal tidal friction. One part of this energy turns into potential energy 
of the moon's position relatively to the earth, and the rest developes heat in the 
interior of the earth. Section 16 contains the investigation of the amount which has 
turned to heat between any two epochs. The heat is estimated by the number of 
degrees Fahrenheit, which the lost energy would be sufficient to raise the temperature 
of the whole earth's mass, if it were all applied at once, and if the earth had the specific 
heat of iron. 

The last column of Table IV., Section 15, gives the numerical results, and it appears 
therefrom that, during the 57 million years embraced by the solution, the energy lost 
suffices to heat the whole earth's mass 1760° Fahr. 

It would appear at first sight that this large amount of heat, generated internally, 
must seriously interfere with the accuracy of Sir William Thomson's investigation 
of the secular cooling of the earth ;t but a further consideration of the subject in the 
next paper will show that this cannot be the case. 

There are other consequences of interest to geologists which flow from the present 
hypothesis. As we look at the whole series of changes from the remote past, the 
ellipticity of figure of the earth must have been continually diminishing, and thus the 
polar regions nmst have been ever rising and the equatorial ones falling ; but, as the 
ocean always followed these changes, they might quite well have left no geological 
traces. 

The tides must have been very much more frequent and larger, and accordingly the 
rate of oceanic denudation much accelerated. 

The more rapid alternations of day and night J would probably lead to more sudden 
and violent storms, and the increased rotation of the earth would augment the violence 
of the trade winds, which in their turn would affect oceanic currents. 

Thus there would result an acceleration of geological action. 

The problem, of which the solution has just been discussed, deals with a spheroid of 

* In my paper " On the Effects of Geological Changes on the Earth's Axis," Phil. Trans. 1877, p. 271, 
I arrived at the conclusion that the obliquity had been unchanged throughout geological history. That 
result was obtained on the hypothesis of the earth's rigidity, except as regards geological upheavals. The 
result at which I now arrive affords a warning that every conclusion must always be read along with the 
postulates on which it is based. 

f ' Nat. Phil.,' Appendix. 

J At the point where the solution stops there are just 1,300 of the sidereal days of that time in the 
year, instead of 366 as at present. 
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constant viscosity ; but there is every reason to believe that the earth is a cooling 
body, and has stiffened as it cooled. We therefore have to deal with a spheroid whose 
viscosity diminishes as we look backwards. 

A second solution is accordingly given (Section 17) where the viscosity is variable; no 
definite law of diminution of viscosity is assumed, however, but it is merely supposed 
that the viscosity always remains small from a tidal point of view. This solution gives 
no indication of the time which may have elapsed, and differs chiefly from the preceding 
one in the facfc that the change in the obliquity is rather greater for a given amount of 
change in the moon's distance. 

There is not much to say about it here, because the two solutions follow closely 
parallel lines as far as the place where the former one left off. 

The first solution was not carried further, because as the month approximates in 
length to the day, the three semi-diurnal tides cease to be of nearly equal frequencies, 
and so likewise do the three diurnal tides; hence the assumption on which the solution 
was founded, as to their approximately equal speeds, ceases to be sufficiently accurate. 

In this second solution all the seven tides are throughout distinguished from one 
another. At about the stage where the previous solution stops the solar terms have 
become relatively unimportant, and are dropped out. It appears that (still looking 
backwards in time) the obliquity will only continue to diminish a little more beyond 
the point it had reached when the previous method had become inapplicable. For 
when the month has become equal to twice the day, there is no change of obliquity ; 
and for yet smaller values of the month the change is the other way. 

This shows that for small viscosity of the planet the position of zero obliquity is 
dynamically stable for values of the month which are less than twice the day, while 
for greater values it is unstable; and the same appears to be true for very large vis- 
cosity of the planet (see the foot-note on p. 500). 

If the integration be carried back as far as the critical point of relationship between 
the day and month, it appears that the whole change of obliquity since the beginning 
is y q" . 

The interesting question then arises — Does the hypothesis of the earth's viscosity 
afford a complete explanation of the obliquity of the ecliptic ? It does not seem at 
present possible to give any very conclusive answer to this question ; for the problem 
which has been solved differs in many respects from the true problem of the earth. 

The most important difference from the truth is in the neglect of the secular changes 
of the plane of the lunar orbit; and I now (September, 1879) see reason to believe 
that that neglect will make a material difference in the results given for the obliquity 
at the end of the third and fourth periods of integration in both solutions. It will 
not, therefore, be possible to discuss this point adequately at present ; but it will be 
well to refer to some other points in which our hypothesis must differ from reality. 

I do not see that the heterogeneity of density and viscosity would make any very 
material difference in the solution, because both the change of obliquity and the tidal 
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friction would be affected pari passu, and therefore the change of obliquity for a given 
amount of change in the day would not be much altered. 

Although the effects of the contraction of the earth in cooling would be certainly 
such as to render the changes more rapid in time, yet as the tidal friction would be 
somewhat counteracted, the critical point where the month is equal to twice the day 
would be reached when the moon was further from the earth than in my problem. I 
think, however, that there is reason to believe that the whole amount of contraction 
of the earth, since the moon has existed, has not been large (Section 24). 

There is one thing which might exercise a considerable influence favourable to change 
of obliquity. We are in almost complete ignorance of the behaviour of semi-solids 
under very great pressures, such as must exist in the earth, and there is no reason to 
suppose that the amount of relative displacement is simply proportional to the stress 
and the time of its action. Suppose, then, that the displacement varied as some other 
function of the time, then clearly the relative importance of the several tides might be 
much altered. 

Now, the great obstacle to a large change of obliquity is the diurnal combined 
effect (see Table IV., Section 15) ; and so any change in the law of viscosity which allowed 
a relatively greater influence to the semi-diurnal tides would cause a greater change of 
obliquity, and this without much affecting the tidal friction and reaction. Such a law 
seems quite within the bounds of possibility. The special hypothesis, however, of 
elastico-viscosity, used in the previous paper, makes the other way, and allows greater 
influence to the tides of long period than to those of short. This was exemplified where 
it was shown that the tidal reaction might depend principally on the fortnightly tide. 

The whole investigation is based on a theory of tides in which the effects of inertia 
are neglected. Now it will be shown in Part III. of the next paper that the effect 
of inertia will be to make the crest of the tidal spheroid lag more for a given height 
of tide than results from the theory founded on the neglect of inertia. An analysis of 
the effect produced on the present results, by the modification of the theory of tides 
introduced by inertia, is given in the next paper. 

On the whole, we can only say at present that it seems probable that a part of the 
obliquity of the ecliptic may be referred to the causes here considered; but a complete 
discussion of the subject must be deferred to a future occasion, when the secular 
changes in the plane of the lunar orbit will be treated. 

The question of the obliquity is now set on one side, and it is supposed that when 
the moon has reached the critical point (where the month is twice the day) the 
obliquity to the plane of the lunar orbit was zero. In the more remote past the 
obliquity had no tendency to alter, except under the influence of certain nutations, 
which are referred to at the end of Section 17. 

The manner in which the moon's periodic time approximates to the day is an 
inducement to speculate as to the limiting or initial condition from which the earth 
and moon started their course of development* 
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So long as there is any relative motion of the two bodies there must be tidal 
friction, and therefore the moon's period must continue to approach the day. It would 
be a problem of extreme complication to track the changes in detail to their end, and 
fortunately it is not necessary to do so. 

The principle of conservation of moment of momentum, which has been used 
throughout in tracing the parallel changes in the moon and earth, affords the means of 
leaping at once to the conclusion (Section 18). The equation expressive of that principle 
involves the moon's orbital angular velocity and the earth's diurnal rotation as its two 
variables ; and it is only necessary to equate one to the other to obtain an equation, 
which will give the desired information. 

As we are now supposed to be transported back to the initial state, I shall hence- 
forth speak of time in the ordinary way ; there is no longer any convenience in 
speaking of the past as the future, and vice versa. 

The equation above referred to has two solutions, one of which indicates that tidal 
friction has done its work, and the other that it is just about to begin. Of the first I 
shall here say no more, but refer the reader to Section 18. 

The second solution indicates that the moon (considered as an attractive particle) 
moves round the ea,rth as though it were rigidly fixed thereto in 5 hours 36 minutes. 
This is a state of dynamical instability ; for if the month is a little shorter than the day, 
the moon will approach the earth, and ultimately fall into it ; but if the day is a little 
shorter than the month, the moon will continually recede from the earth, and pass 
through the series of changes which were traced backwards. 

Since the earth is a cooling and contracting body, it is likely that its rotation would 
increase, and therefore the dynamical equilibrium would be more likely to break down 
in the latter than the former way. 

The continuous solution of the problem is taken up at the point where the moon 
has receded from the earth so far that her period is twice that of the earth's rotation. 

I have calculated that the heat generated in the interior of the earth in the course 
of the lengthening of the day from 5 hours 36 minutes to 23 hours 56 minutes would 
be sufficient, if applied all at once, to heat the whole earth's mass about 3000° Fahr., 
supposing the earth to have the specific heat of iron (see Section 16). 

A rough calculation shows that the minimum time in which the moon can have 
passed from the state where it had a period of 5 hours 36 minutes to the present state, 
is 54 million years, and this confirms the previous estimates of time. 

This periodic time of the moon corresponds to an interval of only 6,000 miles 
between the earth's surface and the moon's centre. If the earth had been treated as 
heterogeneous, this distance, and with it the common periodic time both of moon and 
earth, would be still further diminished. 

These results point strongly to the conclusion that, if the moon and earth were ever 
molten viscous masses, then they once formed parts of a common mass. 

We are thus led at once to the inquiry as to how and why the planet broke up. 

MDCCCLXXIX. 3 z 
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The conditions of stability of rotating masses of fluid are unfortunately unknown, and 
it is therefore impossible to do more than speculate on the subject. 

The most obvious explanation is similar to that given in Laplace's nebular hypo- 
thesis, namely, that the planet being partly or wholly fluid, contracted, and thus rotated 
faster and faster until the ellipticity became so great that the equilibrium was unstable, 
and then an equatorial ring separated itself, and the ring finally conglomerated into a 
satellite. This theory, however, presents an important difference from the nebular 
hypothesis, in as far as that the ring was not left behind 240,000 miles away from the 
earth, when the planet was a rare gas, but that it was shed only 4,000 or 5,000 miles 
from the present surface of the earth, when the planet was perhaps partly solid and 
partly fluid. 

This view is to some extent confirmed by the ring of Saturn, which would thus be a 
satellite in the course of formation. 

It appears to me, however, that there is a good deal of difficulty in the acceptance 
of this view, when it is considered along with the numerical results of the previous 
investigation. 

At the moment when the ring separated from the planet it must have had the 
same linear velocity as the surface of the planet ; and it appears from Section 22 that 
such a ring would not tend to expand from tidal reaction, unless its density varied 
in different parts. Thus we should hardly expect the distance from the earth of 
the chain of meteorites to have increased much, until it had agglomerated to a con- 
siderable extent. It follows, therefore, that we ought to be able to trace back the 
moon's path, until she was nearly in contact with the earth's surface, and was always 
opposite the same face of the earth. Now this is exactly what has been done in the 
previous investigation. But there is one more condition to be satisfied, namely, that 
the common speed of rotation of the two bodies should be so great that the equilibrium 
of the rotating spheroid should be unstable. Although we do not know what is the 
limiting angular velocity of a rotating spheroid consistent with stability, yet it seems 
improbable that a rotation in a little over 5 hours, with an ellipticity of one-twelfth 
would render the system unstable. 

Now notwithstanding that the data of the problem to be solved are to some extent 
uncertain, and notwithstanding the imperfection of the solution of the problem here 
given, yet it hardly seems likely that better data and a more perfect solution would 
largely affect the result, so as to make the common period, of revolution of the two 
bodies in the initial configuration very much less than 5 hours/'" Moreover we obtain 
no help from the hypothesis that the earth has considerably contracted since the shed- 
ding of the satellite, but rather the reverse ; for it appears from Section 24 that if the 
earth has contracted, then the common period of revolution of the two bodies in the 

* This is illustrated by my paper on " The Secular Effects of Tidal Friction," < Proc. Roy. Soc.,' No. 197, 
1879, where it appears that the " line of momentum " does not cut the " curve of rigidity " at a very small 
angle, so that a small error in the data would not make a very large one in the solution. 
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initial configuration must have been slower, and the moon more distant from the earth. 
This slower revolution would correspond with a smaller ellipticity, and thus the system 
would probably be less nearly unstable. 

The following appears to me at least a possible cause of instability of the spheroid 
when rotating in about 5 hours. Sir William Thomson has shown that a fluid spheroid 
of the same mean density as the earth would perform a complete gravitational oscillation 
in 1 hour 34 minutes. The speed of oscillation varies as the square root of the density, 
hence it follows that a less dense spheroid would oscillate more slowly, and therefore a 
spheroid of the same mean density as the earth, but consisting of a denser nucleus and 
a rarer surface, would probably oscillate in a longer time than 1 hour 34 minutes. It 
seems to be quite possible that two complete gravitational oscillations of the earth in 
its primitive state might occupy 4 or 5 hours. But if this were the case, then the solar 
semi-diurnal tide would have very nearly the same period as the free oscillation of 
the spheroid, and accordingly the solar tides would be of enormous height. 

Does it not then seem possible that, if the rotation were fast enough to bring the 
spheroid into anything near the unstable condition, then the large solar tides might 
rupture the body into two or more parts ? In this case one would conjecture that it 
would not be a ring which would detach itself. 

It seems highly probable that the moon once did rotate more rapidly round her own 
axis than in her orbit, and if she was formed out of the fusion together of a ring of 
meteorites, this rotation would necessarily result. 

In Section 23 it is shown that the tidal friction due to the earth's action on the 
moon must have been enormous, and it must necessarily have soon brought her to 
present the same face constantly to the earth. This explanation was, I believe, first 
given by Helmholtz. In the process, the inclination of her axis to the plane of her 
orbit must have rapidly increased, and then, as she rotated more and more slowly, 
must have slowly diminished again. Her present aspect is thus in strict accordance 
with the results of the purely theoretical investigation. 

It would perhaps be premature to undertake a complete review of the planetary 
system, so as to see how far the ideas here developed accord with it. Although many 
facts which could be adduced seem favourable to their acceptance, I will only refer 
to two. The satellites of Mars appear to me a most remarkable confirmation of these 
views. Their extreme minuteness has prevented them from being subject to any per- 
ceptible tidal reaction, just as the minuteness of the earth compared with the sun has 
prevented the earth's orbit from being perceptibly influenced (see Section 19) ; they thus 
remain as a standing memorial of the primitive periodic time of Mars round his axis. 
Mars, on the other hand, has been subjected to solar tidal friction. This case, however, 
deserves to be submitted to numerical calculation. 

The other case is that of Uranus, and this appears to be somewhat unfavourable to 
the theory ; for on account of the supposed adverse revolution of the satellites, and of 
the high inclinations of their orbits, it is not easy to believe that they could have 
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arisen from a planet which ever rotated about an axis at all nearly perpendicular to 
the ecliptic. 

The system of planets revolving round the sun present so strong a resemblance to the 
systems of satellites revolving round the planets, that we are almost compelled to believe 
that their modes of development have been somewhat alike. But in applying the 
present theory to explain the orbits of the planets, we are met by the great difficulty 
that the tidal reaction due to solar tides in the planet is exceedingly slow in its 
influence ; and not much help is got by supposing the tides in the sun to react on 
the planet. Thus enormous periods of time would have to be postulated for the 
evolution. 

If, however, this theory should be found to explain the greater part of the configura- 
tions of the satellites round the planets, it would hardly be logical to refuse it some 
amount of applicability to the planets. We should then have to suppose that before 
fche birth of the satellites the planets occupied very much larger volumes, and possessed 
much more moment of momentum than they do now. If they did so, we should not 
expect to trace back the positions of the axes of the planets to the state when they 
were perpendicular to the ecliptic, as ought to be the case if the action of the satellites, 
and of the sun after their birth, is alone concerned. 

Whatever may be thought of the theory of the viscosity of the earth, and of the 
large speculations to which it has given rise, the fact remains that nearly all the effects 
which have been attributed to the action of bodily tides would also follow, though 
probably at a somewhat less rapid rate, from the influence of oceanic tides on a rigid 
nucleus. The effect of oceanic tidal friction on the obliquity of the ecliptic has already 
been considered by Mr. Stone, in the only paper on. the subject which I have yet 
seen/'" His argument is based on what I conceive to be an incorrect assumption as to 
the nature of the tidal frictional couple, and he neglects tidal reaction ; he finds that 
the effects would be quite insignificant. This result would, I think, be modified by a 
more satisfactory assumption. 



* Ast. Soc. Monthly Notices, March 8, 1867. 
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